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ABSTRACT 



Fleet Numerical Meteorology and Oceanography Center (FNMOC) forecasts the 
atmospheric environment and weather using several meteorological and oceanographic 
models. These models’ forecasting abilities are verified by comparing the model forecast 
against the observational data and model’s analysis. Currently, some models are verified by 
several inconsistent, maintenance-intense, non-standardized, and hard-to-use model 
verification systems designed for a particular model. Some models are not verified because 
there is no model verification system. 

This thesis demonstrates the concept of a single model verification system for all 
FNMOC models to eliminate the inconsistencies and redundancies. The single model 
verification system standardizes the model verifications and provides the ability to verify those 
models which are currently unverified. The prototype used a GUI and web browsers to 
display the model verification statistics. The prototype demonstrates that convenient access 
to the model verification statistics could aid FNMOC users in evaluating the forecast models’ 
performance. 

This thesis identifies and documents the user specified verification requirements for 
several models and implements the most immediate requirements. A complete quantitative 
model verification system for all FNMOC models will be implemented incrementally, as all 
the requirements are identified. 
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I. INTRODUCTION 



Fleet Numerical Meteorology and Oceanography Center (FNMOC) is a U.S. Navy 
organization responsible for preparing and disseminating a wide variety of weather and other 
environmental forecasts. These forecasts support a large and diverse group of civil and 
military users who are located throughout the world. These forecasts are generated by 
meteorological and oceanographic models maintained by FNMOC. The accuracy and 
timeliness of the model forecasts are critical because the model forecasts support operational 
missions. M part of the quality control process, the model forecasts are periodically verified. 
While the verification process and frequency vary by model, generally it is an ongoing 
process. The model verification process is currently accomplished by diverse processes. 
Generally, all verification systems calculate a variety of summary statistics. These summary 
statistics are then examined to identify strengths and weaknesses. The models are then 
modified to improve their forecast performance. These summary statistics are also included 
in the weekly, monthly, quarterly and annual model performance reports. These reports are 
read by many meteorological and oceanographic organizations. 

This research attempts to provide an enhanced model verification system that 
streamlines and standardizes the model verification processes in FNMOC. An enhanced 
model verification system that can be used to verify all FNMOC models would streamline the 
model verification process. Additionally, a verification system which is easily modifiable to 
meet new requirements would provide enhanced capability. Furthermore, a system that is 
easily accessable to the model developers and managements would reduce difficulties in 
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rapidly assessing verification data. 

First, FNMOC needs a common verification system to verify all the model forecasts. 
Ideally this verification system should enable the users to make comparisons between several 
different combinations of verification parameters by controlling the independent variables. 
In other words, a verification system should provide the equivalent comparison statistics for 
the variable of interest. The comparison might be between models for the same geometry, 
between the atmospheric levels, between the forecasting period, etc. For example, if a user 
wants to compare one model’s performance with another model’s performance, the 
verification system should provide the equivalent comparison statistics for the different 
models in the same geometry, same atmospheric level, same forecast periods, etc., thus 
controlling these commonalities and providing only the true comparison statistics of the 
model performances. 

Second, FNMOC needs a flexible verification system that can accommodate 
requirement changes. These changes could be an addition of a new model, change in the 
geometry where the models forecast, addition of new environmental parameters, etc. The 
requirement changes need to be incorporated immediately. Additionally, these changes 
should not require extensive modifications to the applications source code when each change 
is required. 

Third, FNMOC needs an easy-to-use tool to access the verification statistics. It 
should be easy for the model developers and managers to use to generate the needed 
information in graphical format. 
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A. CURRENT MODEL VERIFICATION AT FNMOC 



Currently, there is no single verification system for FNMOC model developers to use 
to access and easily analyze the models’ performance. There are several model verification 
systems for different models developed by individual model developers and model verification 
personnel without adhering to any standards. For example, FNMOC’s global meteorological 
model. Navy Operational Global Atmospheric Prediction System 
(NOGAPS ) , is verified by two verification systems, nogs tat and verobs. Nogstat 
verifies NOGAPS forecast against its analysis and verobs verifies against the observational 
data. FNMOC’s regional meteorological models. Navy Operational Regional 
Atmospheric Predictiori System (NORAPS) and Coupled 
Ocean/Atmosphere Mesoscale Prediction System (COAMPS) , are verified 
against the observational data by verobs . Verobs compares NOGAPS and NORAPS in 

the NORAPS regions - Asia, Continental US (Conus), Europe and Indian Ocean. Verobs 
also provides NORAPS and COAMPS comparisons; however, their geometries do not match 
exactly. The different geometries add undesirable variability in the model comparison. 
Another example is FNMOC’s global wave model, WAM_GLOBAL. It is verified by a 
verification scheme developed by the model developer. There is no known documentation and 
the developer is the only person who has any knowledge about the scheme. There are other 
systems that can only handle a single model or are very difficult to modify. Additionally, there 
are redundancies in the model verification efforts among many of the individual verification 
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systems. Finally, some models do not have a verification system at all. 

This lack of standardization also limits FNMOC’s ability to rapidly respond to user’s 
new requirements. As an operational center responding to a variety of operational users 
ranging from the Navy, Air Force and other civil and government organizations, FNMOC 
fi'equently receives requests to provide modifications to models. The ability to respond is 
limited by the fact that the data to compute the model verification statistics are in various 
formats at various locations. Some data are located in FNMOC's relational database. 
Integrated Stored Information System (ISIS) on the Cray and the Sun 
server. The model forecast data (grid data) and the model analysis are stored in the ISIS 
grid database. The observational data are stored in the ISIS latitude/longitude/time (LLT) 
database in a format different fi’om the grid data. Some data are stored in a private directory 
as text files. All of these different data require different reading schemes. The model 
verification statistics are computed in many different ways since there is no single library to 
compute the statistics. The model statistics files are in different formats - some in text files, 
some in the Fortran formatted files, and in different locations - some on the Cray in the 
operational directory, some in private directories, some on a Sun system. 

Finally, the accessibility to and the display of the data in the current systems do not 
provide the flexibility needed. Accessibility to the verification statistics is very limited in the 
current verification systems. Only users who have in-depth knowledge about the verification 
system can access the model verification data. The display of the model verification data is 
in numerous formats, styles and conventions. Notably, graphical displays with similar 
purpose are prepared using multiple graphical software packages. This results in a wide 
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variety of styles and methods used to display similar model statistics. Clearly, there is a 
strong need for one universal model verification system to verify any FNMOC models and 
standardize the model verification process. 



B. FNMOC MODEL VERIFICATION SYSTEM 



The FNMOC Model Verification System is a single, easy-to-modify and easy-to-use 
model verification system for all FNMOC meteorological and oceanographic models. It was 
developed by modifying the current FNMOC operating verification systems verobs and 
nogs tat. The new model verification system standardizes the model verification in 
FNMOC. It also makes the data and information widely available by leveraging web 
technology. It currently resides on FNMOC’s intranet, but could be placed on FNMOC’s 
internet server to meet the external user’s requirements. 

C. RESEARCH QUESTIONS 

FNMOC model developers need a single, easily modified and easy-to-use model 
verification system that would standardize the model verification. This research attempted to 
address these needs through the following questions. 

Is it feasible to develop a prototype system that demonstrates the concept of one 
universal model verification system to verify all FNMOC models? Will this system assist in 
eliminating redundancies and inconsistencies caused by the present system that uses multiple 
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model verification methods? Is it feasible to implement a common system for those models 
that currently do not have a model verification system? Can this standardization of the model 
verification create a reusable statistics library and consolidate all model statistics in one 
database in a standardized data format? Can this model standardization be created using 
graphics with the same format, style and conventions? 

Second, can a universal system provide the flexibility needed to rapidly and efficiently 
modify the model verification system for all FNMOC models? Will this flexibility meet 
requirements of FNMOC customers? 

Third, can a web-based technology provide the access, and ease of use needed to meet 
the customer and FNMOC personnel needs? 

D. THE OVERVIEW OF THESIS 



Chapter I of this thesis provides an introduction and problem statement. Chapter II 
describes the Model Verification System including the background. Chapter HI provides the 
user requirements. Chapter IV provides the design. Chapter V provides the evaluation and 
Chapter VI provides the recommendations for future research. Appendix A provides the 
values of the SMS environment variables and namelist values for each model using the system. 
Appendix B provides the source code. Appendix C provides the prototype testing data. 

The figures in this thesis use Yourdon’s [Ref 1] graphical modeling notation. A 
process is represented by a rectangular box. A process that is further decomposed is 
represented by a shaded rectangular box. An input or output flow is represented by an arrow. 
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An external entity is represented by a box with a folded upper left hand comer. A decision 
point is represented by a diamond. A data store is represented by a dmm. A library module 
is represented as a rectangular box with double side bars. 

The courier font is used to indicate programs and systems external to the model 
verification system such as ISIS and HTML. The italic is used to indicate the parts of the 
model verification system such as the graphics component. 
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n. FNMOC MODEL VERIFICATION SYSTEM 



The FNMOC Model Verification System is a single model verification system for all 
FNMOC meteorological and oceanographic models. It replaces the currently existing model 
verification systems that are redundant, hard to modify and maintenance intensive. It provides 
an easy-to-use model verification system for the models that do not currently have a 
verification system. It standardizes the model statistics computations by creating a reusable 
statistics library. It standardizes the format and storage of the model statistics by 
consolidating the model statistics in one location using one database management system. It 
also standardizes the display of the model statistics by using the graphics library routines 
created with FNMOC ’s graphics software. This system provides FNMOC model developers 
a easy-to-use tool to aid them in their analysis of model performance. 

It uses web technology to enable the users to make their requests via web browsers 
and receive the resulting graphics as shown in Figure 1. Implementing the GUI via a web 
browser eliminates the requirement that the client system have all the necessary software 
locally, i.e., all the processing is done on the server. 

Context Level Diagram 





user request ^ 






► 


Model 


User 




Verification 




^ graphics 


System 



Figure 1 . Context Level Diagram 
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The model verification system has several components. The operational run (ops run) 
component computes the model statistics twice a day after each model run is completed. The 
statistics computation library (statistics library) component is a reusable library that performs 
the statistics computation for this system as well as other systems (programs) at FNMOC. 
The database for the model statistics (statistics database) component is the single data 
storage for all the model statistics. The graphics component creates the graphs of the model 
statistics. The user interface component receives the user requests, processes them, then 
returns the resulting graphs back to the user via a web browser. The following sections 
describe these components in broad, general terms. More specific descriptions of these 
components are in Chapter IV. 



A. OPERATIONAL RUN 



FNMOC runs the prediction models twice a day. The model verification statistics are 
calculated after each operational run of the models. This component has four major parts. 
First, after the prediction model has successfully completed, the FNMOC operations staff uses 
the Scheduler Monitor Supervisor (SMS) system to invoke a Korn shell job 
script, mverifjob. This operational job script (mverifjob) is run on the same FNMOC Cray 
computer as the model forecasts, model analysis and the observational data. The mverifjob 
also executes the second and third parts of the operational run components, verobs and 
veronal. Additionally, it writes a set of text files that contadn the calculated statistics to the 
appropriate directory and computer. Second, a group of Fortran programs, verobs, verifies 
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the model forecasts against the observational data. Third, a group of Fortran programs, 
veranal, verifies the model forecasts against the model analysis. The model analysis is the 
initial conditions over a set of grid points covering the forecast region. Fourth, a Korn shell 
script, statupd.ksh, adds the computed model statistics to the model statistics database. 



B. STATISTICS LIBRARY 



The statistics library has a collection of reusable general purpose statistics programs. 
These programs are a group of Fortran programs that compute the required statistics for 
model verification. These include commonly used verification statistics such as mean error 
(bias), standard deviation of the error (stdev), and root mean squared error (rms). All the 
programs compute the statistics on the difference between the verifying and predicted sets of 
data. There is also a program to compute the map factors for different earth projections used 
by FNMOC. 

The current FNMOC routines for computing statistics were modified and adapted in 
this component. The inventory of the existing statistics computation routines was obtained 
from the FNMOC Unix Utility Library and the applications under FNMOC configuration 
management (CM). Table 1 shows the existing Fortran subroutines and the statistics 
computed by each of them. 
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Table 1. Inventory ofFNMOC statistics routines 



Routine 


Statistics computed 


ancor.f 


anomaly correlation (anc) for 10 NOGAPS areas 


diflferp.f, dififers.f' 


mean bias, mean, sum, standard deviation (std), root mean square (rms) 


htrms.f 


mean bias, std, rms 


nr ancor.f 


anc 


nrhtrms.f 


mean bias, std, rms 


nrvelrms.f 


mean wind u component, mean wind v component, mean wind speed, rms 
wind speed 


stats.f 


mean bias, rms, std, threat score, probability of detection (pd), false 
alarm, skill score 


velrms.f 


mean wind u component, mean wind v component, mean wind speed, rms 
wind speed for 1 0 NOGAPS areas 



* The statistics computed do not take any projection into account. All the other subroutines 
compute the statistics for the spherical projection. 



These routines are redundant, not in a particularly reusable form, and application 
specific rather than general purpose. For example, the routines ancor . f , htrms . f and 
velrms.f are specific to nogstat. The routines nrancor.f, nrhtrms.f and 
nrvelrms . f are specific to norstat (a no longer active modification of nogstat). 
The routine stats . f is specific to verobs. The routines differs . f and dif f erp . f 
are part ofyet another application Ocean Model Support Program (OMSP). This 
component makes them more reusable and general purpose. 



C. STATISTICS DATABASE 



This component is the data storage for the model verification statistics. It uses 
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FNMOC’s relational database system, Empress. There is a table for each model and a 
generic table structure used to create the table structures for each model’s table. This 
database resides on the database workstation. This component allows easy data access and 
consolidates all the model statistics in a single location. It will have one year’s data on line 
and archive the older data off line. This is to accommodate the model developers’ 
requirements to frequently use the historical data for weekly, monthly, quarterly, seasonal and 
annual analyses and reports. 

D. GRAPHICS 



This component creates the graphs of the model statistics based upon the user’s 
requests. The programs are written with one of FNMOC’s graphics software packages. 
Interactive Data Language (IDL). These programs read the model statistics 
from a text file. The model statistics are retrieved from the database, formatted and written 
to a text file for IDL programs. The IDL program creates a GIF file to display in a Hyper 
Text Markup Language (HTML) page back to the user. 

E. USER INTERFACE 



This component is the gateway to the model verification system via a web browser. 
It gives the users the ability to compose their requests for the model verification statistics they 
want via HTML forms. When the users access the system they are presented with a variety 
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of on-screen choices. These include a list of models, forecast periods (tau), atmospheric 
pressure levels, parameters (air temp, wind speed, etc.), statistics, and verification sources. 
Additionally, they can select the observation types and types of graphs and can also specify 
the time period they want graphically presented. After they make the various selections and 
request the information, their request is processed and the information is returned to their 
screen. Their request is processed dynamically, using a perl script. This is important since 
it eliminates the need to have any pre-created static GIF files. This is an important capability 
since static files are less flexible, use greater storage space, and require a great deal of 
maintenance to keep them current. The perl script parses the user requests, executes the 
script to query the statistics database and retrieve the data, executes the graphics programs 
and returns the graphs to the user via a web browser. 

F. MEETING THE REQUIREMENTS 

Each of the FNMOC Model Verification System components demonstrates the ability 
to eliminate shortfalls in the current system as well as providing additional enhancements. 
First, the operational run component computes the statistics twice daily for continuous 
accumulation of model statistics. It is designed for flexible modification whenever a 
requirement change occurs via SMS variables and the namelist functionality in Fortran 90. 
Additionally, the statistics database standardizes the data format, data storage, data access 
and data location by having all the model statistics in a single database. Also, the statistics 
library standardizes the statistics computation for aU the map projections. The text file format 
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for the intermediate files for the data insertion into the statistics database standardizes how 



the intermediate files need to be written. This is beneficial in the event these files need to be 
read directly or if the model statistics from some other systems need to be added into the 
statistics database. The graphics component creates the standardized graphs of the model 
statistics. The user interface component provides the users an easy-to-use model verification 
system for all FNMOC models. No special training is required to use the model verification 
system if users know how to use a web browser. It is a single point and common way to 
access the model statistics without any knowledge about the system. All these components 
help to standardize the various aspects of the model verification in FNMOC. 

The FNMOC Model Verification System Prototype demonstrates the concept of a 
single, universal model verification system for all FNMOC models. It provides a vehicle for 
a better understanding of the environment and requirements problem being addressed. It 
demonstrates what is actually feasible with the existing technology, and where the technical 
weak spots still exist at FNMOC. Additionally, it is an effective way to ensure the 
requirements accurately reflect the user’s real needs. The prototype demonstrates to the users 
what is functionally feasible and provides an analysis test bed and vehicle to validate and 
evolve the system requirements. [Ref 17] 
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m. REQUIREMENTS 



A. GENERAL REQUIREMENTS 

Generally, the model forecasts are verified against the observed meteorological values 
or the initial conditions over a set of grid points covering the forecast region called the 
model’s analysis. “V. Bjerknes in 1904 recognized that forecasting is fundamentally an 
initial-value problem in mathematical physics and, moreover, that the basic system of 
equations to be solved was already known, at least in general form.” [Ref 10] This research 
addresses a part of the general requirements, verification against the observational values, but 
future efforts could expand to address the entire general requirements in the future. 

B. FNMOC USER REQUIREMENTS 

FNMOC user requirements were obtained through a user survey, review of the current 
model performance summary reports [Ref 13], and a meeting with the model team leaders 
and model verification personnel. The initial user requirements indicated that the models need 
to be verified against the model’s analysis and/or observational data which matches with the 
general requirements. The verification against the observational data appears to be more 
widely and frequently used at FNMOC. Therefore, this research concentrated on the 
verification against observational data initially and will add the verification against the model’s 
analysis in the future. Table 2 below shows the results of the user requirements analysis for 
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each model. The requirements analysis also showed that bias, standard deviation (stdev) and 
root mean square (rms) are used more frequently than any other statistical measures. Users 
indicated that the various model statistics should be stored on-line and immediately accessible 
for up to one year and off-line for longer period. The users indicated the most widely used 
graphics to display these statistics are scatter plots, time-series plots and height-series plots. 
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Table 2. List of user requirements for each models 
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C. ASSUMPTIONS 



Several assumptions were made during this research to focus on meeting the most 
urgent requirements. First, this research concentrates only on verifying the model forecast 
against the observational data since more users indicated they would like to verify against the 
actual data. The assumption here is that the results and methodology can be generalized to 
the verification versus the model analysis as well. Second, the system is purposely designed 
to be as independent from FNMOC’s relational database management system (RDBMS) 
Empress as possible by using the text files between the DBMS and the system in the event 
the DBMS is changed to some other RDBMS. The assumptions is that these text files can 
be interfaced with any RDBMS. Third, WAM_GLOBAL is used as the representative ocean 
model in the initial system rather than all the ocean models. Other ocean models as well as 
the meteorological models will be added later. This assumes that the results and methodology 
can be generalized to other models as well. Fourth, the system only calculates the bias, stdev 
and rms initially and will add the other statistics later. The asumption is that other statisics 
can be added. 

The system will evolve and mature as the users’ inputs from the user prototype 
evaluations are evaluated. More general verifications will be implemented as the full system 
is implemented to meet all the requirements identified during the requirements analysis, 
namely, the verification against the model analysis, the addition of other oceanographic and 
meteorological models, the addition of statistics (anomedy correlation, probability of 
detection, threat scores, etc.), consideration of the map factors in the verification against the 
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model analysis and addition of other graphics for different combinations. 



D. LIMITATIONS OF RESEARCH 



This research demonstrates the concept of a universal model verification system for 
all FNMOC models by use of a prototype system. This prototype is implemented to meet the 
minimum requirements, but provides the capability to expand to meet the additional 
requirements as well as future requirements. The limitations in the system are listed in the 
following paragraphs. 

The model statistics is computed for different observation types separately rather than 
combined. The map factors are not considered in the verification against the observations. 
In other words, when the model forecasts are verified with respect to the observational data, 
all the latitude/longitude points are considered equally weighted. This system creates only the 
time-series plots to narrow the scope of the graphics component. Other types of graphics can 
be added in the future. 

The prototype is bound to the current FNMOC operating systems and software since 
it should use the FNMOC environment. Therefore, this system may not be as portable as it 
could be when there is a significant change within the FNMOC environment, such as a 
different operating system, new DBMS, new graphics software, new standard shell, etc. 
Some of the lower level modules that interface to the ISIS latitude/longitude/time (LLT) 
data have unavoidable coupling to ISIS due to the different observational data type 
structures. 
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rV. DESIGN 



The initial prototype system design has narrow scope to concentrate on meeting the 
immediate and minimum requirements. Additional requirements will be added in future 
versions to test the system’s flexibihty. The immediate requirements are veriJBcation of the 
NOGAPS , NORAPS , COAMPS and WAM_GLOBAL models after the model runs for daily, 
weekly, monthly and quarterly reports. The prototype will present bias, stdev and rms on a 
few specified parameters (e.g., air temp and wave height). The WAM_GLOBAL model is used 
as a representative ocean model in the prototype. 

The prototype of the initial system addresses these requirements to prove the design 
and implementation concept. Additionally, this will demonstrate the feasibility of using a 
single universal model verification system for FNMOC. “The system need not have a finely 
tuned prototype before it is implemented. In fact, one merely asks of the prototype that it 
contains a flexible set of hardware and software tools for continuous redesign, and have 
access to real-time data (with standardized formats).” [Ref 15] The prototype has a GUI to 
allow the users to view the model statistics in graphical form on a web browser. 

The prototype is developed and implemented in FNMOC’s operational environment 
to ensure that it is useful in that environment. “The important thing is that the development 
takes place primarily in the operational office, with the direct involvement of the people who 
are going to use it.” [Ref 15] 
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A. 



OVERALL STRUCTURAL DESIGN 



Figure 2 shows the overall structural design of the FNMOC model verification system. 
It shows the five main components described in Chapter II and their relationships. 
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Figure 2. Level 1 Diagram 



The operational run component computes the model verification statistics twice a day 
by executing either verobs or veranal. Verobs reads the model forecasts and analysis fi’om 
the ISIS grid database and observational data from the ISIS LLT database. Verobs and 
veranal use the statistics computation routines in the reusable statistics library component 
to perform the computation. They write the statistics in the text files in the standardized 
format (one for each model) for transfer to the FNMOC database development workstation 
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where the statistics database component resides. The operational run component then 
inserts the model statistics into the database. The transfer and the database population is not 
done twice a day by FNMOC operations currently as part of the ops run because operations 
is not one of the owners of the statistics database component. Currently, the transfer and 
database populations are done by FNMOC staff about once a day during the week days. This 
is because the statistics database is not part of ISIS currently, and the FNMOC ops run 
does not write to any non-ISIS databases. The statistics database should become a part of 
ISIS after the prototype evaluation. The user interface component presents HTML forms 
to users which they use to make their requests. This component processes the user requests 
by retrieving the model statistics from the statistics database and formats the data for the 
graphics component. The grcq)hics component then creates the graphics to send to the users. 
The following sections describe the detailed design for each module in each of the five 
components. The source code for all the modules is in the Appendbc B. 

B. DESIGN OF OPERATIONAL RUN (OPS RUN) 

The ops run has two korn shell scripts, mverif.job and statupdksh. It also has 
Fortran programs verobs.j90 and veranal.f90 (in the future) along with a Fortran header file 
v_data.h. mverif.job reads the SMS variables and executes the appropriate program based 
on the SMS variable VERIF SOURCE value. FNMOC operations executes this job script 
to compute the model verification statistics twice a day after each model run is completed. 
The model statistics are written to a text file for each model, geometry and date-time-group 
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(dtg). Statupdksh transfers the model statistics text files fi'om one computer to another 
workstation, formats the text file and populates the statistics database using Empress 
Standard Query Language (SQL). Figure 3 shows the design of ops run 
graphically. 
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model-verif.job (ops run) 




Figure 3. Ops Run Diagram 
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1 . 



Design of mverif.job (Cray) 



This system needs to be flexible to incorporate new requirements and changes to the 
existing requirements. The use of SMS environment variables in mverif.job is one of the 
methods this system uses to assure this flexibility. FNMOC uses SMS to run the operational 
runs (ops runs) and uses the SMS environment variables extensively to control the ops runs. 
The SMS variables are set before the mverif.job is executed and passed into the job. Figure 
4 shows the list of the SMS variables and their description. By changing the value of these 
SMS variables, mverif.job is executed to verify all the FNMOC models without changing the 
job itself The values for these SMS variables were set based on the user requirements 
analysis shown in Table 1. The complete SMS value list for each model is included in the 
Appendix A. 



CRDATE - watch date-time-group (e g., 19980131 12) 

ISIS_TABLE - ISIS table name, coincides with the model name (e.g., NOGAPS , 
WAM_GLOBAL) 

GEOMNMl - geometry name defined in ISIS (e.g., global_360xl81) 

GEOMNM2 - second geometry name 

TAUI - beginning forecast time 

TAUE - ending forecast time 

TAUINC - increment of forecast time 

VERIF_SOlJRCE - either ‘obs’ or ‘anal’ 

OPSBIN - directory in which the binaries for the operational runs reside 
ISIS_INIT - ISIS initialization script name 

PROGBIN - directory in which the binary for this system resides, may be same as 
OPSBIN 



Figure 4. SMS Variables and Description 
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mverif.job is a korn shell script that FNMOC operations execute to compute the 
model verification statistics twice a day. This job performs various operational job required 
tasks such as setting up the operational job environment and executing the correct ISIS 
initialization script. It reads the required SMS variables and exports them to the subsequent 
shells. It creates the appropriate sub-directory with the month (mmm, e.g., mar) and year 
(yyyy, e.g., 1998) from CRDATE to write the model verification statistics. There are 
separate sub-directories for each month and year. It copies the model’s namelist files to the 
$TMP directory and executes the appropriate executable file based on the variable 
VERIFSOURCE’s value as shown in Figure 3. If the value of VERIFSOURCE is ‘obs,’ 
verobs is executed and if the value is ‘anal,’ veranal is executed. It then writes the resulting 
model statistics into the standardized text file in the appropriate data directory 
(/a/ops/etc/dynamic/app/mverifimmmyyyy). It finishes the process by completing the job 
accounting information and writing the joblog file to the operational joblog output 
directory (/a/ops/job/out). 

2. Design of statupd.ksh (Sun) 

This is a korn shell script that inserts the model statistics into the appropriate data 
table in the statistics database, stat db. It is executed after mverif.job is finished. It formats 
the transferred text file using the awk . Statupd.ksh determines if new statistics files exist 
in the data directory. It also determines the database table to populate based on the model 
name in the statistics files. It then inserts the statistics into the appropriate data table in 
stat db. 
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3. Design of Fortran Include File, V DATA.H 

All data definitions and constants used throughout the Fortran 90 programs (verobs, 
veronal) and their subroutines are in the Fortran include fiJe called V_DA TA.H. This provides 
for easier maintenance since the modification is isolated in one location rather than in multiple 
locations when changes in a variable occur. 

4. Design of VEROBS.F90, The Main Verification Against Obs Program 

This is the top level Fortran program that verifies the model forecasts with respect to 

the observational data for various models. Verobs.f90 uses the namelist fianctionality. This 
helps to achieve flexibility in addition to the usage of the SMS variables. Like the SMS 
variables, namelist values are set external to the program. Each model can have a namelist 
file for the surface level and a namelist file for the upper levels. The namelist files are 
explained in more detail in subsection a below. Figure 5 shows the structural design of 
verobs module. 
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verobs 




Figure 5. Verobs Design 



Verobs.f90 reads the SMS environment variables described in the ops run. It 
determines if the geometry in GEOMNMl is defined in the ISIS. If it is not defined, the 
geometry information is added. Once the geometry is defined, it determines the earth area 
to read the observations for the model verification. It then reads the model forecast fi'om the 
ISIS grid database and observations fi'om the ISIS LLT database. It interpolates the model 
forecast to the observational data points to compute the bias, stdev and rms using the 



33 



statistics library component. It writes the statistics to a text file. The text filename is 
determined by the SMS variable values for ISIS_TABLE, GEOMNM2 and CRD ATE. The 
modules (Fortran subroutines) called by verobs are described in the following sections starting 
in section b. 

a. Namelist file format for each model 

Each model can have two namelist files, one for the upper environmental levels 
and one for the surface level. The second namelist file is needed because the surface level is 
different for different parameters and level types while the upper levels are common levels for 
all the parameters. For example, the surface level for air temperature is 2.0 meters, while the 
surface level for wind speed is 19.5 meters. Figure 6 shows the list of the namelist items and 
their descriptions. The complete namelist files for each model are included in Appendix A 
along with Table 3. Table 3 shows the namelist items specified by the users for each model. 
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parm - the list of model parameters to verify 

dsetname - the data set name for the model forecast parameter in ISIS grid data 

obs_parm - the corresponding parameter in the observational data in ISIS LLT data 

odsetname - the data set name for the observational parameter in the LLT data 

units - unit of the given parm 

Ivltype - level type of the given parm 

obstype - the observational type to read from LLT data 

level - vertical level to read 

stat - type of statistic to compute (bias, std, rms) 

tval - threshold value for the threat score, probability of detection, false alarm, and 
skill score; it will be an optional item since it is not always be used. 

Figure 6. Namelist Variables and Description 



For each value of parm, there must be a corresponding value for dsetname, 
obs_parm, odsetname, units, Ivltype and obstype. For example, if the parm value is 
‘air_temp,’ there must be values for the dsetname which can be ‘fcst_ops,’ the obs_parm 
which should be ‘air_temp,’ the odsetname which can be ‘frimoc,’ the units which should be 
‘deg_K,’ the Ivltype which can be ‘isbr_lvl’ and the obstype which can be ‘raob_qc.’ 
b. boundary 

This Fortran subroutine determines the geographical area from which to 
extract observations to use in the verification. It computes the minimum and maximum 
latitude and longitude for a given geometry’s area box. These minimum and maximum 
latitude and longitude are required for reading the observational data from ISIS lit database. 
It fills the single dimension arrays, x and y, with the values representing the left, right, bottom 
and top boundaries. It converts x/y values to latitude/longitude values by calling FNMOC 
utility, vxyll. It then finds the minimum and maximum latitude and longitude. 
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c. f2ob 

This Fortran subroutine interpolates the model forecast fields to the 
observation locations. It obtains the number of rows and columns for the geometry defined 
by GEOMNM2 by the ISIS utility getgeom and finds the maximum array dimensions for 
X and y arrays. It converts observations latitude/longitude points to x/y values by calling 
FNMOC utility, vllxy. It then interpolates the model forecast to the observation points by 
calling FNMOC utility, f intip. 

d. intgeom 

This Fortran subroutine interpolates one geometry to another geometry. This 
interpolation is necessary when a model needs to be verified in the geometry other than the 
geometry in which it is stored. For example, the global models such as NOGAPS are 
compared to the regional models such as NORAPS_ASIA in the same Asia region. It gets 
the information on the two geometries defined by GEOMNMl and GEOMNM2 by calling 
the ISIS utilities ggrd and getgeom. It then interpolates the first geometry to the second 
geometry points by calling the FNMOC utility chgeom. 

e lltread 

This Fortran subroutine calls the appropriate LLT read subroutine based upon 
the observation type. This subroutine acts as the middle man between the main program 
verobs and the lit read subroutines for each observation type. There are separate LLT read 
modules for each observation types because the include files and the data structures in ISIS 
are different for different observation type. For example, the observation type 'raob_qc’ is 
read by the raob_qc_read subroutine. The subroutine raob_qc_read uses two ISIS includes 
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files, conraion . inc and RAOB_QC . H. It has to use the data structure TYPE (raob_qc_int) 
and TYPE (raob_qc), and use the ‘raob_qc’ sequence type in the call to the ISIS LET read 
utility Ird. Currently, verobs reads only the observation types, raob_qc, sfc_lnd, sfc_ship, 
sfc_ship_met_qc and alty using raob_qc_read, sfcjndjread, sfc_ship_read, 
sfc_ship_met_qc_read and altyjead, respectively. Other observation types will be added 
as more models and parameters are added to the model verification system. 

C. DESIGN OF STAT LIB 

The requirements analysis indicated the following statistics should be included in the 
initial prototype; Mean error (bias), standard deviations (stdev), root mean square (rms), 
anomaly correlation (anc) for most of the model parameters and threat score, probability of 
detection, false alarm, and skill score for selected model parameters. This research focusd 
on just the first three statistics since they are most often used. The other statistics will be 
added in future versions. Some of the existing routines will be modified to make them more 
general purpose and added as part of the statistics library. The design of each statistics 
module is described in the following sections. 

1. Map Factors 

The basic formula for the statistics are the same, but the statistics routines need to 
take the map projeaions into account to compensate for differences in earth area at different 
latitudes on the earth. This is achieved by applying a map correction factor at each grid point. 
Then the routines calculate the area using the new x- and y- coordinates. The general 
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mathematical formula for the statistics is 



J surface s(A,(p)6a / surface area = (Xs„ WJ / £ W„ 

where s=value at a given i/j grid point, X=longitude, (p=latitude, W„ are weighting factors. 

There are two distinct cases to consider in the weighting fectors or map factors. First, 
the observations are all independent. Therefore, each observation has equal weight, and the 
statistics formula becomes a simple average formula. Second, the grid points in which they 
lie must be weighted according to the relative size of the physical area related to the map 
factors for the relevant projection. In general, 

W„=(AXijAyij)„ 

where Axjj = h,jAA^j, Ayj = ^ A(j^ , ^ and |;i depend on the map projection and AX is the 
difference between two longitude, and A(p is the difference between two latitude of the box 
the s value lies. 

According to FNMOC’s model meta data database table, 
ops_meta_grid_db . grid_reg_geom, FNMOC’s models use polar stereo, spherical, 
lambert, and mercator projections. FNMOC verifies NOGAPS in spherical projection against 
its analysis with the map correction factor currently. Some models are verified without the 
map projection consideration against the observational data with the assumption that each 
observation has the equal weighting factor of 1 . The map factors for each map projection 
used at FNMOC are as follows; [Ref 10, Ref 12] 

Spherical: 

K = COS<Pij 

h,=l 
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Polar Stereo: 

K = hy = O^Xl+siiKpij) 



Mercator: 

= hy = coscpij / cos(po 

where (Po is the latitude at which the projection is “true” and can be obtained from FNMOC 
database attribute geom_parm_l in degrees. 

Lambert: 

h^ = hy = (cos(pij/cos(p,)‘’’^ [(l+sin(Pi) / (l+suKpg)]^ 

where K = In(cos(Pi/cos(p 2 ) ^ ln[tan(-jr/4 - (Pi/2) / tan(7t/4 - (P 2 / 2 )] 

(Pj, (P 2 are standard latitudes of the projection; their values can be obtained from FNMOC 
database attributes geom_parm_l and geom_parm_2 in degrees. 

2. Bias (Mean Error) 

“Error is the simple difference of forecast minus verifying analysis or observation. The 
difference (error) field provides a quick look at a model’s forecast performance or bias. Bias 
or tendency describes whether a synoptic field or feature is under or over-forecast.” [Ref 13] 
The advantage of the simple difference fields is that they are easy to compute and understand. 
They provide a quick look at the model forecast performance. [Ref 13] The formula used 
for bias is 

E(F„ -OJW„/EW, where F is forecast, O is the observation. 

3. Standard Deviations (stdev) 

“Standard Deviation (stdev) is a measure of the scatter or variabihty about the mean 
in a series of observations. Standard Deviation is the positive square root of the variance.” 
[Ref 13] The formula used for stdev is 

sqit [(E((F. - / EWJ - (E((F„ - OJWJ / EWJ^]. 
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4. 



Root Mean Square (rms or rmse) 



“Root Mean Square Error (RMSE) is defined as the positive square root of the mean 
square error (MSE). MSE is the mean square of any residual. RMSE is the also called the 
standard error of estimate.” [Ref 13] “The RMSE is a quadratic score that gives the average 
magnitude of the errors. This statistics gives more weight to large errors than to small errors 
in the average, and is useful when large errors are undesirable.” [Ref 16] The formula used 
for rms is 

sqrt((E(F,-OJ^WJ/EWJ. 

D. DESIGN OF STATISTICS DATABASE 



The statistics database component uses FNMOC’s RDBMS, Empress, to 
consolidate all the model statistics in one location. The database has a table for each model. 
The data in the tables are based upon the model developers’ recommendations. A 
consolidated database will speed up the data retrieval and insertion. It will also accommodate 
modifications if the table structure needs to be modified for a particular model in the future. 
Figure 7 shows the SQL command used to create the generic table model stats. 
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CREATE TABLE MODEL STATS 



(verif_date 


character(10,l) 


not null. 


sample_size 


integer 


not null. 


parm_name 


character(32,l) 


not null. 


unit_name 


character(32,l) 


not null. 


geom_name 


char acter(3 2,1) 


not null. 


lvl_type 


character(24,l) 


not null. 


levell 


float(2) 


not null. 


tau 


integer 


not null. 


stat_type 


character(16,l) 


not null. 


stat_value 


float(2) 


not null. 


verif_source 


character(8,l) 


not null. 


obs_type 


character(24, 1) 


not null) 



Indices: NORMAL (2, 15) parm_id ON (parm_name) 
NORMAL (2, 15) geom_id ON (geom_name) 
NORMAL (2, 15) level_l_id ON (level_l) 
NORMAL (2, 15) taujd ON (tau) 

NORMAL (2, 15) stat_type_id ON (stat_type) 
NORMAL (2, 15) obs_type_id ON (obs_type) 



Figure 7. SQL used to create the generic table 



The basic table attributes are shown in Figure 8. This structure was saved in a file 
stattblhy the Empress command ‘displc^ modeljstats all dump into stattbl;. ’ The tables for 
other models were created using the structure stored in the file stattbl by the Empress 
command ‘create nogaps from stattbl; 

The use of the indices improved the retrieval performance fi’om greater than 3 
minutes to less than 5 seconds. However, the insertion performance was degraded especially 
as the number of records increased. The tradeoff between these two scenarios favored the 
increase in the retrieval performance. The compromise may be to drop the indices, perform 
the insertion, then rebuild the indices every time data are updated according to a local 
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database expert, but this compromise has not been tried for this research. 




E. DESIGN OF GRAPHICS 

The graphics component consists of one IDL program currently in the prototype. 
This IDL program plots the bias, stdev and rms for a single model, single geometry, single 
parameter, single forecast period, single atmospheric pressure level and single observation 
type on one graph. More IDL programs will be added in future enhancements to create the 
graphs for other combinations. The IDL program creates a GIF file for display in web 
browsers. Figure 9 shows the structural design of the graphics component. 
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Figure 9. Graphics Design 



The IDL program reads the data from a text file into the IDL array. It reads various 
environment variables to use within the program. It builds sub-arrays based upon the 
statistics type, e.g., bias, stdev and rms. Future IDL programs will build different sub-arrays 
based on the models, geometries, forecast periods, atmospheric pressure levels, or 
observation types. It formats FNMOC’s 10-digit date-time group (1998032000) into a more 
meaningful date-time format (OOZ 20Mar 98) to use on the x-axis label. It plots the bias and 
rms with different symbols and colors, and stdev as the shaded areas above and below the bias 
as requested by the users. Future IDL programs will create the scatter plots and other plots 
as the requirements change. An example of the current graphs is shown in Figure 13 in the 
next section. 
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F. DESIGN OF USER INTERFACE 



Figure 10 shows the structural design of the user interface component. This 
component ties together the statistics database and graphics components. It has a GUI and 
Common Gateway Interface (CGI ) to give the users access to the model statistics. 
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Figure 10. User Interface Design 



The GUI portion has several html forms to interface with the users. The first 
model verification system web page shown in Figure 1 1 is created by the file \ndex.html. It 
is shown using a Netscape web browser. 
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Figure 1 1 . Home Page 



When the user picks a model by clicking on a selection button, the model page is 
displayed. An example model page is shown in Figure 12 for the WAM_GLOBAL model. This 
page is created by the file wam.html. The model pages have the default values already 
selected, but the users can change the values by simply clicking other choices. The user 
would then click on the ‘Submit the query’ button when (s)he is finished composing the 
request or the ‘Cancel’ button at any time. 
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Figure 12. WAM_GLOBAL Statistics Page 

Once the user submits the query, the CGI portion of this component takes over. 
The CGI portion has several programs to satisfy the user requests dynamically. Since this 
research was implemented in the current operational environment, several accommodations 
had to be made. The FNMOC intranet web server is on a Sun workstation called ‘devul.’ 
As discussed earlier in the statistics database component, the stat_db is on a database 
development workstation called ‘div60-3.’ The model statistics must be retrieved from div60- 
3 and transferred to devul to work in the current configuration. The GIF images are also 
created on div60-3 because it has better performance than devul. 

When the user’s request is submitted, a CGI perl script in the cgi - bin is executed 
to processes the user request. This CGI script parses the user request and starts the 



47 



processing of the user request via runj ob on div60-3. This run j ob involves data retrieval 
from stat db by Empress SQL. It then formats the retrieved data by awk for IDL to read. 
It executes the IDL program described in the previous section which creates a GIF file. It 
then transfers the GIF file to the web server via f tpbatch. Once the GIF file transfers 
successfully, the GIF image is displayed on the web browser. 

Figure 13 shows the result of the query from Figure 12 for WAM_GLOBAL. At this 
point, the user can print the image, save the image in different formats supported by the web 
browser (e.g. GIF, postscript, etc.) or just view the image. 




Figure 13. WAM_GLOBAL Statistics Graph 
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V. EVALUATION 



The evaluation of the model verification system included alpha testing and beta 
testing. The alpha testing was performed by the developer by comparing the computed 
statistics to the existing verification system verobs. The beta testing was performed by 
FNMOC operations and FNMOC model developers (users). The operations ran the model 
verification system and verobs in parallel for four weeks. The users used the prototype to 
access the verification statistics via the Netscape web browser. The addition of a new 
model, COAMPS_SOUTHWEST_ASIA, demonstrated the ease of extending the model 
verification system to include other models. 

A. ALPHA TESTING 



Each unit (program and subroutines) of the ops run components was unit tested using 
the Fortran debugger, TotalView, which enabled the developer to test the binaries 
interactively. The debugger allows the tester to assign values to variables, print values of 
variables, and step through the statements in the calling program and called subroutines. This 
testing was useful in identifying code-level bugs, e.g., the parameter orders in the actual and 
formal parameter lists of the subroutine call to boundary in the main program verobs did not 
match. 

The statistics computed by the statistics library subroutines were compared to the 
statistics computed using Microsot t Excel. A set of numbers from the model forecast 
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and a set of numbers from the observational data were assigned to two arrays. The statistics 
bias, stdev, and rms were computed by the statistics library subroutines and Excel. These 
two sets of statistics matched as shown in Appendix C. The test program for the statistics 
library subroutines, static st, is included in Appendix C. 

The developer compared the model forecast and observational data at the first 15 
latitude/longitude points of the test run of the model verification system to the first 15 
latitude/longitude points of the operational verobs to ensure they matched point-by-point 
as well as the overall computed statistics. An error in the loop index in the statistics library 
routines was discovered and corrected during the alpha testing. 

The test runs included verification of each model in its own area and the verification 
of the global model in the regional model areas. The purpose was to test the interpolation 
involved when a model was verified in an area other than its own. The global model NOGAPS 
was verified in the regional model NORAPS areas (Asia, Continental US, Europe, Indian 
Ocean) and in the COAMPS area (Europe). This allowed the comparison of multiple models 
in the same geometry. A change in one of the COAMPS geometries was discovered during 
the alpha testing. The COAMPS model developer modified the geometry for the model, but 
the verification systems were not aware of this change. The change in the geometry resulted 
in the job running without computing any statistics. This change was implemented in the 
model verification system by changing a SMS variable. 

The statistics database component was tested by inserting data and performing 
queries against the database. The time required to return data from the initial queries was 
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much longer than anticipated. Modifications to the database tables were made which 
improved query performance significantly. 

A set of data from stat_db was saved in a text file to be the test data for the IDL 
program in the graphics component. The graphics component was tested by printing out the 
data, printing out the subarrays aeated from the data and examining the graphs of the model 
statistics. The user interface component was tested using the Netscape web browser. The 
cgi scripts were installed on the FNMOC intranet web server cgi-bin because the current 
web server configuration does not allow the execution of the cgi scripts in cgi-bin 
directory of individual users/developers. The FNMOC intranet webmaster installed the cgi 
scripts in the web server cgi-bin directory. 

B. BETA TESTING 



The purpose of the beta testing was to compare the model verification system to the 
existing verification process and to test the user interface of the prototype. The beta testing 
was performed by operations and the users at FNMOC. The operational beta testing started 
on 18 March 1998 by FNMOC operations to verify NOGAPS, NORAPS, COAMPS and 
WAM_GLOBAL models in their own areas. The verification of NOGAPS in the NORAPS and 
COAMPS areas was added on 25 March 1998. The model verification system was comparable 
to the existing verification system verobs, but it was much easier to add the verification of 
NOGAPS in the NORAPS and COAMPS areas. 
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/ 

An important part of the beta test was the evaluation of the prototypes ease of use. 
Generally, users indicated that the system was very easy to use. Another area examined was 
how easy it was to access the model statistics. Again, users evaluated this area positively. 
Another area evaluated was whether the system provided the required displays. For the most 
part this area was positively evaluated. However, the users wanted more types of graphs such 
as multiple models on a single chart, wanted the number of observations displayed, and better 
performance, i.e., faster return of the images. Based on the Beta testing results, the prototype 
meets the easy-to-use model verification system requirements. 



C. MODIFICATION 



A new model, COAMPS_SOUTHWEST_ASIA, was added to the prototype after the 
initial product implementation. This provided an opportunity to evaluate how easy it would 
be to modify the verification system. It was very easy to add the additional model to the 
model verification system and no design changes were required. Adding this new model 
required only a few modifications. The operations executes the job script mverif.job with 



appropriate SMS variables 


for the 


model 


(See 


Appendbc 


A) 


after 


COAMPS_SOUTHWEST_ASIA 


completes 


its 


run. 


A 


table 


for 



COAMPS_SOUTHWEST_ASIA was created in the stat_db. The html forms were updated 
to reflect the additional model. This demonstrated that the ease of modification requirement 
was satisfied. 
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D. SUMMARY 



This research has resulted in the development and demonstration of an easy-to-use 
model verification system for all FNMOC models. Furthermore, it has demonstrated that 
model verification can be standardized in FNMOC. The model verification standardization 
was achieved by the several components of the model verification system. The reusable 
statistics library standardized the statistics computation, the statistics database standardized 
the data format, data location, data schema, interface to the database, and the graphics 
standardized the graphical display of the model verification statistics. Additionally, it proved 
the concept of a easily modified verification system. Adding a model to the system required 
no design changes and no source code level changes other than adding new modules to handle 
the additional observation types. Finally, the prototype demonstrated that an easy to use 
system could be developed to return data and graphs to the user’s desktop. The use of a GUI 
and a web browser provided the users with an easy-to-use access to the model verification 
statistics. No user training was required. Therefore, the research successfully addressed the 
research questions and has proven the concepts can be implemented in the full model 
verification system. 
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VI. RECOMMENDATIONS AND FUTURE RESEARCH 



This chapter discusses the recommendations for future FNMOC model verification 
systems implementaion based upon the results of this research. It suggests areas of future 
research in fully implementing a single model verification system for FNMOC. 



A. RECOMMENDATIONS 



FNMOC needs to use one verification system to verify all the meterological and 
oceanographical models. This single system will provide many benefits to FNMOC. First a 
single verification system will eliminate the redundancies and inconsistencies of multiple 
verification systems. Additionally, it will eliminate the maintenance-intense verification 
systems needed to maintain data consistency and accuracy. It will also standardize the model 
verification process. An additional benefit is that by using the same verification system the 
model developers will be able to focus more on the model development effort than on the 
verification process. This research proved the feasibility of using a single model verification 
system for all the meteorological and oceanographic models at FNMOC by providing a easily 
modified system. This model verification system demonstrated it can verify all the models and 
standardize the model verifications at FNMOC and became useful to the model 
verification/validation group in FNMOC. This model verification system also provides a 
verification system for the models that do not currently have a verification systems. 

FNMOC can provide users with an easy-to-use tool to access the model verification 
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statistics that meets the users’ needs. Having an easy-to-use tool that is readily assessible via 
the v/eb technology is important for both internal and external users. The internal users could 
be the model developers, model researchers, managers, statisticians, or others in the command 
responsible for FNMOC data accuracy. The external users could be US Navy organizations, 
US Air Force organizations, or other organizations. This research proved that the concept 
of an easy-to-use model verification system via GUI and the web browsers is feasible. It 
provides access to the verification statistics and graphs to anyone who has a web browser. 
Giving users easy access to the model verification statistics makes their efforts to fine tune 
the models and/or evaluate their performance more effective. It is much too difBcult and time 
consuming to access the information in the currently existing model verification systems in 
FNMOC. 

FNMOC needs to document the model verification requirements for all the models. 
This research identified and documented each model’s verification requirements. This will 
be useful for current and future model developers and model users both within FNMOC and 
outside FNMOC. 

FNMOC should use prototyping as much as possible in developing new systems. The 
prototype was a very useful communication tool between the designer/developer of the model 
verification system and the users. It helped to firm up the requirements and identified 
additional requirements that were not identified in the initial requirements analysis, especially 
in the graphics area. 

Finally, FNMOC should build a comprehensive model verification/validation system 
for all the models in FNMOC. The first step in this effort should be building a comprehensive 
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quantitative model verification/validation system that has three subsystems. One subsystem 
should display the model forecast and analysis, another subsystem should display the 
observational data, and a final subsystem should display model verification statistics. This 
research dealt with the third subsystem. Each subsystem should have the same look and feel 
interfaces with a GUI. Addition of subject matter expertise (qualitative aspect) to aid the 
interpretation of the quantative measures will provide a complete model verification/validation 
system for all the models at FNMOC. 

B. FUTURE RESEARCH 



There are many potential extensions for the model verification system to implement 
all of the user requirements identified during the requirements analysis. 

1. Ops Run 

First, the new verobs should be completed by adding more observation types. This 
means adding the Fortran subroutines to read the additional observation types from the ISIS 
LLT database. 

Second, more models should be added to the model verification system. There are 
several additional ocean models to verify, and these ocean models were prioritized for their 
implementation order. The observation types MCSST, Buoy and Bathy need to be read by 
the subroutines buoy_read, bthy_read, respectively, for the OTIS_GLOBAL, 

OTIS_W_ATL and OTIS_W_PAC models. The next models to implement are 
TOPS_GLOBAL, TOPS_W_ATL and TOPS_W_PAC using the same observation types as 
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OTIS models. The users want to compare the OTIS and TOPS models and determine 
whether FNMOC needs to run both models. The observation type SSMI_EDI_ICE needs 
to be read by the subroutine ssmi_edi_ice_read for the PIPS_N_HEM model. 

Third, the verification against the model analysis, veranal should be added. This 
involves several steps. Complete the design, using the existing verification system nogstat 
if feasible. Expand the statistics database and user interface components to reflect the 
veranal statistics. 

2. Statistics Library 

The map factor subroutine in the statistics library component should be revisited and 
completed for all the projections used at FNMOC as described in the design chapter because 
the map factors will be used in veranal. Other statistics computation modules should be 
added to compute the additional statistics such as anomaly correlation, probability of 
detection, threat score, etc. 

3. Statistics Database 

The stat _db should be incorporated into FNMOC’ s database and FNMOC database 
adrninistrators (DBA) should be responsible for administering the stat_db. This will allow the 
operations to write directly into the database as the statistics are computed twice a day. This 
will eliminate the current method of one person transferring and populating the database 
manually. This will also make the latest statistics available to the users, 

4. Graphics 

The graphics component should be expanded by adding more types of graphs. There 
are many different possible graphs based on the combination of models, forecast periods. 
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levels, observation type and statistics. Some of the combinations of interest include multiple 
statistics for a single model, single forecast period, single level, single observation type; 
multiple levels for a single model, single forecast period, single observation type, single 
statistic; multiple forecast periods for a single model, single level, single observation type, 
single statistic; multiple models for a single forecast period, single level, single observation 
type, single statistic. One of the users suggested adding the maps of the geometries used in 
the various FNMOC models and adding the number of observations on the graphs. These are 
very good suggestions and should be added to the system. 

5. User Interface 

The user interface component should be expanded to handle the various types of 
graphs described in the previous paragraph. The html forms should be made to be more 
robust and remove the possibility of user input errors. Currently, there is room for user input 
errors in the html forms in matching the atmospheric levels, observation types and 
parameters because all the levels and observation types (both surface and upper levels) are 
available on the forms. Here is an example scenario to demonstrate the erroneous input. 
Let’s say the user wants to view the graph of NOGAPS for air_temp at the surface level. The 
user clicks air_temp for the parameter, 0 for the level (the correct surface level for air_temp 
is 2.0 meters), raob_qc for the observation type (raob_qc is not a surface level observation 
type). The correct selections would have been air_temp, 2 meters, and sfc_land or 
sfc_ship_met_qc. The users for whom this research is designed are familiar with these facts. 
However, this is a definite shortcoming when this model verification system’s user group 
expands in the future. 
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6. Other Open Issues 

The current hardware configuration in FNMOC has a severe impact on the response 
time in the model verification system. The intranet web server is on an older and slower 
workstation (devul) while the model statistics database is on a different workstation (div60- 
3). This hardware configuration requires a communication (handshake) from devul to div60- 
3. The model verification system then has to use the runj ob script to query the database 
and create a GIF image on div60-3 from devul. Once the GIF image is created, another 
handshake from div60-3 to devul is needed to send the GIF image back to devul from 
div60-3 by FTP. The first handshake takes approximately one minute while the actual data 
query and image generations take approximately seven seconds. It takes another 1 to 2 
seconds to transfer the GI F file. The delays associated with the handshakes, runj ob and 
ftp could be eliminated if the database and web server were on a single, faster workstation. 
The specific hardware requirements for a single server for web, database, and graphics need 
to be identified. Having hardware with the appropriate capability would also ease the 
development and testing process. It would eliminate the need for various routines such as 
runj ob and ftp. This would decrease the time needed to process the user request by 
eliminating the need for the various hardware systems to communicate with each other. 
Finally, this would also help in planning for the future installation of the model verification 
system on the internet for external users. 
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APPENDIX A 



A. SMS VARIABLES 

The following are common SMS variables for all the models; 

ISIS_INlT=/a/ops/isis/db_init/init_ops.ksh 

OPSPATH=/a/ops 

OPSBIN=/a/ops/bin 

CRDATE=$(dtg) 

1. NOGAPS 



TAUI=0 

TAUE=144 

TAUINC=12 

GEOMNMl=global_360xl 8 1 
GEOMNM2=global_3 60x181 
ISIS_TABLE=NOGAPS 
VERIF_SOURCE=obs 

2. NORAPS_ASIA 

TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNM 1 =asia_nest l_appl 
GEOMNM2=asia_nest 1 _appl 
ISIS_TABLE=NORAPS_ASIA 
VERIF_SOURCE=obs 

3. NORAPS_CONUS 

TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNM 1 =conus_nest 1 _appl 
GEOMNM2=conus_nest 1 _appl 
ISIS_TABLE=NORAPS_CONUS 
VERIF SOURCE=obs 
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4. 



NORAPS EUROPE 



TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNM 1 =europe_nest 1 _appl 
GEOMNM2=europe_nest 1 _appl 
ISIS_TABLE=NORAPS_EUROPE 
VERIF_SOURCE=obs 

5. NORAPS_IND_OCN 

TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNM 1 =ind_ocn_nest 1 _appl 
GEOMNM2=ind_ocn_nest 1 _appl 
ISIS_TABLE=NORAPS_IND_OCN 
VERIF_SOURCE=obs 

6. NOGAPS FOR ASIA 

TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNMl=global_360xl81 
GEOMNM2=asia_nest l_appl 
ISIS_TABLE=NOGAPS 
VERIF_SOURCE=obs 

7. NOGAPS FOR CONUS 

TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNMl=global_360xl 81 
GEOMNM2=conus_nestl_appl 
ISIS_TABLE=NOGAPS 
VERIF_SOURCE=obs 

8. NOGAPS FOR EUROPE 
TAUI=0 
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TAUE=48 

TAUINC=12 

GEOMNMl=global_360xl 8 1 
GEOMNM2=europe_nest 1 _appl 
ISIS_TABLE=NOGAPS 
VERIF_SOURCE=obs 

9. NOGAPS FOR EUROPE_NEST2 



TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNMl=global_360xl 8 1 
GEOMNM2=europe_nest2_appl2 
ISIS_TABLE=NOGAPS 
VERIF_SOURCE=obs 

10. NOGAPS FOR EUROPE_NEST3 



TAUI=0 

TAUE=24 

TAUINC=6 

GEOMNMl =global_3 60x181 
GEOMNM2=europe_nest3_appl3 
ISIS_TABLE=NOGAPS 
VERIF_SOURCE=obs 

11. COAMPS_EUROPE 



TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNM 1 =europe_nest2_appl2 
GEOMNM2=europe_nest2_appl2 
ISIS_TABLE=COAMPS_EUROPE 
VERIF_SOlJRCE=obs 

12. COAMPS_EUROPE FOR NEST3 

TAUI=0 

TAUE=24 

TAUINC=6 

GEOMNM 1 =europe_nest3_appI3 
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GEOMNM2=europe_nest3_appI3 

ISIS_TABLE=COAMPS_EUROPE 

VERIF_SOURCE=obs 

13. COAMPS_SOUTHWEST_ASIA FOR NEST2 



TAUI=0 

TAUE=48 

TAUINC=12 

GEOMNM 1 =southwest_asia_nest2_appl 
GEOMNM2=southwest_asia_nest2_appl 
ISIS_TABLE=COAMPS_SOUTHWEST_ASIA 
VERIF_SOURCE=obs 

14. COAMPS_SOUTHWEST_ASIA FOR NESTS 



TAUI=0 

TAUE=24 

TAUINC=6 

GEOMNM 1 =southwest_asia_nest3_appI 
GEOMNM2=southwest_asia_nest3_appI 
ISIS_TABLE=COAMPS_SOUTHWEST_ASIA 
VERIF_SOURCE=obs 

15. WAM_GLOBAL 

TAUI=0 

TAUE=144 

TAUINC=12 

GEOMNMl=global_360xl8 1 
GEOMNM2=global_360xl8 1 
ISIS_TABLE=WAM_GLOBAL 
VERIF SOURCE=obs 
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NAMELIST REQUIREMENTS OBTAINED FROM THE USERS FOR EACH 
MODEL 
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NOGAPS, ALL NORAPS, ALL COAMPS 
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c. 



NON-SURFACE NAMELIST FILES 



1. NOGAPS, all NORAPS and all COAMPS 



&verlst 

parm="air_temp", "geop_ht", "wnd_spd”, 
dsetname="fcst_ops", "fcst_ops", "fcst_ops'', 
obs_parm="air_temp", "geop_ht", "wnd_spd", 
odsetname=''fhmoc'',''fhinoc'V'fhmoc'', 
units="deg_K","gpm'', "m/s", 
lvltype="isbr_lvl","isbr_lvl","isbr_lvl", 
obstype="raob_qc","raob_qc","raob_qc", 
level=1000, 925, 850, 700,500, 400, 300,250, 200, 150,100, 
stats="bias","rms","std", 
tval=35.0 
&end 



D. SURFACE NAMELIST FILES 



1. NOGAPS, all NORAPS and all COAMPS 



&sfclst 

sfc_parm="air_temp","air_temp", "pres","pres", "wnd_spd","wnd_spd", 
sfc_dsetname="fcst_ops","fcst_ops" , "fcst_ops", "fcst_ops", "fcst_ops", "fcst_ops", 
sfc_obs_parm="air_temp","air_temp", "sea_lvl_pres","sea_lvl_pres","wnd_spd", "wnd_spd", 
sfc_odsetname="fhmoc","fhmoc","fhmoc","fhmoc","fhmoc","fhmoc", 
sfc_units="deg_K","deg_K","mb","mb","m/s", "m/s", 
sfc_lvltype="ht_sfc","ht_sfc", "msl","msr',"ht_sfc","ht_sfc", 

sfc_obstype="sfc_lnd","sfc_ship_met_qc","sfc_lnd","sfc_ship_met_qc","sfc_lnd","sfc_ship 

_met_qc", 

sfc_level=2.0,2.0,0.0,0.0, 19.5,19.5, 

sfc_stats="bias","rms","std", 

sfc_tval=35.0 

&end 
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2 . 



WAM GLOBAL 



&verlst 

parm="sig_wav_ht", "peak_wav_per", 

obs_parm="inst_wav_ht_2", "inst_wavjper", 

units="m'', "s", 

lvltype=" surface" , " surface", 

level=0,0, 

stats="bias","rms","std", 
obstype=" sfc_ship " , 
tval=35.0 
&end 
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APPENDIX B SOURCE CODE 



A. OPS RUN 



1. mverif.job 

# QSUB -It 00:30:00 

# QSUB -IT 00:30:00 

# QSUB -Im lOMw 

# QSUB -IM lOMW 

# QSUB -eo 

# QSUB -ko 

# QSUB -o /home/pacek/mverif/etc/mverif.out 

# QSUB -s /bin/ksh 

# QSUB -X 

II — 

# local function to handle the exit code from executing binary 

rNIT_user_exit() { 

case $? in 

0) #Normal exit 

INITJoblog_coniment=MVERIF_COMPLETED_OK 

code=0 

35 

1) #p6cgetenv error in the main program 
IMTJoblog^comment=MVERIF_INFORMATIVE_PFXGETENV_ERROR 
code=0 

>9 

2) #ch2int error in the main program 
IMTJoblog_comment=MVERIF_INFORMATIVE_CH2INT_ERROR 
code=0 

33 

3) #not enough info to continue 

IMTJoblog^comment=MVERIF_INFORMATIVE_NOT_ENOUGH_INFO 

code=0 

33 

4) #output file open error 

INITJoblog_comment=MVERIF_INFORMATIVE_OUTPUT_FILE_OPEN_ERR 

code=0 

33 

5) #error with ggrd call 

INIT Joblog_comment=MVERIF_INF ORMATI VE_GGRD_CALL_ERR 
code=0 

99 

6) #error with getgeom call 
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INIT Joblog_comrnent=MVERIF_INF ORMATI VE_GET GEOM_CALL_ERR 
code=0 

>> 

7) #error with boundary subroutine call 

IMTJoblog_coniment=M\^RIF_INFORMATIVE_BOlJNDARY_CALL_ERR 

code=0 

*) #Abnormal exit, save log file 
IMTJoblog_cominent=MVERIF_F AILED 
code=99 

if [ -f "Stempfile" ]; then 

cp Stempfile $HOME/mverif^etc/o_mverif_err.log 
fi 

esac 

return Scode 

} 



# setup the job environment 

fl^ ^ ■ 

. SOPSBIN/initJob 
set -Sx 

INIT_do_not_notify_operator=0 

IMT_notify_status=EI^ 

INIT_notify_users="pacek" 

. $ISIS_INIT 

^ 

# test to see if env var CRD ATE exist, otherwise take the ops dtg 

- — 

if[-z "SCRDATE"] 
then 

echo "Null CRD ATE, setting to ops dtg." 

CRDATE=$(dtg) 
export CRD ATE 
fi 

... 

# set the required env var and export the SMS and other var 

if- - ■ ■ ■ ■ ■■ 

#PROGBIN=/home/pacek/mverifibin 

PROGBIN=${PROGBIN:-$OPSBIN} 

MVERIF_DIR=$OPSPATH/etc/dynaniic/app/mverif 

#MVERIF_DIR=$HOME/mverif/etc 

echo $MVERIF_DIR 

integer TAUI 

integer TAUE 

integer TAUINC 

export TAUI TAUE TAUINC 

#TEST_DIR=/home/pacek/mverif7test 

TEST_DER=$OPSPATH/etc/static/app/mverif 
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NAMLIST__FILE 1 =$ISIS_T ABLE-verfil 
NAMLIST_FILE2=$ISIS_TABLE-sfcverfil 

print "namlist__filel=" $NAMLIST_FILE1 " namlist__file2=" $NAMLIST_FILE2 
export NAMLIST_FILE1 NAMLIST_FILE2 

II -=r-, = 

# job accounting info 

11 -^-^^-^^^^ --r.- r. - 

ja 



# need to pull out the month and year from CRD ATE 

II ■■ ■ - 

year='echo SCRDATE | cut -cl -4' 
month='echo SCRDATE | cut -c5-6' 
case Smonth in 

01) mon-jan' ;; 

02) mon='feb' ;; 

03) mon='mar’ ;; 

04) mon='apr' ;; 

05) mon='may' ;; 

06) mon-jun' ;; 

07) mon- jul’ ;; 

08) mon- aug' ;; 

09) mon=’sep' ;; 

10) mon-oct' ;; 

11) mon='nov' ;; 

12) mon- dec' ;; 
esac 

il — 

# find the appropriate subdirectory 

- 

subdir=$mon$year 

# need to test for the existance of subdir otherwise the 

# job will crash 

ff -====::.=- 

if [ -d $MVERIF_DIR/$subdir ] 
then 

echo $MVERJF_DIR/$subdir "exists" 
else 

mkdir $MVERIF_DIRy$subdir 
chmod775 $MVERIF_DIR/$subdir 
fi 

II . 

# the output directory 

. ■ -_ = 

OUTDIR=$MVERIF_DIR/$subdir 

echo SISIS^TABLE STAUI STAUE STAUINC 
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II ■ ■ ■■ - 

# determine the current watch and month's stat file name 

II ■ - - = 

cuir_file=$ISIS_TABLE$GEOMNMl SCRDATE 
month_file=$ISIS_TABLE$GEOMNMl$mon$year 
if[[$GEOMNM2 != "]] 
then 

curr_file=$ISIS_TABLE$GEOMNM2$CRDATE 

month_file=$ISIS_TABLE$GEOMNM2$mon$year 

fi 

II . 

# if curr_file exists (ran once) delete so program will not exit 

II . - - - == 

if [ -f $OUTDIR/$cuiT_file ] 
then 

echo "removing stat out file" 
rm $OUTDIR/$curr_file 
fi 

II- — - 

# to use STMPDIR, the scratch area 

II = 

cd STMPDIR 

# - ■ -rz.. - ■= 

# copy the namelist files to the scratch area 

II 

if[[-f$TEST_DIR/$NAMLIST_FILEl ]] 
then 

cp $TEST_DIR/$NAMLIST_FILE1 $NAMLIST_FILE1 
fi 

if [[ -f $TEST_DIR/$NAMLIST_FILE2 ]] 
then 

cp $TEST_DIR/$NAMLIST_FILE2 $NAMLIST_FILE2 
fi 

II - -■=: = ==rT:^r^::r = 

# execute mverif program 

II . 

if [[ $VERIF_SOURCE = 'obs' ]] 
then 

print "calling verobs" 

SPROGBIN/verobs 
elif [[ $ VERIF_SOlJRCE = 'anal' ]] 
then 

print "calling veranal" 

SPROGBIN/veranal 

else 

print "VERIF_SOURCE must be either 'obs' or 'anal'." 
exit 8 
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fi 



H 

# append the current run's stat to the monthly stat file 

II 

#cp $curr_file $OUTDER/$curr_file 
cat $curr_file » $OUTDIR/$month_file 

II rr::rr:^-r:::T .:^-=^:r.z====:= = ^:^,:^:=:Tr:^-T.::=== = ==. 

# close up the job accounting info and clean up 

II - - ___ 

rm -f core 
ja -St 

# START EPILOGUE 

# 

# sees IDENTIFICATION; @(#)mverif.job 1.1 04/24/98 /h/cm/libraiy/mveiif^src/job/mverif.job_v 

# 

# CONFIGURATION HDENTIFICATION: 

# 

# SCRIPT NAME: mverif.job 

# 

# SHELL TYPE: Kom 

# 

# DESCRIPTION: Script that runs MVERIF which computes veriiying stats of 

it the models against the LLT observations or VERANAL which 

it computes verifying stats against the model analysis. 

it 

it COPYRIGHT: (c) 1 998 FLENUMMETOCCEN 

it US. GOVERNMENT DOMAIN 

§ ALL RIGHTS RESERVED 

it 

# CONTRACT NUMBER AND TITLE; NONE 
it 

# REFERENCES; NONE 

# 

# CLASSIFICATION; Unclassified 
it 

it RESTRICTIONS: NONE 
it 

it COMPUTER/OPERATING SYSTEM Cray UNICOS 

# DEPENDENCIES: 

# 

it LIBRARIES OF RESIDENCE: /a/ops/app/mveril/src/job 
# 

it USAGE: qsub mverif.job 
it 

a PARAMETERS; SMS variables needed 
it Name Description 

it 

# ISIS_INIT ISIS init script 
it CRD ATE current run dtg 
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OPSPATH 


ops path 


OPSBIN 


binary directory 


PROGBIN 


test binary directory 


MVERIF_DIR mverif stat output files directory 


GEOMNMl 


geometry name (e.g., conus_nestl_appl) 


GEOMNM2 


geometry name to interpolate to 


ISIS TABLE 


model name (e.g., NORAPS_CONUS) 


TAUI 


starting tau 


TAUE 


finishing tau 


TAunsrc 


tau increment 



# 

# RETURN CODE: 

# 

# FILES: 

# Name Usage Description 

# 

# curr_file IN file that contains the current run's stats 

# month_file IN/OUT file that contains the month's stats 

# 

# DATA BASES: 

# Name Table Usage Description 

# 

# 

# NON-FILE INPUT/OUTPUT: 

# Name Type Usage Description 

# 

# 

# ERROR CONDITIONS: 

# Condition Action 

# 

# data not found stop executing 

# ciuT_file not found cannot append to month_file. 



# sends an email to group 

# ADDITIONAL COMMENTS: NONE 

# 

# MAINTENANCE SECTION 

# 



# EXTERNALS CALLED: 

# Name Description 

# 

# verobs reads ISIS grid data, ISIS LLT data and computes 

# the models’ verifying stats 

# veranal reads model forecast and analysis fi-om ISIS grid 

# and computes the models' verifying stats 

# 

# VARIABLES: 

# Name Description 

# 

# year 4-digit year fi-om SCRDATE 

# month 2-digit month fi-om SCRDATE 

# mon 3-character month (e.g.,jan) 

# NAMELIST_FILE 1 namelist file name for upper levels 
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# NAMELIST_FILE2 namelist file name for surface level 

# 

# METHOD: change directory to the STMPDIR and run the Fortran program 

# verobs or veranal, copy the current run's stat file to the 

# SOPSPATH and append the current run's stat file to the 

# monthly stat file in the SOPSPATH 

# 

# RECORD OF CHANGES: 

# «CHANGE NOTICE» version 1 . 1 (29 Apr 1 998) — Kyongsuk Pace 

# initial submission 

# 

# END EPILOGUE 

2. statupd.ksh 

#!/bin/ksh 

if ■ T ■. ■- = 

# inserts the model stats into the stat_db 

■ := 

DATA_DIR=/home/pacek/data 
DB_DIR=/d/model-stats 
cd SDB^DIR 

//-, - ■= 

§ field separation for Empress 

fl - - - = 

export MS VALSEP- ' 

ff---^ - — - ■ ■ -- 

# is there a new data file? if so, move them 

II— ■ 

if[[ -a$DATA_DIR ]] 
then 

mvSDATA^DIR/*. 

else 

print "There are no data files!" 
fl 

II ;■ :^:= • - - -= 

# determine which model by looking at the files that 

# ends with numeric 0 or 2 

n example: NOGAPSglobal_360xl 81 1998021000 

= 

for OBJ in *[0-9] 
do 

print SOBJ 

awk’/^[0-9]+/ { print$l,$2,$3,$4,$5,$6,$7,$8,$9,$10,$ll,$12 }'\ 

SOBJ > tap 

case SOBJ in 

NOGAPSasia_nestl_appl*) 
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empcmd stat_db "insert into nogaps_asia_nestl_appl jfrom top;";; 
NOGAPSconus_nest 1 _appl*) 

empcmd stat_db "insert into nogaps_conus_nestl_appl from top;";; 
NOGAPSeurope_nest 1 _appl *) 

empcmd stat_db "insert into nogaps_europe_nestl_appl from top;";; 
NOGAPSeurope_nest2_appl2*) 

empcmd stat_db "insert into nogaps_europe_nest2_appl2 from top;";; 
NOGAPSeurope_nest3_appl3 *) 

empcmd stat_db "insert into nogaps_europe_nest3_appl3 from top;";; 
NOGAPSglobal*) 

empcmd stat_db "insert into nogaps_global_360xl81 from top;";; 
NOG APSind_ocn_nest 1 _appl*) 

empcmd stat_db "insert into nogaps_ind_ocn_nestl_appl from top;";; 
NORAPS^ASIA*) 

empcmd stat_db "insert into noraps_asia from top;";; 
NORAPS_CONUS*) 

empcmd stat_db "insert into noraps_conus from top;";; 
NORAPS_EUROPE*) 

empcmd stat_db "insert into noraps_europe from top;";; 
NORAPS_IND_OCN*) 

empcmd stat_db "insert into noraps_ind_ocn from top;";; 

COAMPS_SOUTHWEST_ASIA*) 
empcmd stat_db "insert into coamps_sw_asia from top;";; 
COAMPS^EUROPE*) 

empcmd stat_db "insert into coamps_europe from top;";; 
WAM_GLOBAL*) 

empcmd stat_db "insert into wam_global from top;";; 

*) 

print "no stat table for " $OBJ;; 
esac 

mv $OBJ /home/pacek/backup 
done 



3* v_data,h 

START PROLOGUE 

sees IDENTIFICATION: @(#)v_data.h 1.1 04/24/98 
RECORD OF CHANGES: 

«CHANGE NOTICE» Version 1. 1 (29 Apr 1998) - Kyongsuk Pace 
Initial submission 

END PROLOGUE 

The Fortran include file V_DATA will hold all the data definitions 
used throughout the verobs and its subroutines. This will 
help in modifying at one point if the requirement happens to 
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! change in the future. 

integer :: im 
integer ::jm 

integer :: ijmax ! 1 -dim array max size 

integer ::maxobs Imaxnumofobs 

integer : : maxprm ! max num of parameters to verify 

integer :: maxstat ! max num of stats to verify 

integer :: maxtvl ! max num of threshhold values 

integer :: maxlvl ! max num of levels 

integer :: size 

real :: bad_value ! value for missing or bad data 
real :: check_val ! ISIS missing value checking number 

parameter(im=360) 
parameter(jm= 181) 
parameter(ijmax=im*jm) 
parameter(maxobs=50000) 
parameter(maxprm=20) 
parameter(maxstat= 1 0) 
parameter(maxtvl=20) 
parameter(maxlvl=3 0) 
parameter(size=5000) 
parameter(bad_value=l .E+10) 
parameter(check__val=l .E+9) 

! 

! env var 
! 

character(16) :: crdate_val ! current run dtg value 
character(8) :: taui_val ! starting tau value 
character(8) :: taue__val ! ending tau value 
character(8) :: tauinc_val ! tau increment value 
character(32) :: modelname ! model name value 
character(32) :: geomname ! geometiy name value 
character(24) :: pijnnm ! projection name 
character(16) :: vdtg ! verifying dtg, also for write 
integer :: ngeom ! C pointer for the given geomname 
integer ncols I number of columns 
integer :: nrows ! number of rows 
integer :: itaui ! integer beginning tau 

integer :: itaue ! integer ending tau 

integer :: itauinc ! integer tau increment 

4. verobs.f90 

program verobs 

START PROLOGUE 

CSCCS IDENTIFICATION: @(#)verobs.f901.1 04/24/98 /h/cm/library/mverif/src/main/verobs.f90_v 
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c 

C CONFIGURATION IDENTIFICATION; NONE 
C 

C MODULE NAME: verobs 
C 

C DESCRIPTION; This program verifies the model forecast fields against 
C observations for various models. 

C 

C COPYRIGHT; (c) 1998FLENUMMETOCCEN 

C U.S. GOVERNMENT DOMAIN 

C ALL RIGHTS RESERVED 

C 

C CONTRACT NUMBER AND TITLE: N/A 
C 

C REFERENCES; previous verobs.f 
C 

C CLASSIFICATION: Unclassified 
C 

C RESTRICTIONS: NONE 
C 

C COMPUTER/OPERATING SYSTEM 
C DEPENDENCIES: CrayUNICOS 

C 

C LIBRARIES OF RESIDENCE: /a/ops^in 
C 

C USAGE: Kom shell script mverifjob 
C 

C PARAMETERS: N/A 
C 

C COMMON BLOCKS: N/A 
C 

C FILES; 

C Name Unit Type Attribute Usage Description 

C 

C SMODEL-verfil 10 FORMATTED DIE^CT IN Contains parameter, stat 
C types to compute 

C SMODEL-sfcverfil 10 FORMATTED DIFECT IN Contains parameter, stat 
C types to compute for sfc 

C MODELGEOMDTG 10 FORMATTED DIRECT OUT Contains the computed stats 
C (e.g.,NORAPSconus_nestl_appll996010100) for each run 
C 

C DATABASES: 

C Name Table Usage Description 

C 

C ISIS Grid data Various IN Model Forecasts 

C ISIS LLT data Various IN Observed environmental data 

C 

C NON-FILE INPUT/OUTPUT; N/A 
C 

C ERROR CONDITIONS: 

C CONDITION ACTION 

C 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 



no env variables exit with code 1 
err in char to integer 
conversion exit with code 2 

empty input arrays exit with code 3 
err opening output file exit with code 4 
ggrd error exit with code 5 

getgeom error exit with code 6 
boundary error exit with code 7 
no ISIS Grid data stop executing 
no ISIS LLT data stop executing 

ADDITIONAL COMMENTS: NONE 



..MAINTENANCE SECTION.. 



MODULES CALLED; 

Name Description 



BOUNDARY Computes min/max lat/lon for reading obs from LLT DB 

CH2INT Converts character to integer 

DBSTOP ISIS software that terminates database 

DTGMOD FNOC utility that increments a DTG 

EXIT System call that exits program 

LLTREAD Reads obs from LLT DB 

PXFGETENV Gets enviroimient variable 

GETGEOM Gets geometry arguments to be used by other routines 

GGRD Returns ISIS info, on given geometry 

GRD Reads gridded fcst fields from ISIS 

STRLEN FORTRAN function that returns string length 

F20B Interpolates fcst to obs pts 

UV2DF Converts wind u/v to direction and speed 

LOCAL VARIABLES AND Structures are documented in detail 
STRUCTURES; where they are defined in the code 
within include files. 

INCLUDE FILES; 

Name Description 



COMPILER DEPENDENCIES: empef90 
COMPILE OPTIONS: -f fixed -c 

MAKEFILE; Located at /a/ops/app/mverifra'c/main/Makefile 
UNICOS make 

RECORD OF CHANGES: 

«CHANGE NOTICE» Version 1.1 (29 Apr 1998) — KyongsukPace 
Initial submission 

END PROLOGUE 
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c 

implicit none 
include V_datah' 

4:* **********:):* ^************:»::4::ic*4:********3)c^**:4E******3)c***:4E*:ic3)c*:lc:lc*:lc3lc3)c3)e4e* 

c Local variables 

^ :ic :)c 3ie :|e 3|C 4: :(c * :)c :|c :)c 4e :)e * :)e :ic D|C * :^ ***** 3|e Die :)e :4c ******* ************:«:* * 4c ************ :)e :1c **** ic 

! record type 
TYPE ver_data 
character(32) :: param 
character(24) :: dsetnm 
character(32) :: obsjaram 
character(24) :: odsetnm 
character (3 2) :: unit 
character(24) :: typlvl 
real :: lvl_l 
character(24) :: obs_type 
END TYPE ver_data 

^****************** ****************************************************** 
! array of records 

^***** ************************************************************ ******* 
TYPE (ver_data) verif(maxprm*maxlvl), 

2 sfc_verif(maxprm) 

character(l) :: cnul lone blank space var used to initialize 
character(8) :: seclvl ! secret level returned from GRD 
character( 1 0) : : dtg 110 char long crdate_val 
character(16):: cdtg Idtgmod applied dtg, i.e. dtg-0,12,24, etc. 
character(16):: tdtg Itemp cdtg 

character(16):: stats(maxstat), sfc_stats(maxstat) Istat types 
character(24):: stdesc ! storage description returned from GETGEOM 
character(24):: Ivltype(maxprm), sfc_lvltype(maxprm)!level types 
character(24):; obstype(maxprm), sfc_obstype(maxprm) 
lilt seq type (raob_qc, sfc_ship, etc.) 
character(24):: dsetname(maxprm), sfc_dsetname(maxprm) 
character(24):: odsetname(maxprm), sfc_odsetname(maxprm) 
character(32):; parm(maxprm), sfc_parm(maxprm) Iparameters 
character(32):: obs_parm(maxprm), sfc_obs_parm(maxprm) 
character(32):; units(maxprm), sfc_units(maxprm) 
character(32):: geomname2, outgeomname 
character(60):: outstats ! output filename 
character(40) : : namlist_file 1 _val, namlist_file2_val 
character(80):: title 
character(4) :: nul_geom 

integer :: nstat, nparm, nobs, nobstype, nlevel 
integer :: sfc_nstat, sfc_nparm, sfc_nobs, sfc_nobstype 

integer :: arr_size, sfc_arr_size 

integer :: ktau, Itau, Istat, i, j, k, 1, m, n 
integer :: numchar, istat 
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integer :: id, iseq, status 



integer :: lenMODEL, lenmodelname 

integer : : lenGEOMNM, lengeomname 

integer :: lenGEOMNM2, lengeomname2 

integer :: lenTAUI, lentaui_val 

integer :: lenTAUE, lentaue_val 

integer :: lenTAUINC, lentauinc_val 

integer :: lenCRDATE, lencrdate_val 

integer : : lenNAMLIST_FILE 1 , lennamlist 1 _val 

integer :: IenNAMLIST_FrLE2, lennamlist2_val 

real : : fcst(ijmax), fcstu(ijmax), fcstv(ijmax) 

real fcst2(ijmax), fcstu2(ijmax), fcstv2(ijmax) 

real : : fii(ijmax), fv(ijmax), fdir(maxobs) 

real :: oblat(maxobs), oblon(maxobs), obval(maxobs) 

real newlat(maxobs), newlon(maxobs), newobs(maxobs) 

real :: fob(maxobs) 

real :: newfob(maxobs) 

real :: newfu(ijmax), newfv(ijniax) 

real : : level(maxlvl), sfc_level(maxlvl) 

real : : origx, origy, parml , paim2, parm3 , xintdis, 

2 yintdis, minlat, maxlat, minion, maxlon 

real : : xlvl, level__2, paknul 

real :: tval, sfc_tval 

real : : bias, rms, std, ancor 

real :: flau ! float tau 

double precision :: origlat, origlon 

c Function 

4 : *******************:<:*♦* ************************** *************** 

integer strlen 

^:4c**:4c:4;:)c* ********* **4 e **:<(********************************** ****:)c**:)c:^ :);*** 

c Data Initialization 

^* ***************************************************************** ****** 

data cnul P 7 

data Istat 710/ 
data paknul 71 0.e 1 07 
data level_2 70.07 

^***************** ****************************************** ************* 

c Namelists 

^***** ************ ******************************************************* 
namelist 7verlst7 parm, dsetname, obs_parm, odsetname, 

2 units, Ivltype, obstype, level, stats, tval 

namelist 7sfclst7 sfc_parm, sfc_dsetname, sfc_obs_parm, 

2 sfc_odsetname, sfc_units, sfc_lvltype. 
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2 sfc_obstype, sfc_level, sfc_stats, sfc_tval 

^4: 4: 4: 4: 4: 4: 4: 4; 4: 4: 4: 4: 4: 4c 4: 4: 4: 4: 4: 4: 4: 4: 4: 3 ): 4; 4; 4: 4: ^ 4: 4c 4: 4: :)c 4c 4c 4c 4c 4c 4c 4c 4: 3 ^ :}: 4; 4: 

c initialize some var 

q4c 4:4: 4: 4c 4c 4: 4: 4c 4: 4: **4: 4: *4: **4: ****4: 4: 4: 4: 4: 4:4: 4: 4:4c 4: 4:4: 4c *4: ****4: **********4:*** *4:* 4: *4: 4: 4c 4c 4: **4: 

crdate_val( 1:16) = cnul 
dtg(lrlO) = cnul 
cdtg(l:16) = cnul 
seclvl = ’UNCLASS' 

^ 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c ** 4c 4: 4: 4: 4= * 4: 4: 4: 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4= 4c 4c 4c 4c * 4c 4c 4c 4: 4c 4c 4c 4c 4c 4: 4: 4: 4c 4c ♦ 

c Get the environment variables that are set in the job script, 

c CRDATE, TAUS, TAUE, TAUINC, MODEL, GEOMNM, DATASETl, DATASET2, 
c NAMLIST_FILE 1 and NAMLIST__FILE2. 

^4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

namlist_filel_val = cnul 
namlist_file2_val = cnul 
geomname2 = cnul 

CALL PXFGETENV('ISIS_TABLE', lenMODEL, modelname, lenmodelname, 

2 istat) 

if (istat .ne. 0) then 
write *, "MODEL is unspecified" 

CALLEXIT(l) 
end if 

CALL PXFGETENV('GEOMNMr, lenGEOMNM, geomname, lengeomname, 

2 istat) 

if (istat /= 0) then 

write *, "GEOMNMl is unspecified" 

CALLEXIT(l) 
end if 

CALL PXFGETENV(’GEOMNM2’, lenGEOMNM2, geomname2, lengeomname2, 

2 istat) 

if (istat /= 0) then 

write *, "No second GEOMNM is specified" 
end if 

CALL PXFGETENV(’TAUr, lenTAUl, taui_val, lentaui_val, istat) 
if (istat /= 0) then 
write *, "TAUI is unspecified" 

CALL EXIT(l) 
endif 

CALL PXFGETENV('TAUE’, lenTAUE, taue_val, lentaue_val, istat) 
if (istat /= 0) then 
write *, "TAUE is unspecified" 

CALLEXIT(l) 

endif 

CALL PXFGETENV('TAUINC, lenTAUINC, tauinc^^val, lentauinc^val. 
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2 istat) 

if (istat /= 0) then 
write *, "TAUINC is unspecified” 

CALLEXIT(l) 

endif 

CALL PXFGETENV(’CRDATE’, lenCRDATE, crdate_val, lencrdate_val, 

2 istat) 

if (istat /= 0) then 
write *, "CRD ATE is unspecified" 

CALL EXIT(l) 
endif 

CALL PXFGETENVCNAMLIST__FILE1 lenNAMLIST^FDLEl , 

2 namlist__filel_val, lennamlistl_val, istat) 

if (istat /= 0) then 

write *, "NAMLIST_FILE1 is unspecified" 
endif 

CALL PXFGETENVCNAMLIST_FILE2', lenNAMLIST_FlLE2, 

2 • namlist_file2_val, lennamlist2_val, status) 

if (status /= 0) then 
write *, "no second namlist file." 
endif 

if (istat /= 0 .and. status /= 0) then 
write *, "cannot continue without the namlist files." 

CALLEXIT(l) 

endif 

9t;4;4;:i}c4;*9t:*****************4:* ************** ****4:*********3ic %********* 4:* 4:*** 

c get ISIS info on the given geometry by calling GGRD 
c geom is defined as a SMS env var 
c will add the capabihty to add the geom info in the future 
c when the geom info is not defined in ISIS 

Q***** **************************************************** *************** 

CALL GGRD(geomname, ngeom, istat) 

!if geom is not defined then exit the program 
if (istat .ne. 0) then 

write (0,*("undefined geom, ggrd returns istat =",i5)') istat 
CALLEXIT(5) 
end if 

print *, 'taui=', taui_val, ' taue=', taue__val, 

2 ' tauinc=', tauinc_val 

print *, 'model-, modelname 
print *, 'geomnm=*, geomname 

Q* *********************************************************************** 

c need to ojnvert tau info to integers 
CALL CH2INT(taui_val, itaui, istat) 
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if (istat ,ne. 0) then 

write (0,'("ch2int on TAUI returns istat =",i5)') istat 
CALLEXIT(2) 
end if 

CALL CH2INT(taue_val, itaue, istat) 
if (istat .ne. 0) then 

write (0,'("ch2int on TAUE returns istat - ’,i5)') istat 
CALL EXIT(2) 
end if 

CALL CH2INT(tauinc_val, itauinc, istat) 
if (istat .ne. 0) then 

write (0,’C'ch2int on TAUINC returns istat =*',i5)') istat 
CALLEXIT(2) 
end if 

c Initialize the arrays 

dom=l, maxprm 
parm(m) = cnul 
dsetname(m) = cnul 
odsetname(m) = cnul 
units(m) = cnul 
Ivltype(m) = cnul 
sfcjparm(m) = cnul 
sfc_dsetname(m) = cnul 
sfc_odsetname(m) = cnul 
sfc_units(m) = cnul 
sfc_lvltype(m) = cnul 
obstype(m) = cnul 
sfc_obstype(m) = cnul 
end do 

dom=l, maxlvl 
level(m) = paknul 
sfc_level(m) = paknul 
end do 

don=l, raaxstat 
stats(n) = cnul 
sfc_stats(n) = cnul 
end do 

lOpen and read the NAMLIST_FILE1 into the arrays for parm, 
lunit, lvl_type, level_l, and obs_type 

! determine the number of parameters and number of levels 
♦.number of stats and number of obstype 
♦determine the array size 
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********* *********************V*****3t:****:)::|e* ******** 4:^4: *3):* *****3):* *3): 4:4: * 

if (namlist_filel_val /= cnul) then 
open (unit=lstat,file=namlist_file l_val,form='formatted', 

2 status='old' ,iostat=istat) 
if (istat .eq. 0) read(lstat, nml = verlst) 

^************************************************************************ 
c Finish setup, the arrays, stype, tval, typlvl, units are changed 
c to 1 -dimension and will not change throughout the program so that 
c they can be used for multiple parameters which is a change from 
c the original program. 

^********* ******************************************** ******************* 
c Determine the number of parameters 

^* **************************************************** ******************* 
nparm = 0 

do while (parm(nparm+l) .ne. cnul .and. nparm .It. maxprm) 
nparm = nparm + 1 
end do 

if (nparm .eq. 0) then 

write (0,'("No verification parameters specified")') 
c CALLEXIT(3) 
end if 

^* ******************************************************** *************** 
c Determine the number of levels 

^*********************** ********************************** *************** 
nlevel = 0 

do while (level(nlevel+l) .ne. paknul .and. nlevel .It. maxlvl) 
nlevel = nlevel + 1 
end do 

if (nlevel .eq. 0) then 

write (0,'("No verification levels specified.")') 
c CALLEXIT(3) 
end if 

^********* ****************************************************** ********* 
c Determine the number of stats 

^********* ************************************** ************ ******* ****** 
nstat = 0 

do while (stats(nstat+l) .ne. cnul .and. 

2 nstat .It. maxstat) 

nstat = nstat + 1 
end do 

if (nstat .eq. 0) then 

write (0,'("No statistics are requested.")') 
c CALLEXIT(3) 
end if 

end if !namlist_filel exists 
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c if NAMLIST_FILE2 exists, read it and put into the arrays for 
c the surface 

^***** 3 »:* *:►**** ********* ************************************************ * 5 *: 

if (namlist_file2_val /= cnul) then 
open (unit=lstat, file=nainlist_file2_val, 

2 form-formatted', status- old’, iostat=istat) 

if (istat .eq. 0) then 
read(lstat, nml = sfclst) 
sfc_nparm = 0 

do while (sfc_parm(sfc_nparm+l) .ne. cnul .and. 

2 sfc_nparm .It. maxprm) 

sfc_nparm = sfc_nparm + 1 
end do 
sfc_nstat = 0 

do while (sfc_stats(sfc__nstat+l) .ne. cnul .and. 

2 sfc_nstat .It. maxstat) 

sfc_nstat = sfc_nstat + 1 
end do 

if (sfc;_nstat .eq. 0) then 
write (0,’("No statistics are requested in sfc.")') 
c CALLEXIT(3) 

end if 

end if !if namlist_file2 was read successfully 
end if !if namlist__file2 exists 

^************************************************************************ 
c array size based on the namlists 

Q* ************************************************************* ********** 

arr_size = nparm*nlevel 
sfc_arr_size = sfc_npaim 

^********* ******************************************************** ******* 
!Fill the arrays of records from the parm, unit, lvl_type, 

!level_l and obs_type arrays 

^* ********************************** ************************************* 
1=1 

if (arr_size .gt. 0) then 
do i = 1 , nparm 
doj = l,nlevel 
verifO)%param = parm(i) 
verif(l)%dsetnm = dsetname(i) 
verifO)%obs_param = obs_parm(i) 
verif(l)%odsetnm = odsetname(i) 
verifG)%unit = units(i) 
verifG)%typlvl = Ivltype(i) 
verif(l)%lvl__l = level(j) 
verif(l)%obs_type = obstype(i) 

1 = 1+1 
end do 
end do 
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end if 



!add the sfc stuff to the verif array 

^:|E3|e3ic:|e:)c3|e**:|e*:tE************:|e*********** ********** ***********3|c3|c*3|e4:3|c*4:****3|e sickle 3|<: 

if (sfc_arr_size .gt. 0) then 
arr_size = arr_size + sfc_arr_size 
do i = 1 , sfc_nparm 
verif(l)%param = sfc j>arm(i) 
verif(l)%dsetnm - sfc_dsetname(i) 
verif(l)%obs_param = sfc_obs_pann(i) 
verif(l)%odsetnm = sfc_odsetname(i) 
verif(l)%unit = sfc_units(i) 
verif(l)%typlvl = sfc_lvltype(i) 
verif(l)%lvl_l = sfc_level(i) 
verif(l)%obs_type = sfc__obstype(i) 

1 = 1 + 1 
end do 
end if 

^*** *************************************************************** ****** 
!get geom info 

^*** *********************************************************** ********** 
CALL GET GEOM(ngeom, pijnnm, stdesc, ncols, nrows, 

2 origlat, origlon, origx, origy, xintdis, 

3 yintdis, parml, parm2, parmS, istat) 

if (istat .ne. 0) then 

write (0,'("getgeom returns istat =",i5)’) istat 
CALLEXIT(6) 
end if 

^**** ************************************************************** ****** 
c Determine min/max lat/lon to use in obs read 

Q**** ************************************************************** ****** 

CALL BOUND ARY(ngeom, nrows, ncols, minlat, 

2 maxlat, minion, maxlon, istat) 

if (istat .ne. 0) then 

write (0, '("cannot find min/max lat/lon")') 

CALLDBSTOP 
CALLEXIT(7) 
end if 

^*** ************************************************************* ******** 
c The output file, use modelname, geomname and crdate_val 

Q*** ******************************************************* ************** 

1 = strlen(modelname) 
outstats( 1 :1) = modelname( 1 :1) 
if (geomname2 = cnul) then 
numchar = strlen(geomname) 
k = 1 + numchar 
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outstats(H-l :k) = geomname(l :numchar) 
else 

numchar = strlen(geomname2) 
k = 1 + numchar 

outstats(l+l :k) = geomname2(l :numchar) 
end if 



outstats(k+l :k+10) = crdate_val 
write (O,'(”output file name = ", a50)’) 

2 outstats(l :strlen(outstats)) 

open(unit=lstat,file=outstats(l:strlen(outstats)), 
2 form=Tormatted‘,status-new',iostat=istat) 

if (istat .ne. 0) then 
write (0,'("istat = ", i5)') istat 
write (0,’("Cannot open output file for stats")') 
CALL EXIT(4) 
end if 



******** ****♦***♦***♦♦**♦*♦*******♦♦*♦*♦****♦♦*♦*** 
c Write the header in the output file 

^*** ************************** ****** ****** ********** ****** ******* ******** 
write (Istat, '(" vdtg ", " numobs ", 

2 " param ", " unit ", 

3 " geometiy ", 

4 " level type "," level 1", 

5 "tau","stattype "," statval ", 

6 " v_src ", " obs_type ")') 

^***** ********************************************************** ********* 
c for each parameter and level, read obs fi*om LLT DB 
^**** ********************************************************** ********** 
dtg(l ; 1 0) = crdate_val(l : 1 0) 



do i=l, arr_size 

write (0,’("parm = ",a32)') verif(i)%param 
write (0,'("unit = ",al0)') verif(i)%unit 
write (0,'("type Ivl = ",al0)') verif(i)%typlvl 
write (0,'("lvl 1 = ”,f8.1)') verif(i)%lvl_l 
write (0,'("obs type = ",al5)') verif(i)%obs_type 



2 

3 

4 

5 

6 

7 

8 



CALL LLTREAD(verif(i)%obs_type, verif(i)%obs_param, 
verif(i)%l vl_ 1 , verif(i)%odsetnm, 
dtg, minlat, 

maxlat, minion, 
maxlon, verif(i)%typlvl, 
oblat, oblon, 

nobs, obvai, 

istat) 



write (0,’("nobs = ",i6)') nobs 
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if (istat /= 0 .or. nobs = 0) then 
go to 100 
end if 

^ 3|e * :|c :|E :|e * 3|c He * * 3|c 3|: 3|: :|e 3): 3|c * :|E 3|c :|E 3|c 3|e 3|c 3tc 3|c 9^ ^ :|c 3»; :|c 4; :|c 3|c :4: :)c 4: 3|e 3|c 3|c 3)c 3|c 3|c 3^ 

c for each tau read the fcst, inteipolate the fcst to obs 
c compute the stats and output the stats. 

^3|C 3(: 3|C 3|C 3(C 3f; 3|C 3|C 3): 3(C 3t: 3|e 3^ 3): 3(C 3|C 3|e 3(C 3|C 3|C 3|C 3(C 3|C 3t: 3(C 3t: 3(C 3|C 3(C 3): 3|: 3): 3(: 3t: 3|C 3|C 3|: 3): 3|C 3|C 3|C 3(C 3|C 3|: 3^ 3|: 3t; 3»; 3^ 3)e 3|C 3^ 3|e 3|e 3|C 3f: 3»: 3|C 3|: 3|: 3|C 3|C 3|C 3|c 3t: * 

do Itau = itaui, itaue, itauinc 

! determine the correct dtg table to read 
CALL DTGMOD(dtg, -Itau, cdtg, istat) 
if (istat /= 0) then 
go to 200 
end if 

^3|C 3): 3|: 3|: 3»: 3|C 3)C 3|: 3|C 3|C 3tc 3|C 3|: 3|: 3t: 3|: 3|C 3|: 3(c 3|: 3|: 3|C 3|: 3): 3»: 3|C 3|C 3|: 3|: 3|: 3»: 3): 3|C 3»: 3|c ***************** 4: ***************** * 

c ISIS has tables for 00 and 1 2 only, therefore if we need 
c to read other tau model forecasts, e.g., 3, 6, 9, 1 5, 1 8, 2 1 
c etc. , we need to read 1 2 hour old table 
c if cdtg ends with anything other than 00 or 1 2 then 
c use - 1 2 ISIS table dtg 

^3t: * 3|C 3|C 3|C 3|C 3^ * 3|C 3): * 3^ 3t: 3|C 3t: 3): 3|C 3|C 3|C 3^ 3t: 3»; 3|C 3|C 3^ 3|C 3^ 3|C * * * * * 3^ 3^ 3^ 3)C 3|: 3|e 3|; 3|: 3|: 3|: 3|C * 3|C * 3|: 3(: * 3(: 3(C 3»: 3|C 3|C 3t: 3(C 3»: 3^ 3|C 3(: 3(: 3t: 3»: 3»: 3|: 3(: 3f: 3^ 3|: 3|C * 

if (cdtg(9.i0) /= '00' .and. cdtg(9:10) /= '12') then 
CALL DTGMOD(tdtg, -12, cdtg, istat) 
if (istat /= 0) then 

write *, 'DTGMOD error for tau ', Itau 
go to 200 
end if 
else 

tdtg = cdtg 
end if 

write (0,’("ISIS table dtg =",al0)') cdtg 
write (0,'("tau =",i5)') Itau 
ftau = Itau 

^3|C 3»; 3|C 3|e 3)C 3); 3): 3|C 3|C 3|C 3|C 3»: 3|e 3): 3|C 3|C 3): 3)C 3|: 3»: 3|C 3tc 3tc 3|: 3|: 3t: 3»: 3tc 3): 3)C 3»: 3|C 3|C 3|C 3|C 3|C 3|C 3t: **** 3|C 3»; 3|C H: ************************* * 

c if everything is OK, then read the forecast 
c (gridded data) by calling ISIS GRD 
c ISIS grid does not have wnd_spd, therefore have to 
c read wnd_ucmp, wnd_vcmp then compute wnd__spd 

^3|; 3|C 3|C 3|C 3|; 3|; 3|C 3|C * * 3|C * * 3|C * * 3|C 3|C 3|; 3|C 3|; 3^ * 3|; 3); 3|; 3|e * 3|C 3^ * 3|C 3»: * 3|C 3f: * * * * 3|C 3|E 3tc * 3)C * * 3»: 3)C 3^ 3^ 3|C 3(C 3): * * * * * * * * * * 3|; 3|; 3|C * * 3|C * 3|C 

if (verif(i)%param = 'wnd_spd') then 



Q ■ — T-7 : ==i; 

c read forecast for wnd_ucmp; ISIS grid data 

c ===--’ ■ ' ' — ■ = 

CALL GRD(modelname, geomname, 

2 verif(i)%dsetnm, 'wnd_ucmp', 

3 verif(i)%typlvl, verif(i)%lvl_ 1 , 

4 level_2, cdtg, 

5 ftau, verif(i)%unit. 
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6 fcstu, title, 

7 seclvl, id, 

8 . iseq, status) 

if (status /= 0) then 
go to 200 
end if 



c= 

c 

c= 



2 

3 

4 

5 

6 

7 

8 



read forecast for wnd_vcmp: ISIS grid data 



CALL GRD(modelname, geomname, 
verif(i)%dsetnm, 'wnd_vcmp', 
verif(i)%typlvl, verif(i)%lvl_l. 



level_2, 
ftau, 
fcstv, 
seclvl, 
iseq, 

if (status /= 0) then 
go to 200 
end if 



cdtg, 

verif(i)%unit, 

title, 

id, 

status) 



C 

c first interpolate into another geometry if the second 

c geometry is requested 

C ===T -- — -■■T-" — ■■::-:-rT= = = =u 

if (geonmame2 /= cnul .and. geomname2 /= geomname) then 
CALL INTGEOM(geomname, geomname2, ijmax, fcstu, 

2 fcstu2, istat) 

if (istat = 0) fcstu = fcstu2 

CALL ENTGEOM(geomname, geomname2, ijmax, fcstv, 

2 fcstv2, istat) 

if (istat = 0) fcstv = fcstv2 
CALL GGRD(geomname2, ngeom, istat) 
end if 



C ~~=~~'~:— ' ■ = ■-=== = == - ■ - --- = 

c interpolate forecast into obs points, then 

c convert wnd_ucmp and wnd_vcmp to wnd_spd by calling 

c uv2df 

C - ■ - ■-■■ ■ -- — 

CALL F20B(ngeom, fcstu, oblat, oblon, nobs, fu, istat) 

CALL F20B(ngeom, fcstv, oblat, oblon, nobs, fv, istat) 

^ :fc :(c :|c :(c :ic 3|c :4c :tc :<c :|e :|c 3|c :|c :|c :|c :fc :|c :|c :|c :|c :|c :|c :|c :(c :(c :|e :|c :|c 3|c :|c :|e :fc :|c :fc :|c :|e :(c :|c :|e 3^ * 

c if fcst had any missing value, fZob fills fob with 

c -999 which needs to be disregarded in the conversion 

:|c:|c **♦** + ********♦*♦**** Jjc ****♦*♦*♦**♦ 3jc 3jc 3|e 3je 3|e 3(C 3|C 3|C 3jc 3|C 3jc :(£ 3je 3|e 3(C 3|e + 3je 3|e 3|C 5|C 3|C *51: 3jc * 3jc 3j£ 3jc 

k = 0 

do j=l, nobs 

if (fu(j) /= -999 .and. fv(j) /= 999) then 
k = k-i-l 
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newfli(k) = fu(j) 
newfV(k) = fv(j) 
end if 
end do 

CALL UV2DF(newfli, newfv, nobs, fdir, fob) 

^:|c:4c:lc**********3|e*3ic *********** ***3|e* ********* 3|c9|c3|c3ic3(c3|c:|c3|c:tc*********3|c****4c*****9|e* 

c If the parameter is other than winds, 
c call GRD only once 

^ 4c :|c 3|c :|e :|c :|c :|c ******************************** :|c * :|c ******* ic * :|c 3|c :|e :|e:)e :|c 3|c 4c 3|c :|c :|e:|c :|c :|c :|c 4c :|c :|e :|c 

else 



c= 

c 

c= 



2 

3 

4 

5 

6 

7 

8 



read forecast: ISIS grid data 



CALL GRD(modehiame, geomname, 
verif(i)%dsetnm, verif(i)%param, 
verif(i)%typlvl, verif(i)%lvl_l , 
level_2, cdtg, 
ftau, verif(i)%unit, 

fcst, title, 
seclvl, id, 

iseq, status) 



if (status /= 0) then 
go to 200 
end if 



C ; -:T7=-^rr 1: : ~ 

c first interpolate into another geometry if the second 

c geometry is requested 

Q :ir=-:=7r-=-==rT? ■■■ ■ . 

if (geomname2 /= cnul .and. geomname2 /= geomname) then 
CALL INTGEOM(geomname, geomname2, ijmax, fcst, 

2 fcst2, istat) 

if (istat = 0) fcst = fcst2 
CALL GGRD(geomname2, ngeom, istat) 
end if 

^***** ****************************************************** ******* ****** 
! interpolate fcst into obs pts 

^**** ************************************************************** ****** 
CALL F20B(ngeom, fcst, oblat, oblon, nobs, fob, istat) 
if (istat /= 0) then 
go to 200 
end if 

end if Iparam test 

^***** ************************************************ ******************* 
c if fcst had any missing value, fZob fills fob with 

c -999 which needs to be disregarded in the stat computations 
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^ :|c 3|c :|c 3|c :)c :|c :|c 3|C :)c 3(c 3|c 4: :|c :|e 4: ♦♦ ♦♦ 3|c :)( :)c;4: :)c :)c :)c :)c :)(:)( ***♦*♦♦ :)c :|e :)c :)c :(c 3|c 3»: 

k = 0 

do j=l, nobs 
if (fob(j) /= -999) then 
k = k+ 1 

newfob(k) = fob(j) 
newlat(k) = oblat(j) 
newlonQc) = oblon(j) 
newobs(k) = obval(j) 
end if 
end do 

write (0,’C’writmg the 1st 15 lat,lon,obs/cst.")') 
write (0,’(i5,4fl0.2)’) (njiewlat(n)^ewlon(n)^ewobs(n) 

2 ,newfob(n), n= 1 , 1 5) 

! compute the stats 

!for obs, disregard the projection 

^ + 4: :|e :|c + 3|e :|c :|e 3|c :|c ♦ 3|c :|c :Jc :Je 3|c ♦ 3|c ♦♦♦ a|c * a|c ♦ 3|c 3|c 4: 3|c 3|c 3|c :)c 3|c 3|c ♦*♦♦*♦♦♦♦♦** + ♦*♦*♦***♦ ♦ 

nul_geom = 'NONE' 

^ :fc :|c 4: :|e :|e :|c 3|c 3|c 3|e 3|c :|c 3|c 3|c ^ 3|e 3|c :|c :|e 3|e 3(c :|c 4: 3|c 3|c 3|c :|c 4; :)c :i)e 4: 3|c :|c 4: :|c :i»: 9|c :)e 3|c :|c :)e 3|c :)e 3|c 4; 3)e 3|c :)e 3|c 3|e :)e 3|e 3|c ^ 

c if there was no NAMLIST_FILE 1 , then use the stats 

c fromNAMLIST_FILE2. For ocean models 

^4:4c *4:4:4: ♦ Jit**** 5|6*3|c* ♦♦***♦ ♦ 

if (nstat = 0) then 
nstat = sfc_nstat 
do j = 1, nstat 
stats(j) = sfc_stats(j) 
end do 
end if 

if (geomname2 /= cnul) then 
outgeomname = geomname2 
else 

outgeomname = geomname 
end if 

doj = 1, nstat 

if (stats(j) = ’bias' .and. k /= 0) then 
CALL COMPUTE_BIAS(newfob,newobsJk,nul_geom,bias) 

! write the stats to the output file 
write Ostat, 1000) dtg, k, verif(i)%param, 
verif(i)%unit, 

outgeomname, verif(i)%typlvl, 
verif(i)%lvl_l, Itau, 
stats(j), bias, 
verif(i)%obs_type 

write (0,'(a30, 2al5,fl0.2)') verif(i)%obs_type, 
verif(i)%param, stats(j), bias 



2 

3 

4 

5 

6 

2 
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else if (stats(j) = 'std' .and. k /= 0) then 
CALL COMPUTE_STD(newfob, newobs, k, nul_geom, std) 
write (Istat, 1000) dtg, k, verif(i)%param, 

2 verif(i)%unit, 

3 outgeoniname, verif(i)%typlvl, 

4 verif(i)%l vl_ 1 , Itau, 

5 stats(j), std, 

6 verif(i)%obs_type 

write (0,'(a30, 2al5,fl0.2)’) verif(i)%obs_type, 

2 verif(i)%param, stats(j), std 

else if (stats(j) = ’rms’ .and. k /= 0) then 
CALL COMPUTE_RMS(newfob, newobs, k, nul_geom, rms) 
write (Istat, 1000) dtg, k, verif(i)%param, 

2 verif(i)%unit, 

3 outgeomname, verif(i)%typlvl, 

4 verif(i)%lvl_l , Itau, 

5 statsQ), rms, 

6 verif(i)%obs_type 

write (0,’(a30, 2al5,fl0.2)’) verif(i)%obs_type, 

2 verif(i)%param, stats(j), rms 



c else if (stats(j) = ’ancor' .and. k /= 0) then 

c CALL COMPUTE_RMS(newfob,newobs,k,nul_geom,ancor) 

c write (Istat, 1000) dtg, k, verif(i)%param, 

c 2 verif(i)%unit, 

c 3 outgeomname, verif(i)%typlvl, 

c 4 verif(i)%lvl_ 1 , Itau, 

c 5 stats(j), ancor, 

c 6 verif(i)%obs_type 

c write (0,’(a30, 2al 5/10.2)’) verif(i)%obs_type, 

c 2 verif(i)%param, stats(j), ancor 



! other stat types ... 
end if 

end do Ido j=l ,nstat 
200 continue 
end do lltau loop 
1 00 continue 
end do I i= 1 ,aiT_size loop 



1000 format (al0,i6,lx,2al5,a30,al5,f8.2,i5,lx,al2,f8.2,’ obs',lx,al5) 
close(lstat) 

CALL DBSTOP 
stop ’Normal End* 

CALLEXIT(0) 

end 



5. boundary.f90 

subroutine boundary(ngeom, nrows, ncols, minlat. 
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2 maxlat, minion, maxlon, istat) 

C 

C START PROLOGUE 

C 



C sees IDENTIFieATION: @(#)boundary.f90 1.1 04/24/98 /h/cm/library/mverif/src/sub/boundary.f90_v 

e 

e eONHGURATION IDENTIFieATION: NONE 

e 

e MODULE NAME: boundary 

e 

e DESeRIPTION: This subroutine computes the minimum and maximum 
e latitude and longitude for reading observations 

e from ISIS LLT database. 

e 

e eOPYRIGHT: (c) 1996FLENUMMETOeeEN 

e U.S. GOVERNMENT DOMAIN 

e ALL RIGHTS RESERVED 

e 

e eONTRAGT NUMBER AND TITLE: N/A 

e 

e REFERENGES: NONE 

e 

G GLASSIFIGATION: Unclassified 
G 

G RESTRIGTIONS: NONE 

e 

G GOMPUTER/OPERATING SYSTEM 
G DEPENDENGIES: Gray UNIGOS 

G 

G LIBRARIES OF RESIDENGE: /a/ops/bin 
G 

G USAGE: 

G call boundary(igeom, nrows, ncols, minlat, 

G maxlat, minion, maxlon, istat) 

G 

G PARAMETERS: 



c 

c 


Name 


Type Usage Description 


G 


NGEOM 


INTEGER 


INPUT 


Geometry info. 


G 


NROWS 


INTEGER 


INPUT 


No. of rows 


G 


NGOLS 


INTEGER 


INPUT 


No. of columns 


G 


MINLAT 


REAL 


OUTPUT 


Minimum latitude 


G 


MAXLAT 


REAL 


OUTPUT 


Maximum latitude 


G 


MINLON 


REAL 


OUTPUT 


Minimum longitude 


G 


MAXLON 


REAL 


OUTPUT 


Maximum longitude 


G 


ISTAT 


INTEGER 


OUTPUT 


Return status 


G GOMMON BLOGKS: N/A 






G FILES: N/A 








c 


DATA BASES: N/A 
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c 

C NON-FILE INPUT/OUTPUT: N/A 
C 

C ERROR CONDITIONS: 

C CONDITION ACTION 

C 

C Error return from VXYLL Print err message 
C 



C ADDITIONAL COMMENTS: NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED: 

C Name Description 

C 

C MAXVAL Returns maximum value from an array 
C MINVAL Returns rrimimum value from an array 
C VXYLL Computes arrays of lat/lon from arrays of yJy 

C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES: where they are defined in the code 

C within include files. 

C 

C METHOD: 1 . Fill in the working arrays, x & y, with values for left, 

C right, bottom & top boundaries. 

C 2. Call vxyll to get lat/lon from x & y. 

C 3 . Find the minimum & maximum latitude & longitude. 



C 

C INCLUDE FILES: NONE 
C 

C COMPILBR DEPENDENCIES: 190 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mverif/src/sub/makeverobslib 
C UNICOS make 
C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1.1 (29 Apr 1998) — Kyongsuk Pace 
C Initial submission 
C Originial Programmer: M.A. Rennick 
C 

C END PROLOGUE 

C 

implicit none 

1 ********************** 4c****** *********:<(*** ****** 

! parameters 

I **** ****^* ****** ****4c*****^*******^***4;^:tc*:jc*****^** *:)(**** ******** 

integer, intent(in) :: ngeom ! geometry 
integer, intent(in) :: nrows ! no of rows 
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integer, intent(in) ncols ! no of columns 
real, intent(out) :: minlat ! min lat 
real, intent(out) :: maxlat ! max lat 
real, intent(out) :: minion ! min Ion 
real, intent(out) :: maxlon ! max Ion 
integer, intent(out) :: istat ! return status 

:|c :1c :)e 3)c :|c ^ 4: 3): * ^ 3): 4: * * 4: * * * ^ :|c 

local variables 

3|C 3|C 3|C 3|C 3|C 3|C 3|C 3|C 3|C 3|C 3|C 3|C 3(C 3): 3): 3)C 3}C 3|C 3(C 3»; 3|C :1c 3)C 3|C 3»; 3|C 3lc 3|C 3|C 3|C 3)C 3): 3): 3|C 3|C 3»; 3|C 3|C 3|C 3lc 3}C 3): 3): 3*: 3|C 3|C 3)C 3|C 3*: :4c 3)C 3|C 3)C 3»; 3lc 3)C 3)C 3)C 3|C 3)C 3|C 3|C 3|C 

integer :: len ! array size 

integer :: i ! dummy loop var 

integer : : n ! temp var to hold array value 

real :: x(2*nrows + 2*ncols) ! working array 1 st dimension 

real : : y(2*nrows + 2*ncols) ! working array 2nd dimension 

real :: lat(2*nrows + 2*ncols) 1 latitude array 

real :: lon(2*nrows + 2*ncols) ! longitude array 

3lc 3lc 3lc 3*: 3»: 3*: 3j: 3lc 3lc :1c 3lc :1c 3j: 3lc 3lc 3lc 3lc 3lc 3*: 3»: 34: * 3»: * * 34: * 34: 31: * 3»c 3»: 3lc :4: 3»: ♦ 3*: 3lc 3»: 3»: 3*: ♦ 3»: 3ic 3lc 3lc 3»: 34: 3*: 3»: 3*: 3lc 3*: :1c 3*: :1c :1c 3lc 3*: 3lc :1c :1c 3lc 3lc 

for left boundary 

:1c :4: :1c :► :1c :4: :► :4: :1c :► :1c 34: :4: ♦ 34: 34: * 34: :4: :4: ♦ 34: 34: :1c :►:► :1c :4: :4: :4: :1c + 34: 34: :► :1c 34: :4: 34: :1c 34: :4: 34: :1c ♦ :1c :1c :4: :1c :1c :1c :4: 

do i=l, nrows 

x(i) =1 ! all X value at the left boundary is 1 
y(i) = i ! y value at the left boundary 
end do 

:1c :1c :1c :1c :1c :1c 34: :1c :1c :1c :1c :1c 34: :4 34: :1c :1c :1c 34: 34: * :1c * ^4: 4c :4: 4c :4: * * :4: ♦ * :4: :1c :1c :1c 34: 34: :1c 3lc 34: 34: :4: * * 4: 4: 4: * 4c * * * :4: * * 4c * * * 

for right boundary 

3lc 4: 3lc 3lc 3lc 3lc 4c 4c 4: 4c 4c 3lc 4c 3lc 4c 3lc 4c 4: 4: 4: 4: 4c 4c 4: 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4: 4: 4c 4: 4: 4c 4: 4: 4c 4c 4c 4c 4c 4c 4: 4c 4c 4: 4: 4: 4c 4c 

n = 0 ! initialize the temp var 
do i=nrows+l , 2*nrows 
x(i) = ncols ! all x value at the right boundary 
n =n + 1 

y(i) = n ! y value at the right boundary 
end do 

4: 4: 4c 4c 4c 4c 4c 4: 4: 4: 4: 4: 4: 4: 4: 4c 4: 4: 4: 4: 4: 4: 4c 4c 4: 4c 4c 4: 4c 4: 4: 4: 4: 4c 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4c 4: 4: 4: 4: 4: 4: 4: 4c 4: 4: 4c 4c 4c 4c 4: 4c 

for bottom boundary 

4c 4c 4c 4c 4c 4c 3*: 4c 3*: 3»: 3»: 4: 3»: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4c 4c 4: 4c 4: 4: 4c 4: 4: 4c 4: 4c 4: 4c 4c 4: 4: 4c 4c 4c 4: 4c 4c 4c 4c 4: 4: 4: 4c 4c 4c 4c 4c 4c 4c 4c 

n = 0 ! initialize the temp var 
do i=2*nrows+l, 2*nrows-hicols 

n = n + 1 

x(i) = n ! all X value at the bottom boundary 
y(i) =1 ! y value at the bottom boundary 

end do 

3»: 3>: 3»: 4: 3t: 3»: 3»: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4: 4c 4: 4: 4: 4c 4: 4: 4: 4: 4c 4c 4: 4c 4: 4c 4: 4: 4c 4: 4c 4: 4: 4: 4: 4: 4: 4: 4c 4c 4c 4c 4c 4: 4: 4c 4c 4: 4c 4c 4c 4c 4: 

for top boundary 

3»: 3»: 3»: 4: 4c 4: 4c 4c 4: 4c 4c 4c 4c 4c 4: 4: 4: 4: 4: 4c 4: 4c 4: 4: 4: 4: 4c 4: 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4: 4c 4c 4: 4: 4: 4c 4c 4c 4c 4c 4c 4c 4c 

n = 0 ! initialize the temp var 
do i=2*nrows+ncols+l, 2*nrows+2*ncols 
n = n+ 1 
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x(i) = n ! all X value at the bottom boundary 
y(i) = mows ! y value at the bottom boundaiy 
end do 

4: * :^c :Je 4: * sf: :itc 4: * 4: * * * 4: ♦ :^e 4: ♦ + * * * * * * * * * * * * * * * * * ♦ 4: 4: 4: + * * * ♦ * * * 4: * * * * * * * ♦ * ♦ * * * * 

call vxyll to get the lat/lon from yJy 

len = 2*mows + 2*ncols 

call vxyll(ngeom, len, x, y, 'd’, lat. Ion, istat) 

if (istat .eq. 0) then 
get the min/max lat/lon 

♦ * 4 : * 4 : 4 : 4 : 4 : 4 ; * 4 : 4 : ***♦*♦♦♦♦♦******♦**♦♦♦♦♦♦*♦♦♦♦** ♦♦*♦♦♦♦♦♦*♦♦♦♦♦♦*♦♦*♦ 

minlat = minval(lat) 
maxlal = maxval(lat) 
minion = minval(lon) 
maxlon = maxval(lon) 

else 

write (*,'("vxyll returns istat =",i5)') istat 
end if 

return 

end subroutine boundary 



6. f2ob.f90 

subroutine f2ob(igeom,field,obslat,obslon,nobs,fob, istat) 
C 

C START PROLOGUE 

C 



C sees IDENTff ICATION: @(#)f2ob.f90 1.1 04/24/98 /h/cnvlibrary/mverifysrc/sub/f2ob.f90_v 

e 

e eONFIGURATION IDENTIFieATION: NONE 
G 

e MODULE NAME: Gob 
G 

G DESGRIPTION; This subroutine interpolates forecast field values to 
G observation locations. 

G 

G GOPYRIGHT: (c) 1 998 FLENUMMETOGGEN 

G U.S. GOVERNMENT DOMAIN 

G ALL RIGHTS RESERVED 

G 

G GONTRAGT NUMBER AND TITLE: N/A 
G 

G REFERENGES: NONE 
G 

G GLASSDFIGATION: Unclassified 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



RESTRICTIONS: NONE 

COMPUTER/OPERATING SYSTEM 

DEPENDENCIES: CrayUNICOS 

LIBRARIES OF RESIDENCE: /a/ops^in 

USAGE: 

call f 2 ob(igeom,field,obslat.obslon 4 iobs/ob,istat) 



PARAMETERS: 
Name Type 



Usage 



Description 



IGEOM 

FIELD 

OBSLAT 

OBSLON 

NOBS 

FOB 

ISTAT • 



INTEGER INPUT Geometry info. 

REAL INPUT Forecast array to interpolate 
REAL INPUT Obs. latitude 

REAL INPUT Obs. longitude 

INTEGER INPUT No. of observations 

REAL OUTPUT Fcst interpolated to obs array 
INTEGER OUTPUT Return status 



COMMON BLOCKS: N/A 

FILES: N/A 

DATA BASES: N/A 

NON-FILE INPUT/OUTPUT: N/A 

ERROR CONDITIONS: 

CONDITION ACTION 



Unsuccessful getgeom Print err message and exit 
ADDITIONAL COMMENTS: NONE 
MAINTENANCE SECTION. 



MODULES CALLED: 

Name Description 



EXIT System call that exits program 

FINTRP FORTRAN sub. that interpolates within a field to 
obtain values at an array of points 

GET GEOM Gets geometry arguments to be used by other routines 

IMAXCV FORTRAN fimction to determine max. 1 st dimension 
of array 

JMAXCV FORTRAN function to determine max. 2nd dimension 
of array 

VLLX Y Computes arrays of x/y fi'om arrays of lat/lon 
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C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES: where they are defmed in the code 

C within include files. 

C 

C METHOD: 1. Call getgeom. 

C 2. If successful getgeom, determine max. 1st & 2nd dimensions 
C of the array using imaxcv/jmaxcv. 

C 3. Call vllxy to convert from lat/lon to i j. 

C 4. Call fmtrp to interpolate. 

C 

C INCLUDE FILES: 

C Name Description 

C — - - 

C V_DATA.H contains the common variables for verobs 
C 

C COMPILER DEPENDENCIES: f90 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mverif/src/sub/makeverobslib 
C UNICOS make 

C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1.1 (29 Apr 1998) — Kyongsuk Pace 
C Initial submission 
C Original Programmer: MA Rennick 
C 

C END PROLOGUE 

C 

implicit none 
include V_data.h’ 

I 3|c 3|c ^ :|c 3|c 3|c :|c 3|c 3|c 3|c :|c ************* ****************** 3|c :)c 3tc * :|C 3|c ^ 3|c a|c :)c :|c ^ 3|c 3)c 3|c 3|c ^ a|c 3|c 3)c :)c ^ 3|c :1c 3|c 

! Formal parameters 

|3|e^:)c4::)c:)c:lc:)c:)c:)::)cic:)c:)c:)c4!:)c:)c:)c:)c4::|c:)c:)c:)c4c:)c4c:lc:lc:)c:l:4;:)c4::)c:)c:)c:l::)e:)c:)c:)c:|c:|c:lc:|c:)c:lc:)c4::)c:l::)c:)c:)c:l::l::fc:|c:)c:)c:l::lc:|c4:4::)c:)c*;:lc:)c:|c 

integer, mtent(in) :: igeom ! geom ptr from ggrd 
integer, intent(in) :: nobs ! number of obs to interpolate 
real, intent(in) :: field(ijmax) ! array to interpolate 
real, intent(in) :: obslat(maxobs) ! lat of report 
real, intent(in) :: obslon(maxobs) ! Ion of report 
real, intent(out) :: fob(maxobs) ! interpolated array 
integer, intent(out) :: istat ! status 

|3)c:|c:)c:)c:tc:|c:)c^:)c:)c:|c:l::|c:(c:|c:lc:|c:l::)c:|c:)c:l::l::)c:l::(c:|c:|c:|c:tc:|c:tc:|c:tc:(c:|c:|c:|c:|c:tc:|c:|c:l;:lc:|c:(c:|c:)c:le:lc*:)c:)c^:lc:lc:l::lc:|c:|e:|c:l::tc:)c:|c:(c:tc:|c:l:***:)c 

! Local variables 

I :)c :)c :1c :|c :|c * * * * :|c :|c :|c :|c :|c :|c :|c :1c :|c :|c :«: * :|c :|c :|c :)c :|c :|c :)c :1c :1c * :|c :|c :)c :)c :(c * :|c :|c :|c :|c :1c :(c :|c :|c :1c :|c :|c * :1c * :|c :|c :)c :)c :|e :|c :)c :)c :|c :1c :|c :)c :|c :|c :)c * * * 

integer ::imx,jmx 
character(8) :: dsc 
real :: dcol 

real :: drow 

real :: fx(maxobs) 
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real :: fy(maxobs) 

real ::origi 

real origj 

real : ; parm 1 , parm2, parm3 

double precision :: olat 
double precision :: olon 

3)c 3|c afc 3^ 4c 3^ aK 4c * * * 4: * * * 3)c 3|c 3^ * * * * * * * * * * * * * * * * * * * Jfc 4: 3|C 

Functions 

4c4c4c4c4c4c4c4c4c4c4c4c4c4c*4c4c4c4c4c***4c4c4c*4c*****4c*4c4c*4:4c*4c*4:*4c*********4c4c4c4c***4c***4:4:4c4:*4c4c 
integer imaxcv, jmaxcv 

CALL GETGEOM(igeom, pijnnm, dsc, ncols, nrows, olat, olon, 

2 origi, origj, dcol, drow, parml, pami2, parm3, istat) 

if (ncols*nrows .gt. ijmax) istat = 2 

if successful getgeom, obtain max 1 st and 2nd dimension of the array 
and convert from lat/lon to i,j and interpolate 

4c4c4c4c4c4c4c4c4c4c4c4c4c4c*4c4c4c4c4c**4c4c4c4c4c4c*4c4c4c4c4c4c4c4c4c*4c**4c*****4c4c***4c********4c*4c4c4c*4c4c4c*4c 

if (istat .eq. 0) then 

imx = imaxcv(ncols, nrows, dsc) 
jmx = jmaxcv(ncols, nrows, dsc) 

CALL VLLX Y (igeom,nobs,obslat,obslon,'d',fx,fy,istat) 
if (istat .eq. 0) then 

CALL FINTRP(fx,fy,nobs,field,iinx,imx,jmx, 0,0., 0.,0.,fob) 

CALL FINTRP(6c,fy,nobs,field,imx,imx,jrnx,l, 

2 bad_value, -999., -999., fob) 

else 

write (*,’("vllxy returns istat =",i5)') istat 
end if 

4c4c4c4c***4c4c*4c****4c4c****4c*****4c**4c4c******4c4c4c4c4c4c*4c4c4c4c4c4c4c4c4c4c4c4c*4c4c*******4c4c*4c4c4c 

if unsuccessful getgeom, print the error msg 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4: 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4: 4c 4c 4c 4c 4c 4: 4c * 4: * 4c 4c 4c 4c 4c 4: ** 4c ** 4c * 4c 4c * 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

else 

write (*,'C'getgeom returns:" 

2 /" pijnam = ",a 

3 /" dsc = ",a 

4 /" ncols =",i5 

5 /" nrows =",i5 

6 /" olat =",f8.2 

7 /" olon =",f8.2 

8 r origi = ",f8.2 

9 r origj =",f8.2 

a r dcol =",f8.2 

b P' drow =",f8.2 

c /" parml =",f8.2 
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d /” pann2 =",f8.2 

e /” pann3 = "^.2 

f /" istat = *',i5)') pijnnm,dsc,ncols,nrows,olat 

g ,olon,origi,origj ,dcol,drow,pann 1 ,parm2,pann3 ,istat 

if (istat .eq. 2) then 

write (*,'("ERROR: Current geometry requires array " 

2 "length ge",ilO,"; ijmax - ’,ilO)’) ncols*nrows,ijmax 

CALL EXIT (1) 
end if 

end if ! (if istat = 0) 
return 

end subroutine fZob 



7. intgeoni.f90 

subroutine mtgeom(geomnainel, geomname2, arr_size, 
2 fieldl, field2, istat) 

C 

C START PROLOGUE 

C 



C sees IDENTIFieATION: @(#)intgeom.f90 1.1 04/24/98 /h/cm/library/mverif/src/sub/intgeom.f90_v 

e 

e eONFIGURATION IDENTIFieATION; NONE 

e 

e MODULE NAME; intgeom 

e 

e DESORIPTION; Interpolates one geometry to another geometry for 
e spherical projection. 

e 

e eOPYRIGHT; (c) 1 998 FLENUMMETOeOEN 

e U.S. GOVERNMENT DOMAIN 

G ALL RIGHTS RESERVED 

G 

G GONTRAGT NUMBER AND TITLE; N/A 
G 

G REFERENGES; NONE 
G 

G GLASSIFIGATION; Unclassified 
G 

G RESTRIGTIONS; NONE 
G 

G GOMPUTER/OPERATING SYSTEM 
G DEPENDENGIES; Gray UNIGOS 

G 

G LIBRARIES OF RESEDENGE; /a/ops/bin 
G 

G USAGE; 

G call intgeom(geomname 1 , geonmame2,arr_size, field 1 , field2, istat) 
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c 

C PARAMETERS: 

C Name Type Usage Description 

C 

C GEOMNAMEl CHAR*32 INPUT GEOMETRY TO INTERPOLATE FROM 

C GEOMNAME2 CHAR*32 INPUT GEOMETRY TO INTERPOLATE TO 

C ARR_SIZE INTEGER INPUT ARRAY_SIZE 
C FIELD 1 REAL(ijmax) INPUT ARRAY TO INTERPOLATE FROM 

C FIELD2 REAL(ijmax) OUTPUT INTERPOLATED ARRAY 

C ISTAT INTEGER OUTPUT STATUS 

C 

C COMMON BLOCKS; N/A 
C 

C FILES: NONE 
C 

C DATABASES: $META_GRK)_DB 
C Name Table Usage Description 
C 

c 

C NON-FILE INPUT/OUTPUT: NONE 
C 

C ERROR CONDITIONS: 

C CONDITION ACTION 

C 

C Error return from GGRD Print err message & exit 
C Error return from GETGEOM Print err message & exit 
C Error return from VXYLL Print err message & exit 
C Error return from CHGEOM Print err message & exit 
C 

C ADDITIONAL COMMENTS: NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED: 

C Name Description 

C 

C CHGEOM CHANGES THE GEOMETRY 
C EXIT System call that exits program 

C GET GEOM Gets geometry arguments to be used by other routines 

C GGRD Returns ISIS info, on given geometry 
C IMAXCV Selects the first dimension of a field 

C JMAXCV Selects the second dimension of a field 

C VXYLL Converts i/j to lat/lon 

C 

C LOCAL VARIABLES AND 
C STRUCTURES: 

C 

C Name Type Description 

C 

C CNCOLS INTEGER column count (getgeom, geomnamel) 

C CNROWS INTEGER row count (getgeom, geomnamel) 

C CORIGLAT DOUBLE latitude of origin (getgeom. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



CORIGLON DOUBLE longtitude of origin (getgeom, geomnamel) 



CORIGX 
CORIGY 
CPARMl 
CPARM2 
CP ARMS 
CPRJNAM 
CSTDESC 
CXINTDIS 
CYINTDIS 
FILVAL 
FLAT 
FLON 
FVALI 
FVALO 
FWl 
FW2 
GEOMI 

GGEOM 



REAL X coordinate of origin (getgeom,geomnamel ) 

REAL y coordinate of origin (getgeom, geomnamel) 

REAL geometry parameter#! (getgeom, geomnamel) 
REAL geometry parameter #2 (getgeom, geomname 1 ) 

REAL geometry parameter #3 (getgeom,geomname 1 ) 

CHAR*24 projection name (getgeom, geomnamel) 
CHAR*24 storage description (getgeom,geomname 1 ) 



REAL interval distance between columns (getgeom,geomnamel) 
REAL interval distance between rows (getgeom,geomname 1 ) 

REAL chgeom parameter 
REAL latitude array 
REAL longtitude array 

REAL chgeom parameter 

REAL chgeom parameter 
REAL work array for chgeom 

REAL work array for chgeom 

INTEGER data structure containing the input 
grid definition (chgeom) 

CHAR*3 2 data structure containing the output 
grid definition (chgeom) 



GMODEL CHAR*32 NOGAPS model 

I INTEGER counter 

BFLAGI INTEGER input field flag (chgeom) 

IMAX INTEGER first dimension of an array 
ISFLG INTEGER chgeom flag 

ISTGR INTEGER stagger flag (chgeom) 

IVEC INTEGER vector flag (chgeom) 

IWRP INTEGER wrap flag (chgeom) 

J INTEGER coimter 

JMAX INTEGER second dimension of an array 

LAFLAG INTEGER land average flag (chgeom) 

LAPASS INTEGER number of passes (chgeom) 
LASRCH INTEGER number of points to search (chgeom) 
LAVAL REAL values in the field (chgeom) 

LEN INTEGER total number of element in an array 
LSTATS INTEGER output unit 
NCOLS INTEGER column coimt 
NGEOM INTEGER data structure containing the input 
grid definition (chgeom) 



NRFCST REAL NORAPS forecast field 
NRFCSTU REAL NORAPS wind_u forecast 
NRFCSTV REAL NORAPS vrind_v forecast 
NROWS INTEGER row coimt 
ORIGLAT DOUBLE latitude of origin 
ORIGLON DOUBLE longitude of origin 
ORIGX REAL x coordinate of origin 

ORIGY REAL y coordinate of origin 

PARMl REAL geometry parameter #1 

PARM2 REAL geometry parameter #2 

PARM3 REAL geometry parameter #3 
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C PRJNAM CHAR*24 projection name 

C STDESC CHAR*24 Storage description 

C XINTDIS REAL Interval distance between columns 

C YINTDIS REAL Interval distance between rows 

C 

C METHOD; 

c 1 . Get ISIS infonnation on geomnamel by calling GGRD and 
c GETGEOM. 

c 2. Call GGRD and GETGEOM on the geomname2. 
c 3. Test for 'spherical' projection. 

c 4. Find the first and second dimension of the geomname2. 
c 5. Convert i/j to lat/lon by calling vxyll. 
c 6. Interpolate field 1 to field2 by calling CHGEOM. 

C 

C INCLUDE FILES: NONE 
C 

C COMPILER DEPENDENCIES: f90 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mverif/src/sub/makeverobslib 
C UNICOS make 
C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1.1 (29 Apr 1998) - Kyongsuk 
C Initial submission 
C 

C END PROLOGUE 

C 

implicit none 
c Formal parameters 

character(32), Intent(in) :: geomnamel ! first geometry name 
character(32), Intent(in) :: geomname2 ! second geom name 
integer, Intent(in) :: arr_size ! array size 
real, Intent(in) : : field 1 (arr_size) ! array to interpolate 
real : : field2(arr_size) ! interpolated array 

integer :: istat ! status 

^%*3|c**a|c:4(:|(:4c:4e4c:»::4c3|:3|c:|c4(4c3)c:|c*3t:*******3t::t:3t::4c**************9t:*4c**4:*3|:*****4c**********3|c 

c Local variables 

^**:|c4c*4c3|:9»:****3t:******* ******** ************************4:*4:*********4c*****:|c* 

character*! uv 
character*24 stdesc, cstdesc 
character*24 prjnam, cpijnam 
character* 32 gmodel, ggeom 

integer i, j 

integer Istats, ngeom 
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integer geomi, ncols, nrows 

integer cncols, cnrows, im, jm, imax, jmax, len 

integer ivec, iwrp, istgr, iflagi, laflag, lasrch, lapass, isflg 

real fVali, fvalo, filval, laval 

real origx, origy, parm 1 , parm2, parm3 , xintdis, yintdis 
real corigx, corigy, cpannl, cparm2, cpaim3, cxintdis, cyintdis 

parameter (im = 360) 
parameter (jm = 181) 

real fwl (im,jm), fw2(imjm), flat(im jm), flon(im,jm) 

double precision origlat, origlon, coriglat, coriglon 
parameter (uv = 'd') ! for vxyll ’d' means in degrees 

^ :|e :|c 3|e :|e :1c *** 3|e ** 3|e *** :|e ** 3|: :|e *** 3^ ***************** **** 9ie * 3|e 9)e * 9|<: 9)e 9|<: 3|<: * 

c Functions 

^********* ********************************************************* ****** 
integer imaxcv, jmaxcv 

^********* *************************************************** ************ 
c get ISIS info on the geometry 1 by calling GGRD 
c and GETGEOM 

^* ************************************************** ********************* 
call GGRD(geomnamel, geomi, istat) 
if (istat .ne. 0) then 

write (0,'("ggrd for geomnamel returns istat =",i5)') istat 
call exit(l 1) 
end if 

call getgeom(geomi, cpijnam, cstdesc, cncols, cnrows, coriglat, 

2 coriglon, corigx, corigy, cxintdis, cyintdis, cpannl, 

3 cparm2, cparm3, istat) 

if (istat .ne. 0) then 

write (0,'("getgeom for geomname returns istat =",i5)') istat 
call exit(l 1) 
end if 

if (cpijnam .ne. 'spherical') then 

write(0,'("projection name for geomnamel is not spherical")’) 
callexit(ll) 
end if 

^******************************************** *************** **:lc ********** 

c get ISIS info on the geometry2 by calling GGRD 
^******************** ********************************************** ****** 
call ggrd(geomname2,ngeom,istat) 
if (istat .ne. 0) then 

write (0,'("ggrd for geomname2 returns istat =",i5)’) istat 
callexit(ll) 
end if 
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call getgeom(ngeom, pijnam, stdesc, ncols, mows, origlat, 

2 origlon, origx, origy, xintdis, yintdis, parm 1 , parm2, 

3 parm3, istat) 

if (istat .ne. 0) then 

write (0,'("getgeom for geomname2 returns istat =”,i5)') 

2 istat 

callexit(ll) 
end if 

if (pijnam .ne. 'spherical') then 

write(0,'("projection name for geomname2 is not spherical")') 
call exit(l 1) 
end if 



c — 

c convert the i,j to lat/lon 
c ref: /a/hbraiy/omsp/chgeom/src/sub/chgeom.f 
c flat=aiTay of latitudes from vxyll for new 2-D array geomo 
c flon=array of longitudes from vxyll for new 2-D array geomo 

C ^ 

imax=imaxcv(ncols,mows,stdesc) 
jmax=jmaxcv(ncols,mows,stdesc) 
doj=l, jmax 
do i=l, imax 
fwl(i,j) = float(i) 
fw2(i,j) = floatO) 
enddo 
enddo 



len = imax*jmax 

call vxyll(ngeom, len, fwl, fw2, uv, flat, flon, istat) 
if (istat .ne. 0) then 
write(0,'("VXYLL error")') 
callexit(ll) 
endif 



c set some of the chgeom parameter values (ref: chgeom write-up) 



ivec = 1 
iwrp = 1 
istgr = 0 
iflagi = 0 
fr^ali = 0.0 
fvolo = 0.0 
laflag = 0 
lasrch = 20 
laval = 0.0 
lapass = 0 
filval =0.0 
isflg = 0 
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c change the field 1 to fieldZ 
c ---- = 



call chgeom(fieldl, field 1, cncols, cnrows, geomi, imax, jmax, 

1 ngeom, ivec, iwrp, istgr, iflagi, fvali, fvalo, 

2 laflag, lasrch, laval, lapass, filval, isflg, 

3 fiv 1 , fw2, field2, field2, istat) 
if (istat .ne. 0) then 

wiite(0,'("CHGEOM error in intgeom")') 
call exit(l 1) 
endif 

return 

end subroutine intgeom 



8 . 



uv2df.f90 



subroutine uv2df(u, v, n, dir, spd) 

START PROLOGUE 

sees BDENTffieATION: @(#)uv2df.f90 1.1 04/24/98 /h/cnvlibraiy/mverifisrc/sub/uv2df.f90_v 
eONFIGURATIONIDENTIFieATION; NONE 
MODULE NAME: uv2df 

DESeRIPTION: This subroutine converts u/v components to a 
field of direction/speed (dd/fi). 

eOP YRIGHT: (c) 1 996 FLENUMMETOeOEN 

U.S. GOVERNMENT DOMAIN 
ALL RIGHTS RESERVED 

eONTRAGT NUMBER AND TITLE: N/A 

REFERENGES: NONE 

GLASSIFIGATION: Unclassified 

RESTRIGTIONS: NONE 

GOMPUTER/OPERATING SYSTEM 

DEPENDENGIES: GrayUNIGOS 

LIBRARIES OF RESIDENGE: /a/opsA)in 

USAGE: 

call uv2df(u, V, n, dir, spd) 

PARAMETERS: 
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Description 



C Name Type Usage 

C 

C U REAL INPUT Wind u-comp fcst interpolated to obs 

C V REAL INPUT Wind v-comp fcst interpolated to obs 

C N INTEGER INPUT No. of observations 

C DIR REAL OUTPUT Converted direction array 

C SPD REAL OUTPUT Converted speed array 

C 

C COMMON BLOCKS; N/A 
C 

C FILES; N/A 
C 

C DATA BASES; N/A 
C 

C NON-FILE INPUT/OUTPUT; N/A 
C 

C ERROR CONDITIONS; N/A 
C 

C ADDITIONAL COMMENTS; NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED; N/A 
C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES; where they are defined in the code 
C within include files. 

C Name Type Description 

C 

C I INTEGER Counter 

C R2D REAL 45.0/atan(1.0) 

C 

C METHOD; 

C Convert u/v to dir/spd by using simple trigonometric functions. 

C 

C INCLUDE FILES; NONE 
C 

C COMPEER DEPENDENCIES; f90 
C 

C COMPILE OPTIONS; -f fixed -c 
C 

C MAKEFILE; Located at /a/ops/app/mverif/src/sub/makeverobshb 
C UNICOS make 
C 

C RECORD OF CHANGES; 

C 

C «CHANGE NOTICE» Version 1 . 1 (29 Apr 1 998) - Kyongsuk Pace 
C Initial submission 
C Original Programmer; M.A. Rennick 
C 

C END PROLOGUE 

C 
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implicit none 



c formal parameters 

real u(n), v(n), dir(n), spd(n) 
integer n 

^ :«:*** :<c :(c *** 3|c 4c :tc ** ** ;fc :(c :fc :<c :fc **** ;tc :fc * :«c ic :|c :|c * :fc 3|c *** :K :(c ;fc :tc *********** jf: ******** * 

c local variables 

^****************************************************************** 
integer i 
real r2d 
real badval 

parameter(badval= 1 .E+ 1 0) 
r2d = 45.0/atan(1.0) 
doi= l,n 

if (u(i) .eq. badval) then 
dir(i) = 999.0 
spd(i) = 999.0 
else 

if (u(i) .eq. 0.0) then 
u(i) = 1 .Oe-6 
end if 

dir(i) = 270.0 - r2d * atan2(v(i), u(i)) 
if (dir(i) .gt. 360.0) then 
dir® = dir® - 360.0 
end if 

spd(i) = sqrt(u(i)*u(i) + v(i)*v(i)) 
end if 
end do 

return 

end subroutine uv2df 

9. Iltread.f90 

subroutine lltread(seq_type, param, Ivl, dsetnm, 

2 dtg, minlat, maxlat, minion, 

3 maxlon, typlvl, obslat, obslon, 

4 nobs, obsval, istat) 

C 

C START PROLOGUE 

C 

C sees E)ENTIFieATION: @(#)lltread.f90 1.1 04/24/98 /h/cm/libraiy/mveri£'src/sub/lltread.f90_v 

e 

e eONFIGURATIONDDENTEFieATION: NONE 
e 

e MODULE NAME: lltread 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



DESCRIPTION; This module calls the appropriate latitude- 
longitude-time (lit) read module based upon 
the user specificed lit sequence type(s). 

There are separate read modules for each sequence 
type because of the different include files and 
data structure in the ISIS for the 
different observation types. 

COPYRIGHT: (c) 1 996 FLENUMMETOCCEN 

U.S. GOVERNMENT DOMAIN 
ALL RIGHTS RESERVED 

CONTRACT NUMBER AND TITLE; N/A 

REFERENCES: ISIS LLT User's Manual 

CLASSIFICATION; Unclassified 

RESTRICTIONS; NONE 

COMPUTER/OPERATING SYSTEM 

DEPENDENCIES: CrayUNICOS 

LIBRARIES OF RESIDENCE; /a/ops/bin 

USAGE; 

call lltread(seq_type, param, Ivl, dsetnm, dtg, 
minlat, maxlat, minion, maxlon, typlvl, 
obslat, obslon, nobs, obsval, istat) 



PARAMETERS: 
Name Type 



Usage 



Description 



SEQ_TYPE 

PARAM 

LVL 

DSETNM 

DTG 

MINLAT 

MAXLAT 

MINLON 

MAXLON 

TYPLVL 

OBSLAT 

OBSLON 

NOBS 

OBSVAL 

ISTAT 



char*24 INPUT LLT obs type 
CHAR*32 INPUT Parameter to read eg. air_temp 
REAL INPUT Pressure level 

CHAR*24 INPUT ISIS dataset name 



CHAR* 10 
REAL 
REAL 
REAL 
REAL 
CHAR*24 



INPUT Date Time Group for read 
INPUT Minimum latitude of the area 
INPUT Maximum latitude of the area 
INPUT Minimum longitude of the area 
INPUT Maximum longitude of the area 



INPUT level type 
REAL(maxobs) OUTPUT Observation latitude 
REAL(maxobs) OUTPUT Observation longitude 
INTEGER OUTPUT No of good obs read 
REAL(maxobs) OUTPUT Observed parameter value 
INTEGER OUTPUT Status 



COMMON BLOCKS: N/A 
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C FILES; NONE 
C 

C DATA BASES: NONE 
C 

C NON-FILE INPUT/OUTPUT: N/A 
C 

C ERROR CONDITIONS; 

C CONDITION ACTION 

C 

C Numobs exceeds maxobs Print err message & return 
C (istat = -l) 

C 

C ADDITIONAL COMMENTS; NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED: 

C Name Description 

C 

C RAOB_QC_READ reads ISIS raob_qc lit data 
C SFC_LAND_READ reads ISIS surface land lit data 
C SFC_SHEP_READ reads ISIS surface ship met lit data 
C SFC_SHIP_MET_QC_READ reads ISIS sfc ship met qc lit data 
C ALTY_READ reads ISIS alty lit data 

C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES: where they are defined in the code 

C within include files. 

C 

C INCLUDE FILES: 

C Name Description 

C 

C v_data.h common variables used for verobs 
C 

C COMPILER DEPENDENCIES; f90 
C 

C COMPILE OPTIONS; -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mverif^src/sub/makeverobslib 
C UNICOS make 
C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1.1 (29 Apr 1998) - Kyongsuk Pace 
C Initial submission 
C 

C END PROLOGUE 

C 

implicit none 
include 'v_data.h' 

^ * 4: * * * * 4= * :4c * 4c 4: 4: 4: :4: 4: 4: 4: 4: 4: :4: 4: :iic 4: 4: 4: 4: :4: 4: 4: 4: 4: 4: :4c 4: :4: 4: 4: :4c 4: 4: 
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c Fonnal parameters 

♦***********♦♦**♦♦******* ********** 



character(24), intent(in) :: seq_type 
Character(32), Intent(in) :: PARAM 1 parameter 
Real, Intent(in) :: Ivl ! level 

CHARACTER(24), I]SriENT(I^ ::DSETNM ! Data set name used. 
CHARACTER(IO), INTENT(IN) ::DTG 1 Date Time Group for read. 
REAL, INTENT(IN) :: MTNLAT ! South latitude boundary. 

REAL, INTENT(IN) :: MAXLAT ! North latitude boundary. 

REAL, INTENT(IN) :: MENLON ! West longitude boundary. 

REAL, INTENT(IN) :: MAXLON 1 East longitude boundary. 
CHARACTER(24), INTENT(IN) :: TYPLVL 1 type level, e.g, isbr_lvl 
Real :: obslat(maxobs) ! obs lat 

Real :: obslon(maxobs) ! obs Ion 

Real :: obsval(maxobs) ! obs data value 

integer :: nobs ! numober of obs 

integer :: istat 



! Local variables used as arguments for LLT read subroutines: 

:|c :|c :)c :|c :)c :|c :|c :)e :|c :)c :fc :|c :)c :)c :|c :|c :|c :|e 3|c :|c :)c :|c ^ :|c :|c :|c ^ :)c :]e ^ :)e :|c :)c :|c :)c 4c :|c :|c :|c :|c :)c :)c :)c :|c :)c :)c :|c :|c :|e :)c :|e :|c :)c :|c :)c :)c :)c :|c :|c :^c :)c :fc :)c :|c :|c :)c :(c 

CHARACTER(8) :: VRSNNAM ! Version of ISIS software used. 

CHARACTER(8) :: SECLVL ! 7 character security classification level. 

REAL :: HR ! Hour cited in the report. 

CHARACTER(16) :: MINDTG ! Minmum date and time group to read. 
CHARACTER(16) :: MAXDTG ! Maximum date and time group to read. 
CHARACTER(16) :: NEW_DTG ! -12 DTG if current not found 
REAL : : MINER 1 Minmum hour to read. 

REAL :: MAXHR 1 Maximum hour to read. 

^ 4c 4c 4e * 4e 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4e 4c 4c 4c 4c 4c 4c 4c 4: 4: 4c 4c 4c 4c 4: 4e 4c 4c 4c 4c 4c 4c 4: 4e 4c 4c 4c 4c 4e 4c 4e 4c 4c 4c 4c 4c 4c 4c 4c 4e 4c 4e 4c 4e 4c 4c 4c 4c 4c 4e 

c Local variables 

^4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4t ♦ *4c 4c 4e 4c 4e 4c4c 4c 4c 4c 4c 4e 4c * ♦ 4c* 4c* ♦ 4c 4c4c 4c4c**4c 4e ** ** * ♦ *♦ ♦* ♦ ♦♦ ♦ ♦ ♦ * ♦* ♦♦ * ♦ ♦ ♦ ♦ ♦ 

integer :: i 

REAL :: PLAT ! Point latitude 

REAL :: PLON ! Point longitude 

REAL :: DISTANCE 1 Radius of circle centered 
! atPLAT/PLON 

CHARACTER(5) :: BLOCK_STATION(10) 

CHARACTER(20) :: STATION_NAME(10) ! International blksta # 

Character(20) :: stn_name 



vrsnnam = 

seclvl = UNCLASS' 

^4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c ************************* * 

c initialize the arrays oblat, oblon, obval 

^4c******* ************************************ **************************** 

doi=l,maxobs 
obslat(i) = bad_value 
obslon(i) = bad_value 
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obsval(i) = bad_value 
end do 



rx+**+**+ ♦♦+**♦**********♦+*******++**+*****♦**********++*+* +*+**♦+♦+♦♦++* 



c RAOB_QC LLT (nogaps, noraps, coamps 






if (seq_type = ’raob_qc’) then 
CALL ^OB_QC_READ(vrsnnam, 

2 dtg, param, Ivl, 

3 minlat, maxlat, minion, 

4 maxlon, obslat, obslon, 

5 nobs, obsval, istat) 



dsetnm, seclvl. 



^:)c3|c:|c*ic:)c:)c*:)c* ************************************************************** 

c sfc_lnd (nogaps, noraps, coamps 

^++ + :<c + :*:* + :<:j|c*j|cj|c ****** ********************* ************************ +***** + :(: 

else if (secLtype = 'sfc_lnd*) then 
CALL SFC_LND_READ(vrsnnam, dsetnm, seclvl, 

2 dtg, param, Ivl, 

3 minlat, maxlat, minion, 

4 maxlon, obslat, obslon, 

5 nobs, obsval, istat) 

^* ******************************************************* **************** 
c sfc_ship (wam_global, otis_global 

^*********** *********************** ************************************** 
else if (seq_type = *sfc_ship') then 
CALL SFC_SHIP_READ(vrsnnam, dsetnm, seclvl, 

2 dtg, param, Ivl, 

3 minlat, maxlat, minion, 

4 maxlon, obslat, obslon, 

5 nobs, obsval, istat) 



^******* ******************************************************** ********* 



c sfc_ship_met_qc (nogaps, noraps, coamps 



^************************************************************************ 



else if (seq_type = ’sfc_ship_met_qc’) then 
CALL SFC_SHIP_MET_QC_READ 
2 (vrsnnam, dsetnm, seclvl. 



3 

4 

5 

6 



dtg, 

minlat, 

maxlon, 

nobs. 



param, 

maxlat, 

obslat, 

obsval. 



Ivl, 

minion, 

obslon, 

istat) 



^************************************************************************ 



c alty (wam_global 



^* ******************** *************************************************** 



else if (seq_type = ’alty’) then 
CALL ALTY_READ( vrsnnam, dsetnm, 

2 dtg, param, Ivl, 

3 minlat, maxlat, minion, 

4 maxlon, obslat, obslon. 



seclvl. 
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5 



nobs, 



obsval, istat) 



c other seq_types 
end if 
return 

end subroutine lltread 

10. alty.f90 

subroutine alty_read(vrsnnam, dsetnm, seclvl, dtg, param. 



2 Ivl, minlat, maxlat, minion, maxlon, 

3 lat. Ion, nobs, obs, istat) 

C 

C START PROLOGUE 

C 



C sees IDENTfflCATION: @(#)alty.f90 1.1 04/24/98 /h/cm/libraiy/mveriRsrc/sub/alty.f90_v 

e 

e eONFIGURATION IDENTIFieATION: NONE 
e 

e MODULE NAME: alty_read 

e 

e DESGRIPTION: subroutine to read the alty data and pick 
e out the obs data for the given parameter 

e 

e GOPYRIGHT: (c) 1 998 FLENUMMETOGGEN 

G U. S. GOVERNMENT DOMAIN 

G ALL RIGHTS RESERVED 

G 

G GONTRAGT NUMBER AND TITLE: N/A 
G 

G REFERENGES: NONE 
G 

G GLASSIFIGATION: Unclassified 
G 

G RESTRIGTIONS: NONE 
G 

G GOMPUTER/OPERATING SYSTEM 
G DEPENDENGIES: GrayUNIGOS 

G 

G LIBRARIES OF RESIDENGE: /a/ops/bin 
G 

G USAGE: 

G call alty_read( vrsnnam, dsetnm, seclvl, dtg, param, 

G Ivl, minlat, maxlat, minion, maxlon, 

G lat. Ion, nobs, obs, istat) 

G 

G PARAMETERS: 

G Name Type Usage Description 
G 
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C VRSNNAM CHAR*8 INPUT lit version name 



c 


DSETNM 


CHAR*24 


INPUT data set name 


c 


SECLVL 


CHAR’^8 


INPUT 


classification 


c 


DTG 


CHAR"' 10 


INPUT 


date time group for read 


c 


PARAM 


CHAR'"32 


INPUT 


parameter 


c 


LVL 


REAL INPUT level type 


c 


MINLAT 


REAL 


INPUT 


minimum latitude 


c 


MAXLAT 


REAL 


INPUT 


maximum latitude 


c 


MINLON 


REAL 


INPUT 


minimum longitude 


c 


MAXLON 


REAL 


INPUT 


maximum longitude 


c 


LAT 


REAL(size) 


OUTPUT 


obs latitude 


c 


LON 


REAL(size) 


OUTPUT 


obs longitude 


c 


NOBS 


INTEGER 


OUTPUT number of obs 


c 


OBS 


REAL(size) 


OUTPUT 


obs value 


c 


ISTAT 


INTEGER 


OUTPUT 


return status 



C 

C COMMON BLOCKS; N/A 
C 

C FILES: None 
C 

C DATA BASES: ISIS LLT_DB 
C Name Table Usage Description 

C 

C alty ALTY IN alty obs 

C 

C NON-FILE INPUT/OUTPUT; N/A 
C 

C ERROR CONDITIONS: 

C CONDITION ACTION 

C - 

C DTG error Print err message & return 

C Error return from LRD Print err message 
C Error return from LCLOS Print err message 
C 

C ADDITIONAL COMMENTS; NONE 
C 



C MAINTENANCE SECTION 

C 

C MODULES CALLED; 

C Name Description 

C LCLOS ISIS LLT close 

C LEN_TRIM Determines the length of a string 
C LRD ISIS LLT read 
C TRIM Removes the trailing blanks 
C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES: where they are defined in the code 

C within include files. 

C 

C METHOD: 

C Set seq_type to 'alty' 
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C See raob_qc_read for the rest. 

C 

C INCLUDE FILES: 

C Name Description 

C 

C ALTY.H alty header file 

C 

C COMPILER DEPENDENCIES: 190 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mveril/src/sub/makeverobslib 
C UNICOS make 

C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1.1 (29 Apr 1998) — KyongsukPace 
C Initial submission 
C 

C END PROLOGUE 

C 

implicit none 
include 'v_data.h' 
include 'ALTY.H' 

^** ******************* ************************************** *********** 
c formal parameters 

^********************************************************************** 
character(8), intent(in) :: vrsnnam 
character(24), intent(in) :: dsetnm 
character(8), intent(in) :: seclvl 
character(lO), intent(in) :: dtg 
character(32), intent(in) :: param 
real, intent(in) :: Ivl 
real, intent(in) :: minlat 

real, intent(in) :: maxlat 

real, intent(in) :: minion 

real, intent(in) :: maxlon 

real :: lat(size) 

real :: lon(size) 

integer ::nobs 

real :: obs(size) 

integer :: istat 

^******************** ♦*****♦* 5K****************************** *********** 

c local variables used as arguments for LRD: 

^***:(c*********:^::J:*:(:*:|c:|c:|c:K************************************************ 

character(24) :: seq_type 
real :: hr 

CHARACTER(16) :: MINDTG ! Minmum date and time group to read. 
CHARACTER( 16):: MAXDT G ! Maximum date and time group to read. 

REAL : : MINHR I Minmum hour to read. 
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REAL : ; MAXHR ! Maximum hour to read. 

CHARACTER(16) :: RSN__IN ! Reporting source name. 

REAL : : FCST_IN ! Desired forecast "TAU". 

CHARACTER(24) : : MINUPTM ! Minimum update time. 

CHARACTER(56) :: REMARKS ! Description of data/assoc, record. 
CHARACTER(16) :: RPT_DTG ! Actual date & time group for report. 

REAL : : RPT_HR ! Reported hour read. 

REAL : : RPT_LAT ! Reported latitude read. 

REAL : : RPT_LON ! Reported longitude read. 

CHARACTER( 16):: RPT_RSN ! Reported longitude read. 

REAL : : RPT_FCST ! Reported "TAU” or forecast time. 

CHARACTER(24) :: RPT_CRETM ! Record creation date. 

CHARACTER(24) :: RPT_UDT ! Report’s last update time. 

INTEGER : : BUFFL AG ! 0 => Input is in FBUFF 

! 1 => Input is in IBUFF 

INTEGER : : LLT_ID ! Unique database LLT identifier for each 
! dataset. 

INTEGER :: BLKSEQID ! Unique database LLT block identifier. 

INTEGER : : RECSEQID ! Unique database LLT record identifier. 

TYPE(alty_int):: EBUFF ! integer record structure. 

TYPE(alty):: FBUFF ! Real record structure. 

c Arguments for LCLOS (that get for values). 

CHARACTER(24) :: SEQTYPE_X ! Report type. 

CHARACTER(8) :: VRSNNAM_X ! Version of ISIS software used. 
CHARACTER(24) :: DSETNAM_X ! Data set name used. 

CHARACTER(8) :: SECLVL_X ! 7 character security 
! classification level. 

CHARACTER(1 6) :: DTG__X ! Date Time Group for write. 

^4: **** :|: :|c :|c :):***** ic :|c :);*************** :fe 4c :^3|c ic 3|c 

c Other local variables 

integer :: levels 
integer :: status, i, status2 

seq_type = 'alty' ! Report type 
istat = 0 



Set up date and time group in YYYYMMDDHH format in DTG. 

IF ( LEN_TRIM(DTG) ==10) THEN 
READ CUNIT=DTG(9:10),FMT='(F2.0)',IOSTAT=STATUS) HR 
IF ( STATUS == 0 ) THEN 
IF(HR<12. )THEN 
HR = 0. 

ELSE 
HR= 12. 
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END IF 
ELSE 

WRITE *, ' Cannot read hour DTG(9: 10), 

2 ' from date & time group TRIM(DTG) 

istat = -1 
RETURN 
END IF 
ELSE 

STATUS = 10 

WRITE *, 'alty: Got date and time group 

2 TRIM(DTG),' of length', 

3 LEN_TRIM(PTG), ' but expected length = 1 0.' 
istat = -1 

RETURN 
END IF 

:fc :4c :4c 4c * 4c 4c * * ♦ * 4c ♦ 4c 4: :4c * * 4c 4c * ♦ ♦ ♦ 4c 4c :4c 4c 3|c 9(; 9|c 9|c :<c :4; 9|c 4: :4e 4: :4e :4e 4c :4c 4; 4c 3|e :4c :4c :4e :«e :4c 

Set the input parameters used to get a read-back value. 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

MINDTG =DTG; MAXDTG =DTG 
MINHR =HR; MAXHR =HR+11.999 
RSN_IN ='*' 

FCST_IN = 0.0 ! Report forecast period or Tau (normal = 0.0) 

MINUPTM = '*• 

BUFFLAG = 0 ! Want (both) floating (and integer). 

4t 4c 4: 4t 4t 4t 4c 4: 4c * 4t 4: 4t 4t 4t 4t 4t 4c 4: 4: 4e 4t 4t 4: 4e 4c 4e 4: 4: 4: 4e 4c 4t 4c 4t 4c 4c 4c 4t 4c 4t 4c 4t 4c 4: 4c 4e 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 

get the data from LLT db 

4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4e 4c 4c 4e 4c 4c 4c 4c 4c 4c 4c 4c 4c 4e 4c 4c 4c 4c 4c 4c 4c 

1=0 
nobs = 0 

DO WHILE ( STATUS = 0 ) 

CALL LRD(seq_1ype, vrsnnam, dsetnm, SECLVL, 

2 MINDTG, MINHR, MAXDTG, MAXHR, 

3 MINLAT, MAXLAT, MINLON, MAXLON, 

4 RSN_IN, FCST_IN, MINUPTM, BUFFLAG, 

5 RPT_DTG, RPT_HR, RPT_LAT, RPT_LON, 

6 RPT_RSN, RPT_FCST, RPT_CRETM, RPT_UDT, 

7 LLT_K), BLKSEQID, RECSEQK), 

8 IBUFF, FBUFF, STATUS ) 

IF ( STATUS /= 0 ) THEN 

IF ( STATUS /= 1 00 ) THEN 
! Ignore normal no-more-data return code 
WRITE *, ' Read from ISIS failed. Code = ', STATUS, 
istat = status 
END IF 

ELSE ! successful LRD 
1 = 1 + 1 



^ 4c 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4e 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4t 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

c pick out the relevant info and fill the array 
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^:ic***4c** ************************ ********************************4c:)c4c*:(c3|c:4e:4e*** 

if (param = ’sig_wav_ht’) then 
if (fbuflf % sig_wav_ht < check_val) then 
nobs = nobs + 1 
lat(nobs) = fbufif % crse_lat 
lon(nobs) = fbuflf % crse_lon 
obs(nobs) = fbuflf % sig_wav_ht 
end if 

else if (param = ’wnd_spd’) then 
if (fbuflf % wnd_spd < check_val) then 
nobs = nobs + 1 
lat(nobs) = fbuflf % crse_lat 
lon(nobs) = fbufif % crse_lon 
obs(nobs) = fbuflf % wnd_spd 
end if 

end if ! param 
end if ! OK status 
END DO ! i loop 

IF ( STATUS = 100 ) STATUS = 0 
write *, ’ ’ 

WRITE *, ’ Called LRD I, ’ tunes.’ 

write *, ’ Read nobs, ’ obs of sequence type TRIM(SEQ_TYPE), 

2 

I :f: :<e 4: 9)c 4c 4: 3|c :f: :f; 3|e 4c 3|c :<e 3|c :<e :4c :i|e 4: 4c :i|e :4c :fe :f: :4c * * * 4: * * * * * * * * * * * * * * * * :1c * 4c :4c :4e ::|e :4c 3)c :|c 3)c 3|c 3)c 3)c ^ 

! Close the dataset (now open for reading) again. 

|4c4c4c4c4c4c4c4c4c***4c4c4c4c4c4c4c4c**********4c4c**4c*******4e****4c*4c**4c4c4c4c4c4c*4c4c*4c4c4c**4c**3)c4c4c 

SEQTYPE_X = ’*’; VRSNNAM__X = ’*’ ; DSETNAM^X = ’*’ 

SECLVL_X = ; DTG__X = 

CALL LCLOS(SEQTYPE_X, VRSNNAM_X, DSETNAM_X, 

2 SECLVL_X, DTG__X, STATUS2) 

IF(STATUS2/=0)THEN 
WRITE *, ' Could not close ISIS table. 

2 ’ Error code is STATUS2, 

END IF 

return 

end subroutine alty_read 

11. raob_qc.f90 

subroutine raob_qc_read(vrsnnam, dsetnm, seclvl, dtg, param, 

2 Ivl, minlat, maxlat, minion, maxlon, 

3 lat. Ion, kfinal, obs, istat) 
c 

c START PROLOGUE 
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sees IDENTIFieATION: @(#)raob_qc.f90 1 . 1 04/24/98 /h/cm/library/mverif/src/sub/raob_qc.f90_v 
eONFIGURATIONIDENTIFieATION; NONE 
MODULE NAME: raob_qc_read 

DESeRIPTION: subroutine to read the raob_qc data and pick 
out the obs data for the given parameter 

eOPYRIGHT; (c) 1 996 FLENUMMETOOGEN 

U.S. GOVERNMENT DOMAIN 
ALL RIGHTS RESERVED 

GONTRAGT NUMBER AND TITLE; N/A 

REFERENGES; LLT User's Manual 

GLASSIFIGATION; Unclassified 

RESTRIGTIONS: NONE 

GOMPUTER/OPERATING SYSTEM 

DEPENDENGES: GrayUNIGOS 

LERARES OF RESrDENGE; /a/ops/bin 

USAGE: 

call raob_qc_read(vrsnnam, dsetnm, seclvl, dtg, param, 
level, minlat, maxlat, minion, maxlon, 
lat. Ion, kfinal, obs, istat) 



PARAMETERS; 
Name Type 



Usage 



Description 



VRSNNAM INTEGER INPUT Version name 
DSETNM GHAR*24 INPUT ISIS dataset name 
SEGLVL GHAR*8 INPUT Security level 

DTG GHAR* 1 0 INPUT Date Time Group for read 

PARAM GHAR*32 INPUT Parameter to read eg. air_temp 

LVL REAL INPUT Pressure level 

MINLAT REAL INPUT Minimum latitude of the area 

MAXLAT REAL INPUT Maximum latitude of the area 

MTNLON REAL INPUT Minimum longitude of the area 

MAXLON REAL INPUT Maximum longitude of the area 

LAT Real(size) OUTPUT obs latitude 

LON Real(size) OUTPUT obs longtitude 

KFINAL INTEGER OUTPUT Number of obs 

OBS Real(size) OUTPUT Observed parameter value 

ISTAT INTEGER OUTPUT Status 

GOMMON BLOGKS: N/A 
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c FILES: None 
c 

c DATA BASES: ISIS LLT database 
c 

c Name Table Usage Description 

c 

c raob_qc RAOB_QC IN raob obs 
c 

c NON-FILE INPUT/OUTPUT: N/A 
c 

c ERROR CONDITIONS: 
c CONDITION ACTION 

c 



c DTG error Print err message & return 

c Error return from LRD Print err message 

c Error return from LCLOS Print err message 
c 

c MAINTENANCE SECTION 

c 

c MODULES CALLED: 
c Name Description 

c 

c LCLOS ISIS LLT close 

c LEN_TRIM Determines the length of a string 
c LRD ISIS LLT read 

c TRIM Removes the trailing blanks 

c 

c LOCAL VARIABLES AND Structures are documented in detail 
c STRUCTURES: where they are defined in the code 

c within include files. 



c 

c METHOD: 

c 1. Set seq_type to Vaob__qc’ 
c 

c 2. Get hr from dtg, set mindtg, maxdtg, minhr, maxhr, rsn_in 
c (reporting source name), fcst__in (desired tau), minuptm 
c (min update time), bufflag (0) 
c 

c 3.Seti=0 
c 

c do while (status = 0) 

c call LRD(seq_type, vrsnnam, dsetnam, seclvl, mindtg, 

c minhr, maxdtg, maxhr, minlat, maxlat, 

c minion, maxlon, rsn^in, fcst_in, minuptm, 

c bufflag, rpt_dtg, rtp_hr, rpt_lat, rpt_lon, 

c ipt__rsn, rpt_fcst, rpt_cretm, rpt__udt, llt_id, 

c blkseqid, recseqid, ibuff, fbuff, status) 

c if status /= 0 then 
c if status /= 100 then 

c write error msg 

c end if 

c else 
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c 
c 
c 
c 
c 
c 

c pick out the relevant info and fill the array 

C 

c nobs = 0 

c levels = fbuflf % prof_cnt 

c do j=l, levels 

Q — 

c pick out only the right level 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c endif (right level) 

c enddo (j loop) 

c 

c endif c if status=0 

c enddo c i loop 
c 

c if (status = 100) then status is good 
c close the dataset by calling LCLOS 
c 

c INCLUDE FILES: 
c Name Description 

c 



if (profit % pres = level) then 

lat(j) = fbufif % crsejat 
lon(j) = fbufif % crsejon 
if pann = 'air_temp’ then 
if profit % air_temp /= missing_value then 
obs(j) = prof_t % air__temp 
nobs = nobs + 1 
endif 

elseif parm = 'geop^ht' then 
if prof_t % geop__ht /= missing_value then 
obs(j) = prof_t % geop_ht 
nobs = nobs + 1 
endif 

elseif parm = 'wnd_dir' then 
if profit % wnd_dir /= missing_value then 
obs(j) = profit % wnd_dir 
nobs = nobs + 1 
endif 

elseif parm = 'wnd_spd’ then 
if prof_t % wnd_spd /= missing_value then 
obs(j) = prof_t % wnd_spd 
nobs = nobs + 1 
endif 
endif 



i = i+l 

if duplicate data using rpt__rsn 
print msg 
endif 
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c RAOB_QC.H raob header file 
c common, inc LLT data stnicture 

c v_data.h common variables for verobs 
c 

c COMPILER DEPENDENCIES; f90 
c 

c COMPILE OPTIONS: -f fixed -c 
c 

c MAKEFILE: Located at/a/ops/app/mverifisrc/sub/makeverobslib 
c UNICOS make 
c 

c RECORD OF CHANGES: 
c 

c «CHANGE NOTICE» Version 1.1 (29 Apr 1998) — Kyongsuk Pace 
c Initial submission 
c 

c END PROLOGUE 

c 

implicit none 

include ’v_data.h’ 
include ’common, inc' 
include ’RAOB^QC.H’ 

:|c:ic3ic:tc :|c3ic3iE:)c:ic :1c :1c :1c :|c :1c :1c 

c formal parameters 

^:ic :|c :|c :|c :}c :1c :|c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :|c :1c :1c :1c :1c :1c :1c :1c * :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :|c :1c :1c :|c :1c :1c :1c :1c :|c :|c :1c :1c :(c :1c :|c :1c :|c :1c :|c :|c :1c :1c :|c 

character(8), intent(in) :: vrsnnam 
character(24), intent(in) :: dsetnm 
character(8), intent(in) :: seclvl 
character(lO), intent(in) :: dtg 
character(32), intent(in) :: param 
real, intent(in) :: Ivl 
real, intent(in) :: mini at 

real, intent(in) :: maxlat 

real, intent(in) : : minion 

real, intent(in) :: maxlon 

real :: lat(size) 

real :: lon(size) 

integer :: kfinal 

real ::obs(size) 

integer :: istat 

:1c 4: :► :1c :1c :1c :1c + + :1c :1c :1c :► :1c :► :1c :1c :1c + + + :1c * + ♦ :1c :1c :1c :1c ^c :1c ♦ :1c :1c * :1c ♦♦♦♦ :1c ♦ :1c :1c ♦ :1c :1c :1c + + :!: :1c :1c 

c local variables used as arguments for LRD: 

******** ****************:)c**********:)c***«*******«**:lc******3lc*******:lc:lc:lc 

character(24) :: seq_type 
real :: hr 

CHARACTER(16) :: MINDTG ! Minmum date and time group to read. 
CHARACTER(16) :: MAXDTG ! Maximum date and time group to read. 

REAL : : MINHR ! Minmum hour to read. 

REAL : : MAXHR ! Maximum hour to read. 
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CHARA.CTER(16) :: RSN_IN ! Reporting source name. 

REAL :: FCST_IN ! Desired forecast "TAU". 

CHARACTER(24) :: MINUPTM ! Minimum update time. 

CHARACTER(56) : : REMARKS ! Description of data/assoc, record. 
CHARACTER(16) :: RPT_DTG ! Actual date & time group for report. 

REAL :: RPT_HR ! Reported hour read. 

REAL :: RPT_LAT ! Reported latitude read. 

REAL :: RPT_LON ! Reported longitude read. 

CHARACTER(16) :: RPT__RSN ! Reported longitude read. 

REAL : : RPT_FCST ! Reported ’TAU" or forecast time. 

CHARACTER(24) :: RPT_CRETM ! Record creation date. 

CHARACTER(24) :: RPT_UDT ! Report's last update time. 

INTEGER :: BUFFLAG ! 0 => Input is in FBUFF 
! 1 => Input is in IBUFF 

INTEGER :: LLT_ID ! Unique database LLT identifier for each 
! dataset. 

INTEGER :: BLKSEQID ! Unique database LLT block identifier. 

INTEGER :: RECSEQID ! Unique database LLT record identifier. 

TYPE(raob_qc_int) :: IBUFF 

TYPE(raob_qc) :: FBUFF ! Real record structure. 

c Arguments for LCLOS (that get for values). 

* :4c :4c 3|c * :)c :4c 9): * :4c ;|c 3)c 3^ 9); :)c * 9): :4: :)c :|c :4c 4: ic it % ;4c % :i|c :4c % 4: 3|c :4c :|c :4c 4: :4c 4c :fc 

CHARACTER(24) :: SEQTYPE_X ! Report type. 

CHARACTER(8) :: VRSNNAM_X ! Version of ISIS software used. 
CHARACTER(24) :: DSETNAM_X ! Data set name used. 

CHARACTER(8) :: SECLVL_X ! 7 character security 
! classification level. 

CHARACTER(16) :: DTG_X ! Date Time Group for write. 

^ 4c * 4: 4c * * * 4c 4: * * * * * 4c * * 4: * 4: 4: 4c * 4c * * * * 4c * 4c 4c 4c 4c 4c 4c * 4c * 4c * * * * * * 4: * 4c 4c 4c 4c 4c 4^ 4c 4c * * 4 e * 4c * 4c 4: * 4: 4: * * * 

c Other local variables 

^4c 4c 4: 4c 4c 4c 4c 4c 9f; 4 e 4c 4c 4c 4c 4c 4c 4c 4: 4; 4c 4c 4c 4c 4c 4: 4c 4: 4; 4c 4; 4c 4c 4c 4 e 4: 4: 4c 4c 4c 4c 4c 4c 4c 4c 4: 4 e 9je 4: 9)C 4c 4c 4: 4e 4c 4c 4c 4c 4c 4c 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

integer :: levels, lvl_l , nmatch, k 
integer :: status, i, j, status2 

seq_type = ’raob_qc' ! obs report type 
istat = 0 



Q 4c 4c 9|c 4c 4c 4c 4c 4c 9|c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 9|c 4c 4c 4c 4c 4c 4c 4c 4c 4c 9|c 4c 4c 4c 4c 4c 4c 4c 4c 4c 9|c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 



c Set up date and time group in YYYYMMDDHH format in DTG. 

Q 4c 4c 4c 4c 4c 4c 4c 4c 9|c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 9|c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

IF ( LEN_TRIM(DTG) == 10 ) THEN 
READ (UNIT=DTG(9.iO),FMT='(F2.0)’,IOSTAT=STATUS) HR 
IF ( STATUS == 0 ) THEN 
IF ( HR < 12. ) THEN 
HR = 0. 

ELSE 
HR= 12. 

END IF 
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ELSE 

WRITE *, ' Cannot read hour DTG(9: 10), 

2 ' from date & time group TRIM(DTG) 

istat = - 1 
RETURN 
END IF 
ELSE 

STATUS = 10 

WRITE *,'raob_qc; Got date and time group TRIM(DTG), 

2 'of length’, 

3 LEN_TRIM(DTG), ' but expected length = 10.' 
istat = -1 

RETURN 
END IF 

^ :{c 4: 4; :|c :|c * :|c 4: :|c :|e :|c :|c :|c :|c :|c % :|c % 4; * * 4: 4: :|c :>c :|c :|c :|c :|c :|c :|c :1c :|c :|c :i)c 4: 4: 4: 4c 4: 4c ;)c % 4: % 4: :|c % :)c % 4; :tc :(c :)c :(c 

c Set the input parameters used to get a read-back value. 

q************************************4c4c******4c* ****♦♦♦♦♦♦♦* *♦♦♦*♦♦♦♦♦* 

MINDTG =DTG; MAXDTG =DTG 
MINHR =HR; MAXHR =HR+1 1.999 
RSN_IN =’*’ 

FCST_IN = 0.0 ! Report forecast period or Tau (normal = 0.0) 

MINUPTM = '*’ 

BUFFLAG = 0 ! Want (both) floating (and integer). 

^ 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c ^ ^ ^ ^ % ’it ^ ^ 3|c 4c 4c 4c 4c % 4c 4e ^ ^ ^ ^ ^ ^ ^ ’ic ^ ^ 4e 4c ^ ’(c ^ ^ ^ ^ % 

c get the data from LLT db 

j,********************************************************************* 
1 = 0 

kfinal = 0 ! overall number of matched pressure level and param 

DO WHILE ( STATUS = 0 ) 

CALL LRD(seq_type, vrsnnam, dsetnm, SECLVL, 

2 MINDTG, MINHR, MAXDTG, MAXHR, 

3 MINLAT, MAXLAT, MINLON, MAXLON, 

4 RSN_IN, FCST_IN, MINUPTM, BUFFLAG, 

5 RPT_DTG, RPT_HR, RPT_LAT, RPT_LON, 

6 RPT_RSN, RPT_FCST, RPT_CRETM, RPT_UDT, 

7 LLT_ID, BLKSEQID, RECSEQID, 

8 IBUFF, FBUFF, STATUS) 
ff (STATUS /=0) THEN 

ff ( STATUS /= 100 ) THEN 
! Ignore normal no-more-data return code 
WRITE *, ' Read from ISIS failed. Code = ', STATUS, '.’ 
istat = -l 
END IF 

ELSE ! successful LRD 
1 = 1+1 



^ 4c 4c 4c 4c 4: 4e 4c 4c 4c 4: 4: 4c 4: 4e 4c 4e 4c 4c 4e 4c 4c 4c 4e 4e 4: 4c 4c 4c 4e 4e 4c 4c 4: 4c 4e 4c 4e 4c 4: 4c 4e 4c 4: 4c 4e 4c 4c 4c 4e 4c 4: 4c 4c 4c 4e 4c 4: 4e 4: 4: 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4e 4c 4c 4c 

c pick out the relevant info and fill the array 
c prof is a sub-structure of raob_qc 
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c prof_cnt has the number of levels 

^ :|c 3|c :4c :(c :4c :|c 4c :(c :i)e :<c 3|c 3|c :|c :4c :4c :<c :|c :<c :|c :i|c He * 3|e * :4: 4c :|c :(c * :(c :)c :(c :|c :(c :|c 4: 3K * ** >1: ******* :f: * ’t: ****** :|c 3|c ;|c ^ 

levels = fbuff % prof__cnt 
do j = 1 , levels 

^******* ************************************************* **:(c:(c:(c:(c************ 

c pick out only the right level 

^*** ************************************************************:(; :(c*4; 4; 4; 4;:^ 4; :^:|c 

nmatch = 0 ! number of matched pressure level and parameter 

lvl_l = int(fbufif % profG) % pres / 100) 
if (lvl_l = Ivl) then 

^******* ********************************************************** ********* 
c we want to use the qc flag to discard the bad obs 

c use the qc flag value of 1 for this obs type 

^********** ********* *********** ********** ********************************** 
if (param = ’air_temp’) then 
if (fbuS%prof(j)%air_temp_qc_id = 1 .and. 

2 fbuff%prof(j)%air_temp < check__val) then 

^******* *************************************************** ******* ********* 
c want to fill the lat, Ion, obs arrays filled from 

c 1 to overall number of matched obs without any 

c skipped indices 

c j index is for the level_cnt for a given lat/lon,* 

c k index is for the matched param and pressure level 

c within the j index 

^************************************************************************** 
nmatch = nmatch + 1 
do k = kfmal+1 , kfinal-riimatch 
lat(k) = fbuff % crsejat 
lon(k) = fbuff % crse_lon 
obs(k) = (fbuff % prof0 % air_temp) 
end do 
end if 

else if (param = 'geop_ht') then 
if (fbufl%prof(j)%geop_ht_qc_id = 1 .and. 

2 fbuff%prof(j)%geop_ht < check_val) then 

nmatch = nmatch + 1 
do k = kfinal+1 , kfmal-riimatch 
lat(k) = fbuff % crsejat 
lon(k) = fbuff % crse_lon 
obs(k) = fbuff % prof(j) % geop_ht 
end do 
end if 

else if (param = ’wnd_dir') then 
if (fbuff%prof(j)%wnd_qc_id = 1 .and. 
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2 



fbufi^oprof(j)%dir < check_val) then 
nmatch = nmatch + 1 
do k = kfmal+1 , kfinal+nmatch 
lat(k) = fbdT % crsejat 
lon(k) = fbuff % crse_lon 
obs(k) = fbuff % prof(j) % dir 
end do 
end if 

else if (param = ’wnd_spd') then 
if (fbufi%prof(j)%wnd_qc_id = 1 .and. 

2 fbufi%prof(j)%spd < check_val) then 

nmatch = nmatch + 1 
do k = kfmal+1 , kfmal+nmatch 
lat(k) = fbuff % crse_lat 
lon(k) = fbuff % crsejon 
obs(k) = fbuff % prof(j) % spd 
end do 
end if 

end if ! param 
end if ! right level 
kfinal = kfinal + nmatch 
end do ! levels loop 
end if ! OK status 
end do ! i loop 
if (status = 100) status = 0 
write *, ' ' 

write *, ' Called LRD I, ' times.' 

write *, ' Read kfinal, ' obs of sequence type 

2 TRIM(seq__type), 

Q :4c :4c % :4c :4c :4c :4c :4c :4c :<e :4c :4c :4c 3|e ic :4c 4e :i|c :4c ie 3|e 3|e 3(c 3|c :4c 3|e :4c ic % 4e 3|c 3(c :i|c :i|e 3|c :4c 3)c :i|c :i|c :4c :4c :i)c 3|c :4c :4c :4c 3(c :<c :4c :4c :4c 3(e :<e :<c 

c Close the dataset (now open for reading) again. 

^ 3|c ^ a|c 3|c :4c 3|c ic 3|c 3|c 9|c :4c * 4c :4c 34: 3|c :4c :4c :4c 3|c :4c 3|c :4c :4c :4c »|c :4c :4c ie * 3|c :4c aic :4c 4c * * 31: 91^ :4c :4 e * 3|c 9|c ic ic :4c :4c :4c 3|c :|c ^ 

SEQTYPE_X = ; VRSNNAM_X = ; DSETNAM_X = 

SECLVL_X = ; DTG_X = 

CALL LCLOS(SEQTYPE_X, VRSNNAM_X, DSETNAM_X, 

2 SECLVL_X, DTG_X STATUS2) 
if(STATUS2/=0)then 
write *, ' Could not close ISIS table. 

2 ' Error code is STATUS2, 

end if 

return 

end subroutine raob_qc_read 

12. sfcland.f90 

subroutine sfc_lnd_read(vrsnnam, dsetnm, seclvl, dtg, param. 
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2 Ivl, minlat, maxlat, minion, maxlon, 

3 lat. Ion, nobs, obs, istat) 

C 

C START PROLOGUE 

C 



C sees IDENTEFieATION: @(#)sfcland.f90 1 . 1 04/24/98 /h/cm/bbraiy/mverif/src/sub/sfcland.f90_v 

e 

e eONFIGURATION IDENTIFieATION: NONE 
G 

e MODULE NAME: sfc_Ind_read 
G 

G DESGRIPTION; subroutine to read the sfc_land data and pick out the 
G obs data for the given parameter 

G 

G GOPYRIGHT: (c) 1996 FLENUMMETOGGEN 

G U.S. GOVERNMENT DOMAIN 

G ALL RIGHTS RESERVED 

G 

G GONTRAGT NUMBER AND TITLE: N/A 
G 

G REFERENGES: NONE 
G 

G GLASSIFIGATION: Unclassified 
G 

G RESTRIGTIONS: NONE 
G 

G GOMPUTER/OPERATING SYSTEM 
G DEPENDENGES: GrayUNIGOS 

G 

G LBRARESOFRESIDENGE: /a/ops/bin 
G 

G USAGE: 

G call sfc_land(vrsnnam, dsetnm, seclvl, dtg, param, 

G Ivl, minlat, maxlat, minion, maxlon, 

G lat. Ion, nobs, obs, istat) 

G 

G PARAMETERS: 



c 

c 


Name 


Type Usage 


Description 


G 


VRSNNAM GHAR*8 


INPUT lit version name 


G 


DSETNM 


GHAR*24 


INPUT data set name 


G 


SEGLVL 


GFIAR*8 


INPUT 


classification 


G 


DTG 


GHAR*10 


INPUT 


date time group for re 


G 


PARAM 


GHAR*32 


INPUT 


parameter 


G 


LVL 


REAL INPUT level type 


G 


MINLAT 


REAL 


INPUT 


minimum latitude 


G 


MAXLAT 


REAL 


INPUT 


maximum latitude 


G 


MINLON 


REAL 


INPUT 


minimum longitude 


G 


MAXLON 


REAL 


INPUT 


maximum longitude 


G 


LAT 


REAL(size) 


OUTPUT 


obs latitude 


G 


LON 


REAL(size) 


OUTPUT 


obs longitude 


G 


NOBS 


INTEGER 


OUTPUT number of obs 
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C OBS REAL(size) OUTPUT obs value 

C ISTAT INTEGER OUTPUT status code 
C 

C COMMON BLOCKS: N/A 
C 

C FILES: None 
C 

C DATABASES: ISISLLT_DB 
C Name Table Usage Description 

C 

C sfc_lnd SFC_LND IN surface land obs 

C 

C NON-FILE INPUT/OUTPUT: N/A 
C 

C ERROR CONDITIONS: 

C CONDITION ACTION 

C 

C DTG error Print err message & return 

C Error return from LRD Print err message 
C Error return from LCLOS Print err message 
C 

C ADDITIONAL COMMENTS; NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED: 

C Name Description 

C 

C LCLOS ISIS LLT close 

C LEN_TRIM Determines the length of a string 

C LRD ISIS LLT read 

C TRIM Removes the trailing blanks 

C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES: where they are defined in the code 

C within include files. 

C 

C METHOD: 

C Set seq_type to 'sfc_lnd' 

C See raob_qc.f90 for the rest. 

C 

C INCLUDE FILES: 

C Name Description 

C 

C SFC_LND.H surface land header file 
C V_DATA.H common variables for verobs 

C 

C COMPILER DEPENDENCIES: f90 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mverif/src/sub/makeverobslib 
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C UNICOS make 
C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1 . 1 (29 Apr 1 998) — Kyongsuk Pace 
C Initial submission 
C 

C END PROLOGUE 

C 

implicit none 
include 'v_data.h' 
include 'SFC_LND.H' 

^ 3)c :)c :4c 3|C 3^ :4c 3f: ** :)C * 4: ************* 3)c :4c 3)c 3)c 3)c :4c 3)c 3)c 3)c :>)c 4c :4c * 3): * 3)C * :4c :i^ }|C :4c 4: ^ 

c fonnal parameters 

^ * * 3)C 34 : 34 : 34 : 3|C 3|C 3)C 3)C * 3)C 3)C 3)C * * 3)C * * 3)C * * * 34c * * * * 3)C * * * * * 34c 3)C 4: 34 : 3}C 34 : 34 : 4c 3)C * 3)C * 3}C 34: * * 34; * * 3|C * 34; 34: 3}C * * * * * * * * * 34 c 3)C * 

character(8), intent(in) :: vrsimam 
character(24), intent(in) :: dsetnm 
character(8), intent(m) :: seclvl 
character(lO), intent(in) :: dtg 
character(32), intent(in) :: param 
real, intent(in) :: Ivl 

real, intent(in) :: minlat 

real, intent(in) :: maxlat 

real, intent(in) :: minion 

real, intent(m) :: maxlon 

real :: lat(size) 

real :: lon(size) 

integer ::nobs 

real :: obs(size) 

integer :: istat 

^ 4c 4c 4c 4c 34 c 4c 4c 4c 4c 4c 4c 34 c 4c 4c 4c 4c *** 4c ****** 4c 4c ** 4c 4c * 4c 4c **** 4c * 4c *** 4c ** 4c 4c 4c 4c 4c **** 4c 4c 4c 4c 4c 4c 4c 4c 34 c *** * 

c local variables used as arguments for LRD: 

^** ******************************************************* ************* 
character(24) :: seq__type 
real :: hr 

CHARACTER(16) :: MINDTG ! Minmum date and time group to read. 
CHARACTER(16) :: MAXDTG ! Maximum date and time group to read. 

REAL : : MINER ! Minmum hour to read. 

REAL : : MAXHR ! Maximum hour to read. 

CHARACTER(16) :: RSN_IN ! Reporting source name. 

REAL :: FCST_IN ! Desired forecast "TAU". 

CHARACTER(24) :: MINUPTM ! Minimum update time. 

CHARACTER(56) :: REMARKS ! Description of data/assoc, record. 

CHARACTER( 16) :: RPT_DTG ! Actual date & time group for report. 

REAL :: RPT_HR ! Reported hour read. 

REAL :: RPT_LAT ! Reported latitude read. 

REAL :: RPT_LON ! Reported longitude read. 

CHARACTER(16) :: RPT_RSN ! Reported longitude read. 

REAL RPT_FCST ! Reported "TAU” or forecast time. 

CHARACTER(24) :: RPT_CRETM ! Record creation date. 
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CHARACTER(24) :: RPT_UDT ! Report's last update time. 

INTEGER : : BUFFLAG ! 0 => Input is in FBUFF 

! 1 => Input is in IBUFF 

INTEGER :: LLT_ID ! Unique database LLT identifier for each 
! dataset. 

INTEGER :: BLKSEQID ! Unique database LLT block identifier. 

INTEGER :: RECSEQE) ! Unique database LLT record identifier. 

TYPE(sfc_lnd_int) IBUFF 

TYPE(sfc_lnd) :: FBUFF 1 Real record structure. 

^ :1c :1c :1c :ic :4c :1c :1c * :1c * :fc :fc :1c :1c :4c :|c :)c :)c :)c :4c :4c :4c :4c :1c :1c :1c :1c :ic :1c :t; :4c :1c :4c :1c :1c * 

c Arguments for LCLOS (that get for values). 

^ :|c 4c :1c :1c :1c :1c :1c 4: :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c * :1c :1c :1c :1c * :1c :1c :1c :1c :1c :1c :1c :1c :1c 4; :1c :1c :1c 4: :1c * :1c :1c :1c :1c :1c :1c % 4: 4c 4c :1c 

CHARACTER(24) :: SEQTYPE_X ! Report type. 

CHARACTER(8) :: VRSNNAM_X ! Version of ISIS software used. 
CHARACTER(24) :: DSETNAM_X ! Data set name used. 

CHARACTER(8) :: SECLVL_X ! 7 character security 
! classification level. 

CHARACTER( 16):: DT G_X ! Date Time Group for write. 

^ 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 

c Other local variables 

^ 4; 4e * 4e :le :le sic 9l( :1c ♦ :1c ♦ ♦♦♦♦♦ )le ** :1c ♦ * 

integer :: levels 
integer :: status, i, status2 

seq__type = 'sfc_lnd' ! Report type 
istat = 0 



4; 4( 4c 4c 4c 4c 4c 4c 4c 4c 4c ♦ 4c 4c 4c 4c 4c 4c 4c ♦ 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c :1c 4c 4c 4c 4c 4c * 4c 4c 4c 4c 4c 4c 4c ♦♦ :1c 4c 4c 4c 4c 4c 4c 4c 4c :1c 4c ♦ 4c 4c 4c ♦ :1c 4c 4c 4c 4c 

Set up date and time group in YYYYMMDDHH format in DTG. 

4e4c4e4e4e4c4c4e4e4e4E4e4E4E4e4c4E4e4e4e4e4E4e4e4e4c4e4c4e4E4e4e4e4e4E4E4e4e4e4e4e4e4E4E4E4E4e4e4e4e4e4e4e4e4e4c4e4e4e4e4e4e4E4e4e4e4c4e4c4E 

ff ( LEN_TRIM(DTG) = 10 ) THEN 
READ (UNIT=DTG(9: 10) JMT='(F2.0)',IOSTAT=STATUS) HR 
ff ( STATUS = 0 ) THEN 
IF ( HR < 12. ) THEN 
HR = 0. 

ELSE 
HR= 12. 

END IF 
ELSE 

WRITE *, ' Cannot read hour DTG(9:10), 

2 ' from date & time group TRIM(DTG) 

istat = -1 
RETURN 
END IF 
ELSE 

STATUS = 10 

WRITE *, 'sfcland; Got date and time group 

2 TRIM(DTG), ' of length 

3 LEN_TRIMGDTG), ' but expected length — 10.' 
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istat = -1 
RETURN 
END IF 

Set the input parameters used to get a read-back value. 

*************************** ******** ****** ******************** 

MINDTG =DTG; MAXDTG =DTG 
MINHR =HR; MAXHR =HR+11.999 
RSN_IN =’*’ 

FCST_IN = 0.0 ! Report forecast period or Tau (normal = 0.0) 

MINUPTM = 

BUFFLAG = 0 ! Want (both) floating (and integer), 

get the data from LLT db 

^^:^f^:^f^f::tf************************************************************** 

1 = 0 
nobs = 0 

DO WHILE ( STATUS = 0 ) 

CALL LRD(seq_type, vrsnnam, dsetnm, SECLVL, 

2 MINDTG, MINHR, MAXDTG, MAXHR, 

3 MINLAT, MAXLAT, MINLON, MAXLON, 

4 RSN_IN, FCST_IN, MINUPTM, BUFFLAG, 

5 RPT_DTG, RPT_HR, RPT_LAT, RPT_LON, 

6 RPT_RSN, RPT_FCST, RPT_CRETM, RPT_UDT, 

7 LLT_ID, BLKSEQID, RECSEQID, 

8 IBUFF, FBUFF, STATUS ) 

IF ( STATUS /= 0 ) THEN 

IF ( STATUS /= 1 00 ) THEN 
! Ignore normal no-more-data return code 
WRITE *, ' Read from ISIS failed. Code = ’, STATUS, V 
istat = status 
END IF 

ELSE ! successful LRD 
1 = 1 + 1 



^* +***********************************★♦**********★****************♦*++**** 
c pick out the relevant info and fill the airay 
c check for ISIS missing value 

+ ******************♦**♦*********♦**♦************************ + * *♦♦♦♦♦♦* 
if (param = 'air_temp') then 
if (fbufiF% air_temp < check_val .and. 

2 fbuff % air_temp_qc_id = 1) then 

nobs = nobs -H 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % air_temp 
end if 

else if (param = ’wnd_dir') then 
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if (fbuff % wnd_dir < check_val .and. 

2 fbuff % vvnd_qc_id = 1 ) then 

nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % wnd_dir 
end if 

else if (param = ’sea_lvl_pres’) then 
if (fbuff % sea_lvl_j)res < check_val .and. 

2 fbuff % sea_lvl_pres_qc_id = 1 ) then 

nobs = nobs + 1 
lat(nobs) = fbuff % crsejat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = (fbuff % sea_lvl__pres) / 100.0 
end if 

else if (param = 'wnd_spd') then 
if (fbuff % wnd_spd < check_val .and. 

2 fbuff % wnd_qc_id = 1 ) then 

nobs = nobs + 1 
lat(nobs) = fbuff % crsejat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % wnd_spd 
end if 

end if ! param 
end if ! OK status 

END DO ! i loop 

IF ( STATUS = 100 ) STATUS = 0 



write ’ ' 

WRITE *, ’ Called LRD I, ’ times.’ 

write *, ’ Read nobs, ' obs of sequence type TRIM(SEQ_TYPE), 

2 

:|c :|c % :fc % :1c :)c :|c :|c :fc ** :)c :|c :(:**** :|c :):**** ******:):% :1c :|c ** ** :|c 

Close the dataset (now open for reading) again. 

SEQTYPE_X = ’*’; VRSNNAM_X = ; DSETNAM_X = ’*' 

SECLVL_X = ; DTG_X = 

CALL LCLOS(SEQTYPE_X, VRSNNAM_X, DSETNAM_X, 

2 SECLVL_X DTG_X. STATUS2) 

IF(STATUS2/=0)THEN 
WRITE *, ' Could not close ISIS table. 

2 ' Error code is STATUS2, 

END IF 

return 

end subroutine sfc_lnd_read 
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13. sfcship.f90 



subroutine sfc_ship_read(vrsnnam, dsetnm, seclvl, dtg, param. 



2 Ivl, minlat, maxlat, minion, maxlon, 

3 lat. Ion, nobs, obs, istat) 

C 

C START PROLOGUE 

C 



C sees IDENTIFieATION: @(#)sfcship.f90 1 . 1 04/24/98 /h/cm/libraiy/mverif7src/sub/sfcship.f90_v 

e 

e eONFIGURATION IDENTIFieATION: NONE 

e 

e MODULE NAME: sfc_ship_read 

e 

e DESeRIPTION: subroutine to read the sfc_ship data and pick 
e out the obs data for the given parameter 

e 

e eOPYRIGHT: (c) 1 996 FLENUMMETOeOEN 

e U.S. GOVERNMENT DOMAIN 

e ALL RIGHTS RESERVED 

e 

e GONTRAGT NUMBER AND TITLE: N/A 

e 

e REFERENGES: NONE 

e 

G GLASSIFIGATION: Unclassified 
G 

G RESTRIGTIONS: NONE 

e 

G GOMPUTER/OPERATING SYSTEM 
G DEPENDENGIES: Gray UNIGOS 

e 

G LIBRARIES OF RESIDENGE: /a/ops/bin 
G 

G USAGE: 

G call sfc_ship_read(vrsnnam, dsetiun, seclvl, dtg, param, 

G Ivl, minlat, maxlat, minion, maxlon, 

G lat. Ion, nobs, obs, istat) 

G 

G PARAMETERS: 

G Name Type Usage Description 

G - 

G VRSNNAM GHAR*8 INPUT lit version name 
G DSETNM GHAR*24 INPUT data set name 
G SEGLVL GHAR*8 INPUT classification 

G DTG GHAR*10 INPUT date time group for read 

G PARAM GHAR*32 INPUT parameter 
G LVL REAL INPUT level type 
G MINLAT REAL INPUT minimum latitude 

G MAXLAT REAL INPUT maximum latitude 

G MINLON REAL INPUT minimum longitude 
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C MAXLON REAL INPUT maximum longitude 
C LAT REAL(size) OUTPUT obs latitude 

C LON REAL(size) OUTPUT obs longitude 

C NOBS INTEGER OUTPUT number of obs 

C OBS REAL(size) OUTPUT obs value 

C ISTAT INTEGER OUTPUT return status 

C 

C COMMON BLOCKS: N/A 
C 

C FILES; None 
C 

C DATA BASES: ISIS LLT_DB 
C Name Table Usage Description 

C 

C sfc_ship SFC_SHIP IN surface ship obs 
C 

C NON-FILE INPUT/OUTPUT; N/A 
C 

C ERROR CONDITIONS; 

C CONDITION ACTION 

C 

C DTG error Print err message & return 

C Error return from LRD Print err message 
C Error return from LCLOS Print err message 
C 

C ADDITIONAL COMMENTS: NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED: 

C Name Description 

C 

C LCLOS ISIS LLT close 
C LEN_TRIM Determines the length of a string 
C LRD ISIS LLT read 

C TRIM Removes the trailing blanks 

C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES; where they are defined in the code 
C within include files. 

C 

C METHOD: 

C Set seq_type to 'sfc_ship' 

C See raob_qc_read for the rest. 

C 

C INCLUDE FILES: 

C Name Description 

C 

C SFC_SHIP.H siuface ship header file 
C V_DATA.H common variables for verobs 
C 

C COMPILER DEPENDENCIES: f90 
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c 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mveril/src/sub/makeverobslib 
C UNICOS make 

C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1.1 (29 Apr 1998) — Kyongsuk Pace 
C Initial submission 
C 

C END PROLOGUE 

C 

implicit none 
include 'v_data.h' 
include 'SFC_SHIP.H' 

c formal parameters 

^3|c :|e 4: :|e 4; :|e 4c 4; s|e :)|e 3|c 9|c :|c 3|e * 4: 3|c 4: * :|c 3|c * * :|c * * * :4: :|c * :|c * * :|c * 4: ^ * 3|c * * * 3|c * * :|c :)c :|e :)c :|e :|e * :|e :|e * :|e :|e :|e :|e :(e :|e :|e ;fc :|e :|e ^ 

character(8), intent(in) ;; vrsnnam 
character(24), intent(in) :: dsetnin 
character(8), intent(in) :: seclvl 
character(lO), intent(in) :: dtg 
character(32), intent(in) param 
real, intent(in) :: Ivl 

real, intent(in) :: minlat 

real, intent (in) :: maxlat 

real, intent(in) :: minion 

real, intent(in) :: maxlon 

real :: lat(size) 

real lon(size) 

integer :: nobs 

real :: obs(size) 

integer :: istat 

^4c4:j|c4:4:4:*j|c4:4:4:44:4j|c4:4:4:4:444:44:4:*4:***4:*4:*****4:*4:**4:***4:4:***************4f*4e*** 

c local variables used as arguments for LRD: 

Q+ + + + + + + + + + + + + + + + + + + + + + + + 4 + + + + + + + + + + + + + 

character(24) :: seq__type 
real :: hr 

CHARACTER( 16):: MINDT G ! Minmum date and time group to read. 
CHARACTER(16) :: MAXDTG ! Maximum date and time group to read 
REAL : : MINHR ! Minmum hour to read. 

REAL :: MAXHR ! Maximum hour to read. 

CHARACTER(16) :: RSN_IN ! Reporting source name. 

REAL :: FCST_IN ! Desired forecast "TAU". 

CHARACTER(24) :: MINUPTM ! Minimum update time. 

CHARACTER(56) :: REMARKS ! Description of data/assoc, record. 
CHARACTER(16) :: RPT_DTG ! Actual date & time group for report. 

REAL :: RPT_HR ! Reported hour read. 

REAL :: RPT_LAT ! Reported latitude read. 
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REAL : : RPT_LON ! Reported longitude read. 

CHARACTER(16) :: RPT_RSN ! Reported longitude read. 

REAL :: RPT_FCST ! Reported "TAU" or forecast time. 

CHARACTER(24) :: RPT_CRETM ! Record creation date. 

CHARACTER(24) :: RPT_UDT ! Report's last update time. 

INTEGER :: BUFFLAG ! 0 => Input is in FBUFF 
! 1 => Input is in IBUFF 

INTEGER :: LLT_ID ! Unique database LLT identifier for each 
! dataset. 

INTEGER :: BLKSEQID ! Unique database LLT block identifier. 

INTEGER :: RECSEQID ! Unique database LLT record identifier. 

TYPE(sfc_ship_int):: IBUFF ! integer record structure. 

TYPE(sfc_ship):: FBUFF ! Real record structure. 

^ :|c :4c * 3(C * 9lc 3|C :|c 3jc :)c :|c 3(C * ^ 3|C 3|C 3)C * * 3(C 3jc :(c * :<c * 3|C * jfc :|c * * * ^ * 3(C 3(C 3|C :(c 3fc ^ :)c * :(c 3(C 3(C 3|c * :(c 3ic 

c Arguments for LCLOS (that get for values). 

^:|c :1c :)c c)c 3|c :(c cfc 3(c 3|c 3)c 3|c :|c cic c)c :(c :)c :fc 3)c 3(c :(c :ic cic 3)c 3(c * :|c cic * :{c 3)c 3|c 3(c * * * :fc :|c * * sic * 3)c 3)c 3)c 3|c * * :{c * c)c 3)c 3)c 3(c 3)c * 3(c :|c 3(c 3|c 3(c :fc :(c :fe 3|c 

CHARACTER(24) :: SEQTYPE_X ! Report type. 

CHARACTER(8) :: VRSNNAM_X ! Version of ISIS software used. 
CHARACTER(24) :: DSETNAIVl_X ! Data set name used. 

CHARACTER(8) :: SECLVL_X ! 7 character security 
! classification level. 

CHARACTER(16) :: DTG_X ! Date Time Group for write. 

^ :1c :1c :1c :1c :1c :fc :|e :1c **** :1c :1c :1c :1c :1c :1c :1c :1c :1c :1c * :1c * :1c ***** :fc :1c :1c ** :1c * :1c ** :1c *:!:*** :1c :1c :fc :fc :1c :|c :1c :1c :|c :|c * :1c * 4c 4; :1c :1c :1c :|c :1c 

c Other local variables 

^******** ***************************************************** ********* 
integer :: levels 
integer :: status, i, status2 

seq_type = 'sfc__ship' ! Report type 
istat = 0 



********************************************************************** 

Set up date and time group in YYYYMMDDHH format in DTG. 
********************************************************************** 

IF ( LEN_TRIM(DTG) = 10 ) THEN 
READ (UNIT=DTG(9:10),FMT='(F2.0)MOSTAT=STATUS) HR 
IF ( STATUS == 0 ) THEN 
IF ( HR < 12. ) THEN 
HR = 0. 

ELSE 
HR= 12. 

END IF 
ELSE 

WRITE *, ' Cannot read hour DTG(9: 10), 

2 ' from date & time group TRIM(DT G) 

istat = -1 
RETURN 
END IF 
ELSE 
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STATUS =10 

WRITE *, 'sfcship: Got date and time group 

2 TRIM(DTG), ' of length 

3 LEN_TRIM(DTG), ' but expected length =10.' 
istat = -1 

RETURN 
END IF 



3|c 3|c * :|c :|c :|c :|c ** 9|C :1c :ic :|c * :|c :|c :|c 3|C :|c 3|C 3lc 3|C * 3|C **************************** 3|c 3|c 3|c 3|c 3|c 3|c 3|C * 3|c 3|c 3|c 

Set the input parameters used to get a read-back value. 
********************************************************************* 

MINDTG =DTG; MAXDTG =DTG 
MINER =HR; MAXHR =HR+1 1.999 
RSN_IN ='*' 

FCST_IN =0.0 ! Report forecast period or Tau (normal = 0.0) 

MINUPTM = 

BUFFLAG = 0 ! Want (both) floating (and integer). 

********************************************************************* 
get the data from LLT db 

********************************************************************* 
1 = 0 
nobs = 0 



DO WHILE ( STATUS = 0 ) 

CALL LRD(seq_type, vrsnnam, dsetnm, SECLVL, 

2 MINDTG. MINER, MAXDTG, MAXHR, 

3 MINLAT, MAXLAT, MINLON, MAXLON, 

4 RSN_IN, FCST_IN, MINUPTM, BUFFLAG, 

5 RPT_DTG, RPT_HR, RPT_LAT, RPT_LON, 

6 RPT_RSN, RPT_FCST, RPT_CRETM, RPT_UDT, 

7 LLT_D, BLKSEQID, RECSEQID, 

8 BUFF, FBUFF, STATUS ) 
ff (STATUS /=0) THEN 

IF ( STATUS /= 1 00 ) THEN 
! Ignore normal no-more-data return code 
WRITE *, ' Read from ISIS failed. Code = ', STATUS, 
istat = status 
END IF 

ELSE ! successful LRD 
1 = 1+1 



^****** ********************************************************* *********** 
c pick out the relevant info and fill the array 
c we want to use the qc flag to discard the bad obs when available 
^**** ************************************************************** ******** 



if (param = 'air_temp') then 
if (fbuff % air_temp < check_val) then 
nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
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obs(nobs) = fbufif % air_temp 
end if 

else if (param = 'sea_lvl_pres') then 
if (fbuff % sea_lvl_pres < check_val) then 
nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = (fbuff % sea_lvl_pres) / 100.0 
end if 

else if (param = 'sea_temp') then 
(check pos_qc_id for position error first 
if (fbuff % pos_qc_id = 0 .or. 

2 fbuff % pos_qc_id = 1 ) then 

lsea_temp_qc flag of 0 or 1 is the only obs we want to use 
if ( (fbuff % sea_temp_qc_id = 1 .or. 

2 fbuff % sea_temp_qc_id = 0) .or. 

3 (fbuff % sea_temp < check_val) ) then 
nobs = nobs + 1 

lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % sea_temp 
end if 
end if 

(for inst_wav_per (grid parm ’peak_wav_per') 

else if (param = 'inst_wav_per') then 
if (fbuff % inst_wav_per < check_val) then 
nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % inst_wav_per 
end if 

(for inst_wav_ht_2 (grid parm 'sig_wav_ht') 

else if (param = 'inst_wav_ht_2') then 
if (fbuff % inst_wav_ht_2 < check_val) then 
nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % inst_wav_ht_2 
end if 

else if (param = ’wnd_dir’) then 
if (fbuff % wnd_dir < check_val) then 
nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % wnd_dir 
end if 
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else if (param = ‘wnd_spd’) then 
if (fbuff % wnd_spd < check_val) then 
nobs = nobs + 1 
lat(nobs) = fbufif % crse_lat 
lon(nobs) = fbufif % crse_lon 
obs(nobs) = fbufif % wnd_spd 
end if 

end if ! param 
end if ! OK status 
END DO ! i loop 

IF ( STATUS = 100 ) STATUS = 0 
write *, ' ' 

WRITE ' Called LRD I, ' times.' 

write *, ’ Read nobs, ' obs of sequence type TRIM(SEQ_TYPE), 

2 

:|c 3|c :|c :|c 3|c 3|c :|c 3|c * 3|c 3|c jf; :|c :4c :(c He :|c : 1 c 3|c 4c :|c :(c :>|c :(c * :(c 3(c :(c :(c :|c 3|c :(c :(c :|c :|c :>|c :|c :>|c 3|c t :>|e 

Close the dataset (now open for reading) again. 

♦ 4c 4c 4c ♦ 4c 4c * He He He He 4c 4c 4c 4c 4c 4e 4c 4c 4c 4c 4c 4e 4c 4c ♦ He ♦♦ 4c 4e 4c 4e 4c 4c ***♦ 4c 4c 4c 4c 4c ♦ 4c He ♦ 4c 4c He 4c 4c He ♦ =Jc ♦♦♦♦♦♦♦♦♦♦♦♦ ♦ 

SEQTYPE_X = ; VRSNNAM_X = ; DSETNAM_X = 

SECLVL_X = ; DTG_X = 

CALL LCLOS(SEQTYPE_X, VRSNNAM_X, DSETNAM_X, 

2 SECLVL_X, DTG_X, STATUS2) 

IF ( STATUS2 /= 0 ) THEN 
WRITE *, ' Could not close ISIS table. 

2 ' Error code is STATUS2, 

END IF 

return 

end subroutine sfc_ship_read 

14. ssmetqc.f90 

subroutine sfc_ship_met_qc_read(vrsnnam, dsetnm, seclvl, dtg. 



2 param, Ivl, minlat, maxlat, 

3 minion, maxlon, lat. Ion, 

4 nobs, obs, istat) 

C 

C START PROLOGUE 

C 



C sees IDENTIFICATION: @(#)ssmetqc.f90 1.1 04/24/98 /h/cm/library/mverif7src/sub/ssmetqc.f90_v 
C 

C CONFIGURATION IDENTIFICATION; NONE 
C 

C MODULE NAME: sfc_ship_met_qc_read 
C 

C DESCRIPTION: subroutine to read the sfc_ship_met_qc data and pick 
C out the obs data for the given parameter 
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c 

C COPYRIGHT; (c) 1 996 FLENUMMETOCCEN 

C U.S. GOVERNMENT DOMAIN 

C ALL RIGHTS RESERVED 

C 

C CONTRACT NUMBER AND TITLE: N/A 
C 

C REFERENCES: NONE 
C 

C CLASSIFICATION: Unclassified 
C 

C RESTRICTIONS; NONE 
C 

C COMPUTER/OPERATING SYSTEM 
C DEPENDENCIES: CrayUNICOS 

C 

C LIBRARIES OF RESIDENCE; /a/ops/bin 
C 

C USAGE: 

C call sfc_ship_met_qc_read(vrsnnam, dsetnm, seclvl, dtg, 

C param, Ivl, minlat, maxlat, 

c minion, maxlon, lat. Ion, 

C nobs, obs, istat) 

C 

C PARAMETERS; 



c 

f 


Name 


Type Usage 


Description 


c 


VRSNNAM CHAR*8 


INPUT lit version name 


c 


DSETNM 


CHAR*24 


INPUT data set name 


c 


SECLVL 


CHAR*8 


INPUT 


classification 


c 


DTG 


CHAR* 10 


INPUT 


date time group for rei 


c 


PARAM 


CHAR*32 


INPUT 


parameter 


c 


LVL 


REAL INPUT level type 


c 


MINLAT 


REAL 


INPUT 


minimum latitude 


c 


MAXLAT 


REAL 


INPUT 


maximum latitude 


c 


MINLON 


REAL 


INPUT 


minimum longitude 


c 


MAXLON 


REAL 


INPUT 


maximum longitude 


c 


LAT 


REAL(size) 


OUTPUT 


obs latitude 


c 


LON 


REAL(size) 


OUTPUT 


obs longitude 


c 


NOBS 


INTEGER 


OUTPUT number of obs 


c 


OBS 


REAL(size) 


OUTPUT 


obs value 


c 


ISTAT 


INTEGER 


OUTPUT return status 



C 

C COMMON BLOCKS: N/A 
C 

C FILES: None 
C 

C DATA BASES: ISIS LLT_DB 
C Name Table Usage Description 

C 

C sfc_ship_met_qc SFC_SHIP_MET_QC IN surface ship met qc obs 

C 
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C NON-FILE INPUT/OUTPUT: N/A 
C 

C ERROR CONDITIONS; 

C CONDITION ACTION 

C 

C DTG error Print err message & return 

C Error return from LRD Print err message 
C Error return from LCLOS Print err message 
C 

C ADDITIONAL COMMENTS: NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED: 

C Name Description 

C 

C LCLOS ISIS LLT close 

C LEN_TRIM Determines the length of a string 
C LRD ISIS LLT read 

C TRIM Removes the trailing blanks 

C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES: where they are defined in the code 

C within include files. 

C 

C METHOD: 

C Set seq_type to 'sfc_ship_met_qc' 

C See raob_qc_read for the rest. 

C 

C INCLUDE FILES; 

C Name Description 

C 

C SFC_SHIP_MET_QC.H surface ship header file 
C V_DATA.H common variables for verobs 
C 

C COMPILER DEPENDENCIES: f90 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/appAnverif/src/sub/makeverobslib 
C UNICOS make 

C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1. 1 (29 Apr 1998) - Kyongsuk Pace 
C Initial submission 
C 

C END PROLOGUE 

C 

imphcit none 

include V_data.h' 

include 'SFC_SHIP_MET_QC.H' 
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c formal parameters 

^ :|e :|c ic * 3|c * 3|c * * 4; :t: :)c 4c * 9)C :4c 9^ * 3|: :)c :|c * ic * 4c :|c * 4c Sf: * * * * * * 4: 9ic :|c * :|c 4: * :|c * * 3|C % :1c ic 4c ic * * 4c * 3|c * 3|c :|C 

character(8), intent(in) :: vrsnnam 
character(24), intent(in) :: dsetnm 
character(8), intent(in) :: seclvl 
character(lO), mtent(in) :: dtg 
character(32), intent(in) :: param 
real, intent(in) :: Ivl 

real, intent(in) :: minlat 

real, intent(in) :: maxlat 

real, intent(in) :: minion 

real, intent(in) :: maxlon 

real :: lat(size) 

real :: lon(size) 

integer :: nobs 

real :: obs(size) 

integer :: istat 



^4c4c4c*4c**4c4c4c4c4c*4c*4c*4c4c***c4c4c4c*4c4c4c*4c4c4c4c4c*4c4c4c4c4c4c4c4c4c**4c4c4c*4c4c*4c4c*4c4c**4c4c4c***4c4c4c 

c local variables used as arguments for LRD: 

character(24) :: seq__type 
real :: hr 

CHARACTER(1 6) : : MINDTG ! Minmum date and time group to read. 
CHARACTER(16) :: MAXDTG ! Maximum date and time group to read. 

REAL : : MINHR ! Minmum hour to read. 

REAL : : MAXHR ! Maximum hour to read. 

CHARACTER(16) :: RSN_IN ! Reporting source name. 

REAL : : FCST_IN ! Desired forecast "TAU". 

CHARACTER(24) :: MINUPTM ! Minimum update time. 

CHARACTER(56) :: REMARKS ! Description of data/assoc, record. 

CHARACTER( 16):: RPT_DT G ! Actual date & time group for report. 

REAL :: RPT_HR ! Reported hour read. 

REAL :: RPT_LAT ! Reported latitude read. 

REAL :: RPT_LON ! Reported longitude read. 

CHARACTER( 16):: RPT_RSN ! Reported longitude read. 

REAL : : RPT_FCST ! Reported "TAU" or forecast time. 

CHARACTER(24) :: RPT_CRETM ! Record creation date. 

CHARACTER(24) :: RPT_UDT ! Report's last update time. 

INTEGER : : BUFFLAG ! 0 => Input is in FBUFF 

! 1 => Input is in IBUff 

INTEGER : : LLT_ED ! Unique database LLT identifier for each 

! dataset. 

INTEGER : : BLKSEQED ! Unique database LLT block identifier. 

INTEGER :: RECSEQID ! Unique database LLT record identifier. 



TYPE(sfc_ship_met_qc_int) :: IBUFF 
TYPE(sfc_ship_met_qc) :: FBUFF ! Real record structure. 

^♦♦♦♦♦♦♦♦4c**********************************4c**************** ♦♦♦♦♦♦♦♦♦♦ 

c Arguments for LCLOS (that get for values). 
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^ 3|C 3|C 3|C3|C* + :)c 3j: * 3j: *:►***:): :|c :►*********:►* :Jc :Jc *** :|c ♦**♦♦*:►****:):**** 3jc ♦ 3|C3|C3JC 3|C ♦ 

CHARACTER(24) :: SEQTYPE_X ! Report type. 

CHARACTER(8) :: VRSNNAM_X ! Version of ISIS software used. 
CHARACTER(24) :: DSETNAM_X ! Data set name used. 

CHARACTER(8) :: SECLVL_X ! 7 character security 
! classification level. 

CHARACTER(16) :: DTG_X ! Date Time Group for write. 

^ 3|C 3|C :4c :)c * :1c ^ :)c :)c :(c * :)c * * * :|c * 3|C 3): 3|C * * 3|C :|c * * * * 3)C 3)C :)c :)c 3|C 3}C * :|c 3)C * * :4c :|c 4: * 3|C 3)C 3|c * 3|C :(c :)c 

c Other local variables 

^*:Jc** + + :|c**:|c*:Jc:Jc**:4c*:*: ********** *:|c*^**:|c:4c**:Jc:|c*:Jc:|c* :►*:►**♦* ♦*:►***:►:► *:(c+:}c:Jc^c:Jc^c^c^c* 

integer :: levels 
integer :: status, i, status2 

seq_1ype = 'sfc_ship_met_qc' ! Report type 
istat = 0 



♦ *♦*♦**♦*♦♦♦♦*♦****♦**********♦♦***♦*♦*♦♦ Jlc**************************** 

Set up date and time group in YYYYMMDDHH format in DIG. 

IF ( LEN_TRIM(DTG) = 10 ) THEN 
READ (UNIT=DTG(9:10),FMT='(F2.0)',IOSTAT=STATUS) HR 
F ( STATUS = 0) THEN 
F ( HR < 12. ) THEN 
HR = 0. 

ELSE 
HR= 12. 

ENDF 

ELSE 

WRITE *, ' Cannot read hour DTG(9: 10), 

2 ’ from date & time group ', TRIM(DTG) 

istat = -1 
RETURN 
ENDF 
ELSE 

STATUS = 10 

WRITE *, 'ssmetqc: Got date and time group ', 

2 TRIMPTG), ' of length ', 

3 LEN_TRIM(DTG), ' but expected length = 10.' 
istat = -1 

RETURN 

ENDF 

****************************************** 3»C3|C*54C54C54c*3»C***3|e*****3|C********* 

Set the input parameters used to get a read-back value. 
********************************************************************* 

MINDTG =DTG; MAXDTG =DTG 
MINHR =HR; MAXHR =HR-:-l 1.999 
RSN_IN =■*' 

FCST_IN = 0.0 ! Report forecast period or Tau (normal = 0.0) 

MINUPTM = 

BUFFLAG = 0 ! Want (both) floating (and integer). 
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get the data from LLT db 

**************************************************l{f*l{f***********iHH^it::i: 
1 = 0 
nobs = 0 

DO WHILE ( STATUS = 0 ) 

CALL LRD(seq_type, vrsnnam, dsetnm, SECLVL, 

2 MINDTG, MINHR, MAXDTG, MAXHR, 

3 MINLAT, MAXLAT, MINLON, MAXLON, 

4 RSN_IN, FCST_IN, MINUPTM, BUFFLAG, 

5 RPT_DTG, RPT_HR, RPT_LAT, RPT_LON, 

6 RPT_RSN, RPT_FCST, RPT_CRETM, RPT_UDT, 

7 LLT_ID, BLKSEQID. RECSEQID, 

8 BUFF, FBUFF, STATUS ) 

IF ( STATUS /=0) THEN 

IF ( STATUS /= 100 ) THEN 
! Ignore normal no-more-data return code 
WRITE *, ' Read from ISIS failed. Code = STATUS, 
istat = status 
END IF 

ELSE ! successful LRD 
I = I-M 



c pick out the relevant info and fill the array 
c we want to use the qc flag to discard the bad obs 

^ ♦ 4: :► :Jc ♦ ♦ * ♦ * it * ♦ * ♦ 4: * * ♦ ie * * ♦ ♦ * ♦ * * ♦ * ♦ * * ♦ * * * ♦ * * 4: ♦ :je * * 4: 4: 4: * ♦ * % 4: :Je 4: :jc 4: % % % :Jc 4: % 4: % * :jc :jc :Jc :|c :Jc 

if (param = ’air_temp') then 
if (fbuff % air__temp < check__val .and. 

2 fbuff % air_temp_qc_id = 1 ) then 

nobs = nobs + 1 
lat(nobs) = fbuff % crsejat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % air_temp 
end if 

else if (param = ’sea_lvl_pres') then 
if (fbuff % sea_lvl_pres < check__val .and. 

2 fbuff % sea_l vl_pres_qc__id = 1 ) then 

nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = (fbuff % sea_lvl_jpres) / 100.0 
end if 

else if (param = ’sea_temp’) then 
if (fbuff % sea__temp < check__val .and. 

2 fbuff % sea_temp_qc_id = 1 ) then 

nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
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obs(nobs) = fbuff % sea_temp 
end if 

else if (param = 'wnd_dir') then 
if (fbuff % wnd_dir < check__val .and. 

2 fbuff % wnd_qc_id = 1 ) then 

nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % wnd_dir 
end if 

else if (param = 'wnd_spd') then 
if (fbuff % wnd_spd < check^val .and. 

2 fbuff % wnd_qc_id = 1 ) then 

nobs = nobs + 1 
lat(nobs) = fbuff % crse_lat 
lon(nobs) = fbuff % crse_lon 
obs(nobs) = fbuff % wnd_spd 
end if 

end if ! param 
end if ! OK status 
END DO ! i loop 

IF ( STATUS = 100 ) STATUS = 0 
write *, ' ' 

WRITE ’ Called LRD I, ’ times.' 

write *, ’ Read nobs, ' obs of sequence type TRIM(SEQ_TYPE), 

2 

4c 9|e :(c :(c :|c 4: ic :(c 4; :4c :4c :4c :4c :4c :(c 4; :(c 

Close the dataset (now open for reading) again. 

4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c * 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4e 4e 4c 4e 4c 4c 4c 

SEQTYPE_X = ’*’; VRSNNAM_X = ; DSETNAM_X = ’*' 

SECLVL_X = ; DTG_X = 

CALL LCLOS(SEQTYPE_X, VRSNNAM_X, DSETNAM_X, 

2 SECLVL_X, DTG_X, STATUS2) 

BF(STATUS2/=0)THEN 
WRITE *, ' Could not close ISIS table. 

2 ' Error code is STATUS2, 

END IF 

return 

end subroutine sfc_ship_met_qc_read 
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B. STAT LIB 



1. Compute bias module 

SUBROUTINE COMPUTE_BIAS (array 1, array2, 

2 arr_size, geoniname, bias) 

C 

C START PROLOGUE 

C 

C sees IDENTIFieATION: @(#)find-bias.f90 1.1 04/24/98 /h/cm/library/mverif/src/sub/find-bias.f90_v 

e 

e eONFIGURATION IDENTIFieATION: NONE 

e 

e MODULE NAME: compute_bias 

e 

e DESeRIPTION: subroutine to compute the bias (mean error) 

e 

e eOPYRIGHT: (c) 1 998 FLENUMMETOeOEN 

e U.S. GOVERNMENT DOMAIN 

G ALL RIGHTS RESERVED 

G 

G GONTRAGT NUMBER AND TITLE: N/A 
G 

G GLASSMGATION: Unclassified 
G 

G RESTRIGTIONS: NONE 
G 

G GOMPUTER/OPERATING SYSTEM 
G DEPENDENGIES: GrayUNIGOS 

G 

G LIBRARIES OF RESIDENGE: /a/ops/bin 
G 

G USAGE: 

G call compute_bias( array 1 , array2, arr_size, geoniname, bias) 

G 

G PARAMETERS: 

G Name Type Usage Description 

G 

G arrayl REAL(360*181) INPUT first array 

G array2 REAL(360*181) INPUT second array 

G arr_size INTEGER INPUT array size 
G geomname GHAR*32 INPUT geometry name 

G bias REAL OUTPUT computed bias 
G 

G GOMMONBLOGKS:N/A 
G 

G FILES: None 
G 

G ERROR GONDITIONS: 

G GONDITION AGTION 
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C 

c 

C ADDITIONAL COMMENTS: NONE 
C 



C MAINTENANCE SECTION 

C 

C MODULES CALLED: 

C Name Description 
C 

C FIND_MAP_FACT determine the map factor for the geometry 
C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES: where they are defined in the code 

C within include files. 



C METHOD: 

C 

C INCLUDE FILES: NONE 
C 

C COMPILER DEPENDENCIES: f90 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mverifisrc/sub/makestatlib 
C UNICOS make 

C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1 . 1 (29 Apr 1 998) — Kyongsuk Pace 
C Initial submission 
C 

C END PROLOGUE 

C 

! formula used: 

! (array 1 - array2) / arr_size 
! [(array 1 - array2) * xmap_factor * ymap_factor] 

! / [4*pi*a_square] 

implicit none 



integer:: im, jm, imjm 
parameter(im = 360) 
parameter(jm =181) 
parameter(imjm = im * jm) 



I :fc3)c:)c**:)c**:ic****:)c***:)c:)e:)c:)c:)c*:)c3)c:tc:)c:)c3)c:)c:)c:)c:)c:)c:)c:ic3)c*:tc:)c:ic:ic*:ic:)c:ic3)c*:ic:ic ********* 

!formal parameter 

|:ic:)c:|c:)c3)c:)c*:)c*:)c:)c:i::)c:)c*:)c:)c:)c:)c*:)c************************************** 



real, intent(in) 
real, intent(in) 
integer, intent(in) 
character(32) 



: : array 1 (imjm) 
:: array2(imjm) 
:: arr_size 
: : geomname 
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real 



bias 



! local var 

ft* **t*******3t:************************3t:3t:**** %*************** 

integer :; i 
integer : : lengeom 
real :: sum, dif, sumx 

real :: xmap_factor(imjm), ymap__factor(imjm) 
real :: pi, a__square 
integer strlen 

pi = 2.0 * asin(l.O) 
a_square = (6.375e-K)6) ** 2 
sum = 0. 
dif=0, 

if (geomname(l :4) /= ’NONE’) then 
write *, ’’calling find-map-factor from find-bias for ”, 

2 geomname 

CALL FIND_MAP__FACTOR(geomname, xmap_factor, ymap_factor) 
do i = 1, arr_size 
dif = (arrayl(i) - array2(i)) 

2 * xmap_factor(i) * ymap__factor(i) 

sum = sum + dif 
end do 

bias = sum / (4*pi*a__square) 
else 

doi = 1, arr_size 
dif = array l(i) - array2(i) 
sum = sum + dif 
end do 

bias = sum / arr_size 
end if 

return 

END SUBROUTINE COMPUTE_BIAS 

2. Compute_rms 

SUBROUTINE COMPUTE_RMS (array!, array2, 

2 arr__size, geomname, rms) 

C 

C START PROLOGUE 

C 

C sees IDENTIFICATION: @(#)find-nns.f90 1 . 1 04/24/98 /h/cm/libraiy/mverifrsrc/sub/find-rms.f90_v 
C 

C CONFIGURATION IDENTIFICATION: NONE 
C 

C MODULE NAME: compute__rms 
C 
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C DESCRIPTION; subroutine to compute the rms 
C 

C COPYRIGHT; (c) 1 998 FLENUMMETOCCEN 

C U.S. GOVERNMENT DOMAIN 

C ALL RIGHTS RESERVED 

C 

C CONTRACT NUMBER AND TITLE; N/A 
C 

C CLASSIFICATION; Unclassified 
C 

C RESTRICTIONS; NONE 
C 

C COMPUTER/OPERATING SYSTEM 
C DEPENDENCIES; CrayUNICOS 

C 

C LIBRARIES OF RESIDENCE; /a/ops/bin 
C 

C USAGE; 

C call compute_rms( array 1 , array2, arr_size, geomname, rms) 

C 

C PARAMETERS; 

C Name Type Usage Description 

C 

C array 1 REAL(360*181) INPUT first array 

C array! REAL(360*181) INPUT second array 

C arr_size INTEGER INPUT array size 
C geomname CHAR*32 INPUT geometry name 

C tins REAL OUTPUT computed rms 
C 

C COMMON BLOCKS; N/A 
C 

C FILES; None 
C 

C ERROR CONDITIONS; 

C CONDITION ACTION 

C 

c 

C ADDITIONAL COMMENTS; NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED; 

C Name Description 

C 

C FIND_MAP_FACT determine the map factor for the geometry 
C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES; where they are defined in the code 

C within include files. 

C METHOD; 

C 

C INCLUDE FILES; NONE 
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c 

C COMPILER DEPENDENCES: 190 
C 

C COMPEE OPTIONS: -f fixed -c 
C 

C MAKEFEE: Located at /a/ops/app/mveriPsrc/sub/makestatlib 
C UNICOS make 

C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1.1 (29 Apr 1998) -- KyongsukPace 
C Initial submission 
C 

C END PROLOGUE 

C 

1 3jc ** 3|S 3|C ****** :|c ic * ic ie :|c 3|e + + + + + ic ♦ 3|e 3js ******* *5|C* *5»: :je 3jc **** * 

!sqrt[(array 1 -array2)**2/arr_size] 

! sqrt [(array 1 -array2) * *2 *xmap_factor*ymap_factor] 

! / [4*pi*a_square] 

I ******************************************************** 

implicit none 

integer:: im, jm, imjm 
parameter(im = 360) 
parameter(jm =181) 
parameter(imjm = im * jm) 

I ******************************************************** 

Iformal parameters 

j ********************************************* *********** 

real, intent(in) :: array 1 (imjm) 

real, intent(in) :: array2(imjm) 

integer, intent(in) :: arr_size 

character(32), intent(in) :: geomname 
real ::rms 

I******************************************************** 

! local var 

j ******************************* ************************* 

integer : : i 

real :: xmap_factor(imjm), ymap_factor(imjm) 
real : : sum, dif, sumx, difx 
real :: pi, a_square 

pi = 2.0 * asin(l.O) 
a_square = (6.375e+06) ** 2 
difx = 0. 
sumx = 0. 

if (geomname(l :4) /= ’NONE') then 
CALL FIND_MAP_FACTOR(geomname, xmap_factor, ymap_factor) 
do i= 1, aiT_size 
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difx = (array l(i) - array2(i))**2 
2 * xmap_factor(i) * ymap_factor(i) 

sumx = sumx + difx 
end do 

rms = sqrt(sumx) / (4*pi*a_square) 
else 

do i = 1, arr_size 
difx = (array 1 (i) - array2(i))**2 
sumx = sumx + dific 
end do 

rms = sqrt(sumx / arr_size) 
end if 

return 

END SUBROUTINE COMPUTE_RMS 

3. Compute_std 

SUBROUTINE COMPUTE_STD (arrayl, array2, 

2 arr_size, geomname, std) 

C 

C START PROLOGUE 

C 

C sees IDENTIFICATION; @(#)fmd-sld.f90 1.1 04/24/98 /h/cm/library/mverif/src/sub/find-std.f90_v 
C 

C CONFIGURATION IDENTinCATION: NONE 
C 

C MODULE NAME: compute_std 
C 

C DESCRIPTION: subroutine to compute the std 
C 

C COPYRIGHT: (c) 1998 FLENUMMETOCCEN 

C U.S. GOVERNMENT DOMAIN 

C ALL RIGHTS RESERVED 

C 

C CONTRACT NUMBER AND TITLE; N/A 
C 

C CLASSIFICATION: Unclassified 
C 

C RESTRICTIONS: NONE 
CC COMPUTER/OPERATING SYSTEM 
C DEPENDENCIES; Cray UNICOS 

C 

C LIBRARIES OF RESIDENCE: /a/ops/bin 
C 

C USAGE: 

C call compute_std( array 1 , array2, arr_size, geomname, std) 

C 

C PARAMETERS: 

C Name Type Usage Description 
C 
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C airayl REAL(360*181) INPUT first array 

C m&y2 REAL(360*181) INPUT second array 

C arr_size INTEGER INPUT array size 

C geomname CHAR* 3 2 INPUT geometry name 
C std REAL OUTPUT computed std 
C 

C COMMON BLOCKS: N/A 
C 

C FILES: None 
C 

C ERROR CONDITIONS: 

C CONDITION ACTION 

C 

c 

C ADDITIONAL COMMENTS: NONE 
C 

C MAINTENANCE SECTION 

C 

C MODULES CALLED: 

C Name Description 
C 

C FIND_MAP_FACT determine the map factor for the geometry 
C 

C LOCAL VARIABLES AND Structures are documented in detail 
C STRUCTURES: where they are defined in the code 

C within include files. 

C METHOD: 

C 

C INCLUDE FILES: NONE 
C 

C COMPILER DEPENDENCIES: 190 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mverif/src/sub/makestatlib 
C UNICOS make 

C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1 . 1 (29 Apr 1 998) — Kyongsuk Pace 



C Initial submission 
C 

C END PROLOGUE. 

C 



!sqrt[{(arrayl-array2)**2/arT_size} - 
! {((array 1 -array2)/arr_size)**2}] 

! sqrt[{ (array 1 -array2) **2 *xmap_factor*ymap_factor } 

! -{((array 1 -array2)*xmap_factor*ymap_factor)**2}] 
! / [4*pi*a_square] 

implicit none 
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integer:: im, jm, imjm 
parameter(im = 360) 
parameter(jm = 181) 
parameter(iinjm = im * jm) 

Iformal parameter 

I********************************************************* 

real, intent(in) : : array 1 (imjm) 

real, intent(in) :: array2(irajm) 

integer, intent(in) :: air_size 
character(32), intent(in) :: geomname 
real :: std 

I********************************************************* 

! local var 

{********************************************************* 
integer :: i 

real :: xmap_factor(imjm), ymap_factor(imjm) 
real :: sum, dif, sumx, difx 
real :: pi, a_square 

pi = 2.0 * asin(l.O) 

a_square = (6.375e-K)6) ** 2 

sum = 0. 

dif=0. 

dific = 0. 

sumx = 0. 

if (geomname(l :4) /= ’NONE’) then 
CALL FIND_MAP_FACTOR(geomname, xmap_factor, ymap_factor) 
do i = 1 , arr_si 2 e 
dif = (array l(i) - array2(i)) 

2 * xmap_factor(i) * ymap_factor(i) 

dific = (array l(i) - array2(i))**2 
2 * xmap_factor(i) * ymap_factor(i) 

sum = sum + dif 
sumx = sumx + difx 
end do 

std = sqrt(sumx / 4*pi*a_square - (sum / (4*pi*a_square))**2) 
else 

do i = 1, arr_size 
dif = arrayl(i) - array2(i) 
sum = sum + dif 
dific = (array 1 (i) - array2(i))**2 
sumx = sumx + dific 
end do 

std = sqrt(sumx/arr_size - (sum/arr_size)**2) 
end if 

return 
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END SUBROUTINE COMPUTE STD 



C. GRAPHICS 



1. Plot_data.pro 

pro plot_data 

; setup color chart 

.^:3|e:|c3|e3|c3|c3|e:)c3|c3|e:|c:)c3|c:)c3te:jc:|E4::t:3|c3|c3|c:4::t:3|c3|c3|e3|c3|c:|c:|c:|c3|c3|c:4:4E3|c)|c3|c4:3|c:|c3|c3|c:|c3|c:|c:ic3|c3|c3|c:f;:|c4:4:3|c:|c:|c^ 

J 

COMMON colors, r_orig, g__orig, b_orig, r_curr, g_curr, b_curr 

maxcol = !D.N_COLORS 

r_cuTT = bindgen(maxcol) 

g__curr = r_curr 

b_curr = r_curr 

wht,gm,gry,wht,tur,blu,ylw,pur,blk,red 
r_curr=[255, 0,211,255, 0, 0,255,55, 0,255] 
g_curr= [255,100,211,255,200, 0,255, 0, 0, 0] 
b_curr=[255, 0,211,255.230,255, 0,55, 0, 0] 

.*5|C*S|C5|£**S|£*3|C*i|C5|£** + 3|c4:** **♦♦♦****♦*♦♦♦**♦♦*♦* ♦♦*♦*♦*♦*♦♦*♦ 

; to graph on the screen, must have the DISPLAY env set 

.*******i|(:* + **j|tJ|£s|c + * + ♦♦ + + **♦ + * + *♦ + ***♦ + *:►♦ + + ****:►**:►****♦ 

;TVLCT, r_curr, g_curr, b_curr 

;integers 
nColors = 0 
nHeader = 0 
nRecords = 0 
nmatch = 0 
fcstPer = 0 

;floats 
sng = 0.0 
level 1 = 0.0 
minrange = 0.0 
maxrange = 0,0 

; strings 
filename = " 
str = ” 

strFormat = " 
strHeader = ” 
strLegend = " 
param = " 
obType = " 
geomName = ” 
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modelName = 
name = " 
yname = " 



;structures 

datLine = {dat, v_dtg;' nobs:0, parm:' units:' $ 
geom:' typlvl:’ lvl_l :0.0, tau:0, $ 
typstat:' stat_val:0.0, v_src:' obs__type:' '} 

strFmt = '(alO, 2x, i5, a20, al5, a30, al5, f8.2, i5, al5, f8.2, alO, a25)’ 

; get the env vars and data filename 

param = GETENV('PARM_NAME') 
levell =GETENV('LVL_r) 
fcstPer = GETENVCPCSTPER') 
obType = GETENV('OBST YPE') 
geomName = GETENV('GEOM_NAME') 
modelName = GETENVCMODEL') 
fflename = GETENV('FILENAME') 

; determine the number of records in the data file 
data=READ_ASCn(fileName, count=nRecords) 
print, "record count = ", nRecords 

;array declarations 

fields = replicate(datLine, nRecords) 

bias = replicate(datLine, nRecords) 

std = replicate(datLine, nRecords) 

rms = replicate(datLine, nRecords) 

stdl = FLTARR(nRecords) 

std2 = FLTARR(nRecords) 

strDTG = STRARR(nRecords) 

lonDTG = FLTARR(nRecords, /NOZERO) 

numObs = INTARR(nRecords) 

openr, 1 0, filename 

for n=0, nRecords- 1 do begin 
READF, 10, datLine, FORMAT=strFmt 
fields[n] = datLine 
endfor 

close, 10 

fields = fields(SORT(fields[*].v_dtg)) 

y 

; get the sub arrays 

.:|c***********:|c******:4c ***4:***** ************************* 

5 

bias = fields[WHERE((STRTRIM(fields[*].pann) EQ param) and $ 
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(fields[*].lvl_l EQ level 1) and $ 

(STRTRIM(fields[*].geom) EQ geomName) and $ 

(fields[*].tau EQ fcstPer) and $ 

(STRTRJM(fieIds[*].obs_type) EQ obType) and $ 
(STRTRIM(fields[*].typstat) EQ ’bias’), rnnatch)] 
start_date = JULDAY(STRMID(bias[0].v_dtg,4,2), $ 

STRMID(bias[0].v_dtg,6,2), $ 

STRMID(bias[0] . v_dtg,0,4)) 

end_date = JULDAY(STRMID(bias[nmatch-l].v_dtg,4,2), $ 
STRMID(bias[nmatch-l].v_dtg,6,2), $ 

STRMID(bias[nmatch-l ].v_dtg,0,4)) 
start_time = FIX(STRMID(bias[0].v_dtg,8,2)) / 24. 
end__time = (end_date - start_date) $ 

+ (PLOAT(STRMID(bias[nmatch-l].v_dtg,8,2))) / 24. 

for n=0, nmatch-1 do begin 
end_datel = JUEDAY(STRMID(bias[n].v_dtg,4,2), $ 

STRMID(bias[n].v_dtg,6,2), $ 

STRMID(bias[n] . v_dtg,0,4)) 
end_timel = (end_datel - start_date) $ 

+ (FLOAT(STRMID(bias[n].v^dtg,8,2))) / 24. 
lonDTG[n] = (end_timel - start_time) + stait_time 
numObs[n] = bias[n].nobs 
endfor 

std = fields[WHERE((STRTRIM(fields[*].parm) EQ param) and $ 

(fields[*].lvl_l EQ levell) and$ 

(STRTRIM(fields[*].geom) EQ geomName) and $ 

(fields[*].tau EQ fcstPer) and $ 

(STRTRIM(fields[*].obs_type) EQ obType) and $ 
(STRTRIM(fields[*].typstat) EQ ’std’), nmatch)] 
rms = fields[WHERE((STRTRIM(fields[*].pann) EQ param) and $ 

(fields[*].lvl_l EQ levell) and $ 

(STRTRIM(fields[*].geom) EQ geomName) and $ 

(fields[*].tau EQ fcstPer) and $ 

(STRTRIM(fields[*].obs_type) EQ obType) and $ 
(STRTRIM(flelds[*].typstat) EQ ’rms’), rnnatch)] 

for n=0, rnnatch- 1 do begin 
stdl[n] = bias[n].stat_val + std[n].stat__val 
std2[n] = bias[n].stat_val - std[n].stat_val 
endfor 

dummy = LABEL_D ATE(DATE_FORMAT = '%HZ %D%M %Z’, oflfset=start_date) 
i = nmatch - 1 

• ^:|C3fC3fC3fC3|c4E3fC3|C3|e3fC3|e3{C3jC3|C3fC3fC3fC3ie3jC3iC3{e3fC3fC3fc:)C3ie3|e3je3jC3iC3{C3fC3fC3{C3fC3fC3|C3fC3fc:|:3|C3iC3{c4!4^3iC3|e3ic3)c:)c4!4^9tc 
9 

; plot the data 

.:)C:]e3|C3fC3|C3|C3|C3(C3|;3|C3(C3|C3fC3fC3|C3|C3fC3|C3|C3iC3fC3|C3fC3ie3|e3|C3|C3|C3|C3|C3fC3fC3|C3fC3fC3(C3fC3|C3|C3|C:|:3|C3|C:)C3|C:)o|C3|C:)c:|C3fc4^:)c:|C 

J 

; title string 

• ♦ ♦ ♦ 3|e * 3|c 3|c * ♦ * * 3|c :|c 3|c 3|c 3|c sf: * ♦ ♦ ♦ 3|e sj: sj: s(: 3|c :J: s|c 3|c 3|e * 3|c 3je + 3|c + :1s sf: 3|c ♦ :|c sj: 3|c sj: s|: ♦ :|c + + + 

9 
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if STRTRIM(bias[0].typlvl) EQ 'isbrjvl' then $ 
name = modelName + + geomName + + param + '!C* $ 

+ STRTRIM(STRJNGCbias[0].lvl_l)) $ 

+ 'mb, ’ + STRTRIM(STRING(fcstPer)) + ' hrs, ' + obType $ 

else $ 

name = modelName + + geomName + + param + '!C $ 

+ STRITUM(STRING(bias[0].lvl_l)) $ 

+ 'm, ' + STRTRIM(STRING(fcstPer)) + ' hrs, ' + obType 

yname = bias[i]. units ; y-axis title string 
minrange = MIN(bias[0;i].stat_val) 
if (MIN(std2[0:i]) LT minrange) then $ 
minrange = MIN(std2[0:i]) 
if (MIN(rms[0:i].stat_val) LT minrange) then $ 
minrange = MIN(rms[0:i].stat_val) 
maxrange = MAX(bias[0:i].stat_val) 
if (MAX(stdl [0:i]) GT maxrange) then $ 
maxrange = MAX(stdl [0:i]) 
if (MAX(rms[0:i].stat_val) GT maxrange) then $ 
maxrange = MAX(rms[0:i].stat_val) 

.****^********* ******** ***** * 5 ):** * + *:►*** 

j 

; for debugging 

.****************************************************** 
;print, "std" 

;print, std[0:i].stat_val 
;print, "rms" 

;print, rms[0:i].stat_val 

;print, "graph range lies between ", minrange, " and ", maxrange 
;print, "lonDtg" 

;print, lonDTG[0:i] 

;print, "bias" 

;print, bias[0:i] 

;print, "stdl " 

;print, stdl[0:i] 

;print, "std2" 

;print, std2[0:i] 



lX.MINOR = - 1 ;suppress minor tick marks 

!Y.MARGIN(1) = 3 ;top margin 

.****************************************************** 

j 

; to create a post script file 

******************************************************* 

j 

;set_plot, 'PS' 

;psfile = filename + '.ps' 

;device, /color, filename=psfile 

******************************************************* 
; to create a gif file 

******************************************************* 
set_plot, 'Z' 
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psfile = filename + '.gif 



PLOT, lonDTG[0;i], bias[0:i].stat_val, $ 

YRANGE = [minrange, maxrange], $ 

TITLE = name, PSYM = -2, SYMSIZE = 1., $ 

XTITLE = 'Forecast Date', $ 

YTITLE = yname, $ 

XGRIDSTYLE = 1, YGRIDSTYLE = 1, $ 

XTICKLEN = 1 .0, YTICKLEN =1.0,$ 

XTICKFORMAT = "label_date", $ 

XCHARSIZE = 0.7, $ 

MAX_VALUE = 30, $ 

COLOR = 5, XSTYLE = 2, /DEVICE, /NODATA 
for n=l,i do begin 

POLYFILL, [lonDTG(n-l), lonDTG[n-l:n], lonDTG(n)], $ 

[bias[n-l].stat_val, stdl [n-I :n], bias[n].stat_val], COLOR=2 

endfor 

for n=l ,i do begin 

POLYFILL, [lonDTG(n-I), lonDTG[n-l :n], IonDTG(n)], $ 

[bias[n-l].stat_val, std2[n-l :n], bias[n].stat_val], COLOR=2 

endfor 

; overplot the bias 

OPLOT, lonDTG[0:i], bias[0;i].stat_val, PSYM = -2, SYMSIZE = 1., S 
LINE = 0, COLOR = 5, MAX_ VALUE = 30 

, overplot the stdl 

OPLOT, lonDTG[0;i], stdl [0:i], PSYM = -6, SYMSIZE = 1 ., $ 

LINE = 0, COLOR = 9, MAX_ VALUE = 30 

; overplot the std2 

OPLOT, IonDTG[0:i], std2[0;i], PSYM = -6, SYMSIZE = 1., S 
LINE = 0, COLOR = 9, MAX_ VALUE = 30 

; overplot the ims 

OPLOT, lonDTG[0:i], ims[0:i].stat_val, PSYM = -4, SYMSIZE = 1., $ 
LINE = 0, COLOR = 7, MAX_ VALUE = 30 

; add the legends 

XYOUTS, 0.8, 0.16, '!5* - bias', color = 5, /NORMAL ; legend 
XYOUTS, 0.8, 0. 13, '!MB - std', color = 9, /NORMAL ;legend 
XYOUTS, 0.8, 0.1, '!MV - rms', color = 7, /NORMAL ;legend 

image = TVRDQ 
WRITE_GIF, psfile, image 

T)EVICE, /close_file 
SET_PLOT, 'X' 

END 
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D. 



USER INTERFACE 



1. Index.html 



<html> 

<!- 

Author: Susie Pace 

Date: 25 March 1998 

File URL: model_reports/mveri^index.htnil 
— > 

<headxtitle>FNMOC Model Statistics Display </title></head> 

<body bgcolor="#l 91 970" TEXT="#F5F5DC" LINK="#00FF7F" 
VLINK="#CCCC66" ALINK="#FF0000"> 

<FONT SIZE=+1> 

<hlxcenter>Pick a model to see the statistics.</center></hl> 

<UL> 

<A HREF="nogaps.html"> NOGAPS</A></P> 

<AHREF="noraps_asia.html">NORAPS_ASIA</A></P> 

<AHREF="noraps__conus.html">NORAPS_CONUS</A></P> 

<AHREF="noraps__europe.html">NORAPS__EUROPE</A></P> 

<AHREF="noraps_ind__ocn.html">NORAPS_IND_OCN</Ax/P> 

<A HREF ="coamps_europe.html ">CO AMP S_EUROPE</Ax/P> 
<AHREF="coamps_swasia.html">COAMPS_SOUTHWEST_ASIA</Ax/F^ 
<A HREF="wam.html">WAM_GLOBAL</A></P> 

<AJL> 

<HR> 

<CENTER><FONT size="-l "><I> Send Comments Or Suggestions To 
Susie Pace: <A HREF="mailto:pacek@fiimoc.navy.mil"> 
pacek@ftimoc.navy.mil</A> 

<BR> Last Update Was On March 25, 1 998 </I></FONT></CENTER> 

</BODY> 

</HTML> 

2. Nogaps.html 



<html> 

<!- 

Author: Susie Pace 

Date: 25 March 1998 

File URL: model_reports/mverif/nogaps.html 
— > 



<headxtitle>NOGAPS Verification Display </title><yhead> 

<body bgcolor="#191970" TEXT="#00ff7f LINK="#00FF7F" VLINK="#CCCC66" 
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ALINK="#FFOOOO"> 

<form method=GET action="http ://devu 1 /cgi-bin/space3 .pl"> 

<hl><center>Make your selections to see the statistics. </center></hl> 

<h3 >Model : </h3 > 

<input type=radio name=model value="nogaps" checked>nogaps 
<h3>Geometry:</h3> 

<input type=radio name=geometry value="global_360xl 8 1 " checked>global_360xl 8 1 
<input type=radio name=geometry value="asia_nestl_appl"> asia_nestl_appl 
<input type=radio name=geometiy value="conus_nestl_appr’> conus_nestl_appl 
<input type=radio name=geometry value="europe_nestl_appr’>europe_nestl_appl 
<input type=radio name=geometiy value="europe_nest2_appl2"> europe_nest2_appl2 
<input type=radio name=geometry value="europe_nest3_appl3*'>europe_nesG_appl3 
<input type=radio name=geometiy value="ind_ocn_nestl_appr'>ind_ocn_nestl_appl 
<!— <input type=radio name=geometry value="southwest_asia_nest2_appr’>southwest_asia_nest2_appl 
<input type=radio name=geometry value="southwest_asia_nest3_appl">southwest_asia_nest3_appl 
— > 



<h3>Parameters:</h3> 

<input type=radio name=parameter value="air_temp"> air_temp 
<input type=radio name=parameter value="geop_hf checked> geop_ht 
<input type=radio name=parameter value="pres"> pres 
<input type=radio name=parameter value="wnd_spd”> wnd_spd 



<h3>Taus:</h3> 



<input type=radio name=tau value="0"> 0 
<input type=radio name=tau value="12"> 12 
<input type=radio name=tau value="24" checked> 24 
<input type=radio name=tau value="36"> 36 
<input type=radio name=tau value="48"> 48 
<input type=radio name=tau value="60"> 60 
<input type=radio name=tau value="72"> 72 
<input type=radio name=tau value="84"> 84 
<input type=radio name=tau value="96”> 96 
<input type=radio name=tau value="108"> 108 
<input type=radio name=tau value=" 1 20"> 1 20 
<input type=radio name=tau value="l 32"> 1 32 
<input type=radio name=tau value=" 1 44 "> 1 44 



<h3>Levels: </h3> 

<input type=radio name=level value="0">0 
<input type=radio name=level value=”2">2 
<input type=radio name=level value=" 1 9.5">1 9.5 
<input type=radio naine=level value=” 1 000 "> 1 000 
<input type=radio name=level value="925">925 
<input type=radio name=level value="850">850 
<input type=radio name=level value="700”>700 
<input type=radio name=level value="500" checked>500 
<input type=radio name=level value="400">400 
<input type=radio name=level value="300">300 
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<input type=radio name=level value="250''-"*250 
<input type=radio name=level value="200">200 
<input type=radio name=level value=" 1 50">1 50 
<input type=radio name=level value="100”>100 

<!~<h3>Statistics:</H3> 
bias: <input type=checkbox name="bias"> 

stdev: <input type=checkbox name="stdev"> 

rms: <input type=checkbox name=”nns"> 

— > 



<h3>Obs types:</h3> 

<input type=radio name=obstype value="raob_qc" checked>raob_qc 

<input type=radio name=obstype value="sfc_lnd">sfc_lnd 

<input type=radio name=obstype value='’sfc_ship_met_qc">sfc_ship_met_qc 

<h3>Graph Type:</h3> 

<!— <input type=radio name=graph value="scatter plot">scatter_plot--> 
<input type=radio name=graph value="time series" checked>time_series 

<h3>Period: Enter the beginning and ending DIG: e.g., 1 99803 1 81 2</h3> 
<input type=text name=beginning maxlength= 1 0 value=" 1 9980325 1 2"> 
<input type=text name=ending maxlength=l 0 value=" 1 998040500"> 

<input type=submit> 

<input type=reset value="Cancer’> 

</form> 

</bodp^ 

<yhtml> 

3. Noraps_asia.html 

<html> 

<!- 

Author; Susie Pace 
Date: 07 November, 1 997 

File URL: model_reports/mveriL^noraps_asia.html 
— > 



<head><title>NORAPS_ASIA Verification Display </title></head> 

<body bgcolor="#191970" TEXT="#00FF7F" LINK="#00FF7F" VLINK="#CCCC66" 
ALINK="#FF0000"> 

<form method=GET action="http://devul/cgi-bin/space3.pl"> 

<hlxcenter>Make your selections to see the statistics. </center></hl> 

<h3>Model :</h3> 

<input type=radio name=model value="noraps_asia" checked>noraps_asia 
<h3>Geometry:</h3> 

<input type=radio name=geometry value="asia__nestl_appl" checked> asia_nestl_appl 
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<h3>Parameters:</h3> 

<input type=radio name=parameter value="air_temp"> air_temp 
<input type=Tadio name=parameter value="geop_ht” checked> geop_ht 
<input type=radio name=parameter value="pres"> pres 
<input type=radio name=parameter value="wnd_spd"> wnd_spd 

<h3>Taus:</h3> 

<input type=radio name=tau value="0"> 0 
<inputtype=Tadioname==tau value- '12"> 12 
<input type=radio name=tau value="24" checked> 24 
<input type=radio name=tau value="36"> 36 
<input type=radio name=tau value="48"> 48 

<h3>Levels: </h3> 

<input type=radio name=level value="0">0 
<input type=radio name=level value="2">2 
<input type=radio name=level value=" 1 9.5">1 9.5 
<input type=radio name=level value=” 1000">1000 
<input type=radio name=level value="925">925 
<input type=radio name=level value="850">850 
<input type=radio name=level value='700">700 
<input type=radio name=level value="500” checked>500 
<input type=radio name=level value="400”>400 
<input type=radio name=level value=”300”>300 
<input type=radio name=level value="250">250 
<input type=radio name=level value="200”>200 
<input type=radio name=level value=" 1 50">1 50 
<input type=radio name=level value="100">100 

<!— <h3>Statistics:</H3> 
bias: <input type=checkbox name="bias"> 

stdev: <input type=checkbox name="stdev''> 

rms: <input type=checkbox name="rms"> 



<h3>Obs types:</h3> 

<input type=radio name=obstype value="raob_qc" checked>raob_qc 

<input type=radio name=obstype value="sfc_lnd">sfc_lnd 

<input type=radio name=obstype value="sfc_ship_met_qc">sfc_ship_met_qc 

<h3>Graph Type:</h3> 

<!-<input type=radio name=graph value="scatter plot">scatter_plot— > 
<input type=radio name=graph value="time series" checked>time_series 

<h3>Period: Enter the beginning and ending DIG: e.g., 199803 18 12<7h3> 
<input type=text name=beginning maxlength=10 value-' 199803 1812"> 
<input type=text name=ending maxlength= 1 0 value=" 1 9980330 1 2"> 



<input type=submit> 

<input type=reset value="Cancel"> 

</form> 
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</bcxiy> 

<;/htnil> 



4. Noraps_conus.html 



<html> 

<!-- 

Author: Susie Pace 

Date: 25 March 1 998 

File URL: model_reports/mveri^noraps_conus.html 
— > 



<headxtitle>NORAPS_CONUS Verification Display </titlex/head> 

<body bgcolor=”#191970" TEXT="#00FF7F" LINK="#00FF7F" VLINK=”#CCCC66" 
ALINK=”#FF0000"> 

<foim method=GET action=”http://devul/cgi-bin/space3.pr’> 

<hlxcenter>Make your selections to see the statistics. </center></hl> 

<h3>Model: </h3> 

<input type=radio name=model value=”noraps_conus" checked> noraps_conus 
<h3>Geometry:</h3> 

<input type=radio name=geometry value="conus_nestl_appr' checked> conus_nestl_appl 
<h3>Parameters : </h3> 

<input type=radio name=parameter value="air_temp’*> air_temp 
<input type=radio name==parameter value="geop_ht” checked> geop__ht 
<input type=radio name=parameter value="pres”> pres 
<input type=radio name=parameter value="wnd_spd"> wnd_spd 

<h3>Taus:<7h3> 

<input type=radio name=tau value="0"> 0 
<input type=radio name=tau value=” 1 T> 1 2 
<input type=radio name=tau value=”24" checked> 24 
<input type=radio name=tau value=”36"> 36 
<input type=radio name=tau value=''48"> 48 

<h3>Levels: </h3> 

<input type=radio name=level value="0">0 
<input type=radio name=level value="2">2 
<input type=radio name=level value=" 1 9.5">1 9.5 
<input type=radio name=level value=" 1 000"> 1 000 
<input type=radio name=level value="925”>925 
<input type=radio name=level value="850">850 
<input type=radio name=level value='700">700 
<input type=radio name=level value='’500" checked>500 
<input type=radio name=level value="400">400 
<input type=radio name=level value=’'300'*>300 
<input type=radio name=level value="250">250 
<input type=radio name=level value="200”>200 
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<input type=radio name=level value="l50">150 
<input type=radio name=level value="100">100 

<! — <h3>Statistics :</H3> 
bias: <input type=checkbox name="bias"> 

stdev: <input type=checkbox name="stdev"> 

rms: <input type=checkbox name="rms"> 



<h3>0bs types:</h3> 

<input type=radio name=obstype value="raob_qc” checked>raob_qc 

<input type=radio name=obstype value="sfc_lnd">sfc_lnd 

<input type=radio name=obstype value="sfc_ship_met_qc">sfc_ship_met_qc 

<h3>Graph Type:</h3> 

<!— <input type=radio name=graph value=”scatter plot">scatter_plot-> 
<input type=radio name=graph value="time series" checked>time_series 

<h3>Period: Enter the beginning and ending DIG: e.g., 1998031 81 2</h3> 
<input type=text name=beginning maxlength=10 value="l 99803 181 2"> 
<input type=text name=ending maxlength=10 value-' 1 9980330 12"> 

<input type=submit> 

<input type=reset value="Cancel"> 

</form> 

</body> 

</html> 



5. Noraps_europe.html 



<html> 

<!- 

Author: Susie Pace 

Date: 25 March 1998 

File URL: model_reports/mveri£^noraps_europe.html 
— > 



<head><title>NORAPS_EUROPE Verification Display </titlex/head> 

<body bgcolor="#191970" TEXT="#00FF7F" LINK="#00FF7F" VLINK="#CCCC66" 
ALINK="#FF0000"> 

<fonn method=GET action="http://devul/cgi-bin/space3.pl"> 

<hl><Center>Make your selections to see the statistics.</center></hl> 

<h3>Model: </h3> 

<input type=radio name=model value="noraps_europe" checked> noraps_europe 
<h3>Geometry:</h3> 

<input type=radio name=geometry value="europe_nestl_appl" checked>europe_nestl_appl 
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<h3>P arameters: </h3> 

<input type=radio name=parameter value="air_temp"> air_temp 
<input type=radio name=parameter value=”geop_ht" checked> geop_ht 
<input type=radio name=paraineter value="pres"> pres 
<input type=radio name=parameter value="wnd_spd"> wnd__spd 

<h3>Taus:</h3> 

<input type=radio naine=tau value="0"> 0 
<input type=radio name=tau value=" 1 2"> 1 2 
<input type=radio name=tau value="24" checked> 24 
<input type=radio name=tau value="36"> 36 
<input type=radio name=tau value=”48"> 48 

<h3>Levels: </h3> 

<input type=radio name=level value="0">0 
<input type=radio name=level value="2">2 
<input type=radio name=level value=" 19.5 ">19.5 
<input type=radio name=level value=" 1 000 "> 1 000 
<input type=radio name=level value="925">925 
<input type=radio name=level value="850">850 
<input type=radio name=level value="700">700 
<input type=radio name=level value="500" checked>500 
<input type=radio name=level value="400">400 
<input type=radio name=level value="300">300 
<input type=Tadio name=level value="250">250 
<input type=radio name=level value="200">200 
<input type=radio name=level value=" 1 50">1 50 
<input type=radio name=level value- '100 ">100 

<!— <h3>Statistics:</H3> 
bias: <input type=checkbox nanie="bias"> 

stdev: <input type=checkbox name="stdev"> 

rms: <input type=checkbox name="rms"> 



<h3>Obs types:</h3> 

<input type=radio name=obstype value="raob_qc" checked>raob_qc 

<input type=radio name=obstype value="sfc_lnd">sfc_lnd 

<input type=radio name=obstype value="sfc_ship_met_qc">sfc_ship_met_qc 

<h3>Graph Type:</h3> 

<!— <input type=radio name=graph value="scatter plot">scatter_plot— > 
<input type=radio name=graph value="time series" checked>time__series 

<h3>Period: Enter the beginning and ending DIG: e.g., 1 99803 1 8 1 2</h3> 
<input type=text name=beginning maxlength=l 0 value-’ 1 99803 1 8 1 2"> 
<input type=text name=ending maxlength=l 0 value=" 1 9980330 1 2"> 

<input type=submit> 

<input type=reset value="Cancel"> 

</form> 
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</body> 

</htrol> 



6. Noraps_ind_ocn.html 



<html> 

<!- 

Author: Susie Pace 

Date: 25 March 1998 

File URL: model_reports/mverif^noraps__ind_ocn.html 



<head><title>NORAPS_IND_OCN Verification Display </title></head> 

<body bgcolor="#l 91 970" TEXT="#00FF7F" LINK="#00FF7F" VLINK="#CCCC66" 
ALINK="#FF0000"> 

<form method=GET action="http://devul/cgi-bin/space3.pr'> 

<hl><center>Make your selections to see the statistics.</center></hl> 

<h3>Model: </h3> 

<input type=radio name=model value="noraps_ind_ocn*' checked>noraps_ind_ocn 
<h3>Geometry:</h3> 

<input type=radio name=geometiy value="ind_ocn_nestl_appl" checked> 
ind_ocn_nestl _appl 

<h3>Parameters:</h3> 

<input type=radio name=parameter value="air_temp"> air_temp 
<input type=radio name=parameter value="geop_ht" checked> geop_ht 
<input type=radio name=parameter value="pres"> pres 
<input type=radio name=parameter value="wnd_spd"> wnd_spd 



<h3>Taus:</h3> 

<input type=radio name=tau value="0"> 0 
<input type=radio name=tau value=" 1 2"> 1 2 
<input type=radio name=tau value="24" checked> 24 
<input type=radio name=tau value="36"> 36 
<input type=radio name=tau value=’'48"> 48 

<h3>Levels: </h3> 

<input type=radio name=level value="0">0 
<input type=radio name=level value="2">2 
<input type=radio name=level value=" 1 9.5"> 1 9.5 
<input type=radio name=level value=" 1 000"> 1 000 
<input type=radio name=level value="925">925 
<input type=radio name=level value="850">850 
<input type=radio name=level value='700">700 
<input type=radio name=level value="500" checked>500 
<input type=radio name=level value="400">400 
<input type=radio name=level value="300">300 
<input type=radio name=level value="250">250 
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<input lype=radio name=level value="200”>200 
<input type=radio name=level value=" 150">1 50 
<input type=radio name=level value="100">100 

<!-<h3>Statislics:</H3> 
bias: <input type=checkbox name="bias"> 

stdev: <input type=checkbox name="stdev"> 

nns: <input type=checkbox naine="nns''> 

— > 



<h3>Obs types:</h3> 

<input type=radio name==obstype value="raob_qc" checked>raob_qc 

<input type=radio name=obstype value="sfc_lnd">sfc_lnd 

<input type=radio name=obstype value=”sfc_ship_met_qc">sfc_ship_met_qc 

<h3>Graph Type:</h3> 

<!-<input type=radio name=graph value=”scatter plot">scatter_j5lot— > 

<input type=radio name=graph value=''time series" checked>time_series 

<h3>Period: Enter the beginning and ending DIG: e.g., 1998031 81 2</h3> 
<input type=text name=beginning maxiength=l 0 value=" 1 99803 1 8 1 2"> 
<input type=text name=ending maxlength=10 value=" 199803 301 2 "> 

<input type=submit> 

<input type=reset value="Cancel"> 

</form> 

</body> 

</html> 

7. Coamps_europe.html 

<html> 

<!- 

Author: Susie Pace 

Date: 25 March 1 998 

File URL: model_reports/mverifi^coamps_europe.html 
— > 



<head><title>COAMPS_EUROPE Verification Display </titlex/head> 

<body bgcolor="#191970" TEXT="#00FF7F" LINK="#00FF7F" VLINK="#CCCC66" 
ALINK="#FF0000"> 

<form method=GET action=''http://devul/cgi-bin/space3.pl"> 

<hlxcenter>Make your selections to see the statistics.</center></hl> 

<h3>Model: </h3> 

<input type=radio name=model value="coamps_europe" checked> coamps_europe 



<h3>Geometiy : </h3> 

<input type=radio name=geometry value="europe_nest2_appl2" checked> 
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europe_nest2_appl2 

<input type=radio name=geometry value="europe_nest3_appl3">europe_nest3_appl3 
<h3>Paraineters:</h3> 

<input type=radio name=parameter value="air_temp"> air_temp 
<input type=radio name=parameter value="geop_ht" checked> geop_ht 
<input type=radio name=parameter value="pres"> pres 
<input type=radio name=parameter value=”wnd_spd"> wnd_spd 

<h3>Taus;<;/h3> 

<input type=radio name=tau value="0"> 0 
<input type=radio name=tau value=" 1 T> 1 2 
<input type=radio name=tau value="24" checked> 24 
<input type=radio name=tau value="36"> 36 
<input type=radio name=tau value="48"> 48 

<h3>Levels: </h3> 

<input type=radio name=level value="0">0 
<input type=radio name=level value="2">2 
<input type=radio name=level value=" 1 9.5">1 9.5 
<input type=radio name=level value- '1000”> 1000 
<input type=radio name=level value="925">925 
<input type=radio name=level value="850">850 
<input type=radio name=level value="700">700 
<input type=radio name=level value="500" checked>500 
<input type=radio name=level value="400">400 
<input type=radio name=level value="300">300 
<input type=radio name=level value="250">250 
<input type=radio name=level value="200">200 
<input type=radio name=level value=" 1 50"> 1 50 
<input type=radio name=level value-' 100"> 100 



<!— <h3>Statistics:</H3> 
bias: <input type=checkbox name="bias"> 

stdev: <input type=checkbox name="stdev"> 

rms: <input type=checkbox name="nns"> 

— > 



<h3>Obs types:</h3> 

<input type=radio name=obstype value="raob_qc” checked>raob_qc 

<input type=radio name=obstype value="sfc_lnd">sfc_lnd 

<input type=radio name=obstype value="sfc_ship_met_qc">sfc_ship_met_qc 

<h3>Graph Type:</h3> 

<!— <input type=radio name=graph value=" scatter plot">scatter_plot— > 
<input type=radio name=graph value="time series" checked>time_series 

<h3>Period: Enter the beginning and ending DTG: e.g., 1 99803 181 2</h3> 
<input type=text name=beginning maxlength=l 0 value=" 1 99803 1 8 1 2"> 
<input type=text name=ending maxlength= 1 0 value=" 1 99803 30 1 2 "> 

<input type=submit> 
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<input type=reset value=”Cancer'> 



</form> 

<A)ody> 

</html> 



8. Coamps_southwest_asia.html 



<html> 

<!- 

Author: Susie Pace 

Date: 25 March 1998 

File URL: model_reports/mverif/coamps_swasia.html 



<headxiitle>COAMPS__SOUTHWEST_ASIA Verification Display </title></head> 
<body bgcolor="#191970" TEXT="#00FF7F” LINK=”#00FF7F" VLINK="#CCCC66” 
ALINK="#FF0000"> 

<form method=GET action="http://devul/cgi-bin/space3.pl"> 

<hlxcenter>Make your selections to see the statistics. </centerxyhl> 

<h3>Model: </h3> 

<input type=radio name=model value="coamps_sw_asia” checked> 
coamps_southwest_asia 

<h3>Geometiy : </h3> 

<input type=radio name=geometiy value="southwest_asia_nest2_appl" checked> 
southwest_asia_nest2_appl 

<input type=radio name=geometry value="southwest_asia_nest3_appr> 
southwest_asia_nest3_appl 

<h3>Parameters:</h3> 

<input type=radio name=parameter value="air_temp"> air_temp 
<input type=radio name=parameter value="geop_ht" checked> geop_ht 
<input type=radio name==parameter value="pres"> pres 
<input type=radio name=parameter value="wnd_spd"> wnd_spd 

<h3>Taus:</h3> 

<input type=radio name=tau value="0"> 0 
<input type=radio name=tau value="12"> 12 
<input type=radio name=tau value="24" checked> 24 
<input type=radio name=tau value="36"> 36 
<input type=radio name=tau value='’48"> 48 

<h3>Levels: </h3> 

<input type=radio name=level value="0">0 
<input type=radio name=level value="2">2 
<input type=radio name=level value=" 1 9.5">1 9.5 
<input type=radio name=level value=" 1 000”>1 000 
<input type=radio name=level value="925">925 
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<input type=Tadio name=level value="850">850 
<input type=radio name=level value='700">700 
<input type=radio name=level value="500" checked>500 
<input type=radio name=level value="400">400 
<input type=radio name=level value="300">300 
<input type=radio name=level value="250">250 
<input type=radio name=level value="200">200 
<input type=radio name=level value=" 1 50">1 50 
<input type=radio name=level value="100">100 

<!— <h3>Statistics:</H3> 
bias; <input type=checkbox name="bias"> 

stdev: <input type=checkbox name="stdev"> 

nns: <input type=checkbox name="rTns"> 

— > 



<h3>Obs types:</h3> 

<input type=radio name=obstype value=’'raob_qc“ checked>raob_qc 

<input type=radio name=obstype value="sfc_lnd">sfc_lnd 

<input type=radio name=obstype value="sfc_ship_met__qc">sfc__ship_met_qc 

<h3>Graph Type:</h3> 

<!~<input type=radio name=graph value="scatter plot'’>scatter_plot— > 
<input type=radio name=graph value="lime series" checked>time_series 

<h3>Period: Enter the beginning and ending DIG: e.g., 199803 181 2</h3> 
<input type=text name=beginning maxlength= 1 0 value=" 1 99803 1 8 1 2"> 
<input type=text name=ending maxlength=10 value-' 1 9980330 12"> 

<input type=submit> 

<input type=reset value="Cancel"> 

</form> 

</body> 

</html> 



9. Warn GIobaI.html 



<html> 

<!- 

Author: Susie Pace 

Date: 25 March 1998 

File URL: model_reports/mverif/wam.html 
--> 



<headxtitle>WAM_GLOBAL Verification Display </titlex/head> 

<body bgcolor="#191970" TEXT="#00FF7F" LINK="#00FF7F" VLINK="#CCCC66" 
ALINK="#FF0000"> 

<fomi method=GET action="http://devul/cgi-bin/space3.pr> 

<hlxcenter>Make your selections to see the statistics.</center></hl> 
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<h3>Model: </h3> 

<input type=radio name=model value=”wam_global" checked> wam_global 
<h3>Geometiy : </h3> 

<input type=radio name=geometiy value="global_360xl81 " checked> 
global J60xl 81 

<h3>Parameters:</h3> 

<input type=radio name=parameter value="sig_wav_ht" checked> sig_wav_ht 
<input type=radio name=parameter value="peak_wav_per"> peak_wav_per 

<h3>Taus:</h3> 

<input type=radio name=tau value="0"> 0 
<input type=radio name=tau value=” 1 2"> 1 2 
<input type=radio name=tau value="24" checked> 24 
<input type=radio name=tau value="36"> 36 
<input type=radio name=tau value="48"> 48 
<input type=radio name=tau value="60"> 60 
<input type=radio name=tau value='72"> 72 
<input type=radio name=tau value="84"> 84 
<input type=radio name=tau value="96"> 96 
<input type=Tadio name=tau value=" 1 08"> 1 08 
<input type=Tadio name=tau value=” 1 20"> 1 20 

<h3>Levels: </h3> 

<input type=radio name=level value="0" checked>0 



<! — <h3>Statislics:</H3> 
bias: <input type=checkbox name="bias"> 

stdev: <input type=checkbox name=”stdev"> 

rms: <input type=checkbox name="nns"> 



<h3>Obs types:</h3> 

<input type=radio name=obstype value="sfc_ship" checked>sfc_ship 
<!— <input type=radio name=obstype value="alty">alty— > 



<h3>Graph Type:</h3> 

<!— <input type=radio name=graph value=" scatter plot">scatter_j)lot— > 
<input type=Tadio name=graph value="time series" checked>time_series 

<h3>Period: Enter the beginning and ending DIG: e.g., 199803 181 2</h3> 
<input type=text name=beginning maxlength= 1 0 value=" 1 99803 1 8 1 2"> 
<input type=text name=ending maxlength=10 value-' 199803301 2 "> 

<input type=submit> 

<input type=reset value="Cancel"> 

<Jfonn> 

</body> 

</html> 
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10. Procform.pl 



#!/usr/local/bin/perl 

# space3.pl - Try to batch in ninjob to retreive the data 

# and run the DDL programs to create the graphs 

#Use the form libraiy 
require "space2.pl"; 

########## STEP 1 : Get and decode the input from the form 

#get the data from the form 
&ReadParse(*input); 



STEP 2: Process the information from the form 



Sprocid = $$; 

'export procid'; 

#split $query_string into name=value pairs 
local(*FormData) = if #make an alias for the arg 
$query_string=$ENV{'QUERY_STRING’}; 

foreach $name_value (split('&’, $query_string)) { 

#translate any plus signs in the pair string into spaces: 
$name_value — tr/+//; 

#split the name=value pair into a separate name and value: 
(Sname, Svalue) = split $name_value); 

#translate escaped hex numbers back to 8-bit char: 

Sname — s/%(..)/pack("C", hex($l))/eg; 

Svalue — s/%(..)/pack("C", hex(Sl))/eg; 

#convert the names and values into named, assigned var; 
if (defined(SFormData{Sname})) { 

SFormData{Sname} .= "\0Svalue"; 

} else { 

SFormData{ Sname} = Svalue; 

} 

if (Sname eq ’model') { 

Smodel = Svalue; 

'export model'; 

#exit 0 if (Smodel eq ") 

} elsif (Sname eq 'geometry') { 

Sgeometry = Svalue; 

' export geometry' ; 

#exit 0 if (Sgeometry eq ") 
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} elsif (Sname eq 'parameter') { 
Sparameter = Svalue; 

' export parameter' ; 

#exit 0 if (Sparameter eq ") 

} elsif (Sname eq 'tau') { 

Stau = Svalue; 

'export tau'; 

#exit 0 if (Stau eq ") 

} elsif (Sname eq 'level') { 
Slevel = Svalue; 

'export level'; 

#exit 0 if (Slevel eq ") 

} elsif (Sname eq 'stats') { 

Sstats = Svalue; 

#' export stats'; 

} elsif (Sname eq 'obstype') { 
Sobstype = Svalue; 

'export obstype'; 

#exit 0 if (Sobstype eq ") 

} elsif (Sname eq 'graph') { 
Sgraph = Svalue; 

' export graph' ; 

#exit 0 if (Sgraph eq ") 

} elsif (Sname eq 'beginning') { 
Sbeginning = Svalue; 

'export beginning' ; 

#exit 0 if (Sbeginning eq ") 

} elsif (Sname eq 'ending') { 
Sending = Svalue; 

'export ending'; 

#exit 0 if (Sending eq ") 

} else { 

print "no matching env var"; 

} 



#process the request by runjob 

'/a/ops/bin/runjob -h div60-3 -u pacek -t sun -d /home/pubs43/tmp -j get_data.ksh -e "pid=Sprocid 
MODEL=S model BDTG=Sbeginnmg EDTG=Sending PARM_NAME=Sparameter LVL_l=Slevel 
FCSTPER=Stau OBSTYPE=Sobstype GEOM_NAME=Sgeometry" > Sprocid.out 2>&1' ; 



STEP 3: Reply by outputting a new 



} 
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#$Image = "http://l 52.80. 1 3.201/'-pacek/gify$procid.gif'; 
SImage = "/home/pubs43/tmp/$procid.gif'; 



#!!!!! This portion worked for WAM ill!!# 

# 

# Output the MIME- type header, followed by two newlines: 
print "Content-type: image/gif\n\n"; 

# 

# loop until the gif file is transferred to the web server 

# when the gif file arrives, display on the web browser 

# 

while (I (-e SItnage)) { 
sleep 5; 

} 

open (IMAGE, SImage); 
print <IMAGE>; 
close (IMAGE); 

'nn SImage'; 

#print "<HTML><HEAD><TITLE> Model Verification Reply </riTLE></HEAD>\n"; 

#print "</BODY></HTML>\n"; 

# 

#!!!!! Down to here !!!!!# 

#Trying to handle the timed out error 
#if (-e SImage) { 

# print "Content-type: image/gif\n\n"; 

# open (IMAGE, SImage); 

# print <IMAGE>; 

# close (IMAGE); 

#this causes the runjob to be batched in multiple time and PBD change, 

#end up in an infinite loop 
#} else { 

# print "refi-esh: 1 0; \n"; 

# print "Content-type: text/html\n\n"; 

# print "<HTML><HEAD><TITLE>Model Statistics Display<VriTLE></HEAD>\n"; 

# print "<BODYBGCOLOR=NAVYTEXT=WHITE>"; 

# print "<H1 ><CENTER>Please be patient, the image is being created</CENTER></Hl >\n" 

# print "Processing is done when a graph is displayed\n"; 



END OF space3.pl ########## 



11. Get_data.ksh 



#!^in/sh 
/bin/ksh «’EOT' 

. /a/ops/isis/db_init/isis_init_ofs.ksh 
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#set -X 
cd/d/tmp 

export DB_DIR=/ci/model -stats/ 
cd $DB_DIR 

export FILENAME=$pid 
export TMPFN=$pid"tmp" 

SLASH=_ 

if [[ SMODEL = nogaps ]] 
then 

MODEL="$MODEL$SLASH$GEOM_NAME" 

fi 

print SMODEL $GEOM_NAME SFCSTPER $PARM_NAME SOBSTYPE $LVL_1 SBDTG SEDTG 

empcmd $DB_DIR/stat_db 'select bypass_lock verif_date,sample_size,pann_name, 
unit_name,geom_name,lvl_type,level_l convert to decimal(8,2), 
tau,stat_type,stat_value convert to decimal(8,2),verif_soxirce, 
obs_type from SMODEL where verif_date range SBDTG to SEDTG and 
pann__name="SPARM_NAME" and geom_name="SGEOM_NAME" and tau="SFCSTPER" 
and obs_type="SOBSTYPE" and level_l=”SLVL_r into "STMPFN";' 

awkV^[0-9]/\ 

(printf ("%- 1 2s%-5d%-20s%- 1 5s%-3 0s%- 1 5s%-8.2P/o-5d%- 1 5s%-8.2f%- 1 0s%-25s\n", \ 
S1,S2,S3,S4,S5,S6,S7,S8,S9,S10,S1 1,S12)}' STMPFN > SFILENAME 

#print "sourcing the IDL setup script" 

. /h/idl/idl_5/bin/idl_setup.ksh 

#print SPARM_NAME SLVL_1 SFCSTPER SOBSTYPE SGEOM_NAME 
idl -rt=plot_data 

if[[ -fSFILENAME.gif]] 
then 

ftpbatch -h devul -s "cd /home/pubs43/tmp/; put SFILENAME.gif' 
fi 

rm Spid* 

exit 0 
EOT 
exit 0 
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APPENDIX C 



A. STATTEST.F90 

program stattest 
C 

C START PROLOGUE 

C 

C sees IDENTIFICATION: @(#)stattest.I90 1.1 04/24/98 
C 

C CONFIGURATION IDENTffICATION: NONE 
C 

C MODULE NAME: stattest 
C 

C DESCRIPTION: Program to test the stat lib 
C 

C COPYRIGHT: (c) 1998 FLENUMMETOCCEN 

C US. GOVERNMENT DOMAIN 

C ALL RIGHTS RESERVED 

C 

C CONTRACT NUMBER AND TITLE: N/A 
C 

C REFERENCES: NONE 
C 

C CLASSIFICATION: Unclassified 
C 

C RESTRICTIONS; NONE 
C 

C COMPUTER/OPERATING SYSTEM 
C DEPENDENCIES: CrayUNICOS 

C 

C LIBRARIES OF RESIDENCE: /a/opsA>in 
C 

C USAGE; 

C stattest 
C 

C PARAMETERS; 

C Name Type Usage Description 

C 

c 

C PARAMETERS; 

C Name Type Usage Description 

C 

c 

C COMMON BLOCKS; N/A 
C 

C FILES; N/A 
C 

C DATA BASES; N/A 
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c 

C COMPILER DEPENDENCIES: 190 
C 

C COMPILE OPTIONS: -f fixed -c 
C 

C MAKEFILE: Located at /a/ops/app/mverif/test/maketest 
C UNICOS make 
C 

C RECORD OF CHANGES: 

C 

C «CHANGE NOTICE» Version 1.1 (29 Apr 1998) — Kyongsuk Pace 
C Initial submission 
C 

C END PROLOGUE 

C 

implicit none 

integer size 
parameter(size= 1 0) 

real :: stat 
real :: arrl(size) 
real :: arr2(size) 

character(4) :: geom 

geom(l ;4) = 'NONE' 

arr2 = (Z288.6, 304.8, 301.4, 293.2, 293.2, 

2 297.6,295.8,291.0,285.2,287.8/) 

arrl = (Z289.29, 302.26, 302.14, 294.92, 294.92, 

2 296.25, 295.37, 292.32, 285. 1 5, 288.69/) 

CALL COMPUTE_BIAS(arrl ,arr2,size,geom,stat) 
write (0, '("bias = ", 18.2)') stat 

CALL COMPUTE_STD(arrl ,arr2,size,geom,stat) 
write (0, '("std = ", 18.2)') stat 

CALL COMPUTE_RMS(arr 1 ,arr2,size,geom,stat) 
write (0, '("rms = ", 18.2)') stat 

stop 'Normal End' 

end 



B. STATISTICS BY STATISTICS LIBRARY SUBROUTINES 



Script started on Wed Feb 11 1 7:2 1 ;26 1 998 
[Zdev/ttyp017] 

s{j91 }/home/pacek/mveriEtest>tstattest 
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bias = 0.27 

std= 1.31 
nns= 1.34 
STOP Normal End 

STOP executed at line 28 in Fortran routine 'STATTEST' 

CPU: 0.010s, Wallclock:0.011s, 10.6% of 8-CPU Machine 

Memory HWM; 1595052, Stack HWM: 803866, Stack segment expansions: 0 

(j91 }yhome/pacek/mverif/test 

Script finished on Wed Feb 11 17:21:34 1998 



C. STATISTICS BY MICROSOFT EXCEL 



array 1 


! array 2 diff 


288.6 


289.29 


304.8 


302.26 


301.4 


302.14 


293.2 


294.92 


293.2 


294.92 


297.6 


296.25 


295.8 


295.37 


291 


292.32 


285.2 


285.15 


287.8 


288.69 


sum 


-2.71 


bias 


-0.27 



diff sqr 


-0.69 


0.4761 


2.54 


6.4516 


-0.74 


0.5476 


-1.72 


2.9584 


-1.72 


2.9584 


1.35 


1.8225 


0.43 


0.1849 


-1.32 


1.7424 


0.05 


0.0025 


-0.89 


0.7921 


17.9365 




sum/count 



count 

stdev 

rms 



1.3115 sqrt(sum of diff sqr/count - bias*bias) 

1 .3392 sqrt(sum of diff sqr/count) 
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